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VOORWOORD 


Dit boek is bedoeld voor studenten іп de studierichtingen есо- 
nomie, bedrijfskunde en bestuurskunde, Behalve voor studenten 
aan universitaire opleidingen is het boek ook zeer geschikt voor 
studenten aan een Н.Е.А.0. ої H.T.S. 


Bij de behandeling van de stof is gekozen voor een intuïtieve, 
niet-wiskundige benadering; Aan de economische betekenis van 

de bij lineaire programmering gebruikte begrippen en aan de be- 
tekenis van deze begrippen voor de besluitvorming wordt ruime 
aandacht gegeven. In de praktijk is gebleken, dat deze benade- 
ring ook bij het onderwijs aan studenten, die in hun eindexamen- 
pakket van de middelbare school geen wiskunde kozen, goede re- 
sultaten oplevert. 

Voorkennis op het gebied van de lineaire algebra is bij deze 
benadering niet vereist. 


In hoofdstuk 1 wordt de modelcyclus behandeld. De modelcyclus 
wordt geïllustreerd met een voorbeeld op het gebied van de voor- 
raadbeheersing. Eventueel kan men hoofdstuk 1 zonder verlies van 
continuiteit overslaan. 


In de hoofdstukken 2, 3, 4 еп 5 wordt de methodiek van de line- 
aire programmering aan de hand van voorbeelden uiteengezet. Deze 
hoofdstukken vormen de kern van het boek en zouden in iedere 


cursus, waarin dit boek wordt gebruikt, grondig moeten worden 
behandeld. 


In de hoofdstukken 6, 7 en 8 komen onderwerpen als het duale 
probleem, de gevoeligheidsanalyse (in hoofdstuk 3 al op elemen- 
taire wijze behandeld) en het probleem van het hanteren van 
meervoudige doelstellingen aan de orde. 


In hoofdstuk 9 еп 10 worden twee toepassingen van lineaire pro- 
grammering behandeld. Daarbij wordt gebruik gemaakt уап ІМОРТІ. 


LINEAIRE PROGRAMMERING 
ALS HULPMIDDEL BIJ 
DE BESLUIT VORMING 


ІМОРТІ is een interactief programma voor LP-problemen уап be- 
perkte omvang, dat door het Instituut voor Ontwikkeling van het 
Wiskunde Onderwijs (І.0.М.0.) te Utrecht speciaal voor onder- 
wijsdoeleinden is vervaardigd. De benodigde programmatuur wordt 
door het І.0.М.0. gratis ter beschikking gesteld aan universi- 
teiten, Н.Е.А.0.'8 еп H.T.S.'en. 


Het verwerven van vaardigheid in het formuleren van LP-modellen, 
in het oplossen van het model met behulp van een programma als 
IMOPTI en in het interpreteren van de output, beschouwt de 
schrijver als de belangrijkste leerdoelstelling bij een cursus 
waarbij dit boek zou kunnen worden gebruikt. 


Gaarne spreek ik mijn dank uit aan mijn collega Drs. S.K.Th. 
Boersma, die mij bij het schrijven van dit boek enkele waarde- 
volle adviezen heeft gegeven en aan de studenten van de Inter- 
faculteit Bedrijfskunde van de Rijksuniversiteit te Groningen, 
die door het stellen van vragen tijdens de colleges indirect 
aan de vormgeving van dit boek hebben bijgedragen. 


Roden, april 1979 S.W. Douma 


VOORWOORD BĲ DE TWEEDE DRUK 


Het is verheugend te kunnen constateren, dat dit boek zijn weg 
heeft gevonden naar verschillende opleidingen, zodat binnen 
korte tijd een tweede druk noodzakelijk werd. Deze tweede druk 
15 nagenoeg gelijk aan de eerste druk; twee vraagstukken (4-5 
en 5-2) zijn gewijzigd en enkele nieuwe vraagstukken zijn toe- 
gevoegd (2-7, 5-3 en 8-1). Tenslotte zijn de drukfouten, die de 
eerste druk nog ontsierden, hersteld. 


Roden, januari 1982 S.W. Douma 


VOORWOORD BĲ DE DERDE DRUK 


Omdat zowel docenten als studenten gunstig hebben gereageerd op 
de eerste en tweede druk van dit boek, is de tekst van de derde 
druk bijna geheel gelijk aan die van de tweede druk. Wel deden 
verschillende gebruikers het verzoek een groter aantal vraagstuk- 
ken en uitwerkingen van vraagstukken in het boek op te nemen. Van 
alle vraagstukken, die in de tweede druk voorkwamen (27 stuks) 
zijn daarom nu ook de uitwerkingen opgenomen. Daarnaast zijn 15 
nieuwe vraagstukken toegevoegd. 


Sedert het verschijnen van de eerste druk in 1979 heeft de соп- 
putertechnologie een enorme ontwikkeling doorgemaakt. Veel onder- 
wijsinstellingen beschikken thans over eigen microcomputers. 
Voorts is thans een LP-pakket (OPTIMIZER) in Nederland verkrijg- 
baar, dat zeer geschikt is voor onderwijsdoeleinden en dat gebrui- 
kersvriendelijk is. Dit pakket is verkrijgbaar bij SOFTKEY B.V. 

te Deventer. Met behulp van OPTIMIZER kunnen studenten ervaring 
opdoen in het formuleren van LP-modellen en in het oplossen daar- 
van op een microcomputer. In hoofdstuk 9 is naast de output van 
IMOPTI nu ook de output van OPTIMIZER opgenomen. 


De schrijver spreekt de hoop uit, dat door bovengenoemde aanvul- 
lingen de gebruikswaarde van het boek nog is toegenomen. 


Roden, april 1984 S.W. Douma 


l. DE AARD VAN DE KWANTITATIEVE 
BESLUITVORMINGSTECHNIEKEN 


1.1. GESCHIEDENIS 


De geschiedenis van de kwantitatieve besluitvormingstechnieken 
(Operations Research (Am.), Operational Research (Eng.), Opera- 
tionele Research (Ned.), Management Science) gaat terug tot het 
begin van de Tweede Wereldoorlog 1). Gedurende de tweede wereld- 
oorlog werden voor het eerst kwantitatieve besluitvormingstech- 
nieken op grote schaal en met succes toegepast op militaire 
problemen van operationele aard. Na de oorlog begonnen de grote- 
re organisaties eveneens de nieuwe kwantitatieve technieken toe 
te passen. 


1.2. TOEPASBAARHEID 


Men kan onderscheid maken tussen operationele en strategische 
beslissingen. Operationele en strategische beslissingen ver- 
schillen op meerdere punten van elkaar zoals is te zien in on- 
derstaande tabel: 


Operationele Strategische 
beslissingen beslissingen 
Tijdsduur die bij probleemstel- 
ling in beschouwing wordt 
genomen kort lang 
Is een eenmaal genomen beslis- 
sing gemakkelijk reversibel? ja nee 
Draagt de beslissing een repe- 
terend karakter ja nee 
Deel van de organisatie waarop 
de beslissing betrekking heeft klein groot 


Doelstellingen gegeven ja nee 
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Voorbeelden van beslissingen уап operationele aard worden bijv. 
gevonden in de korte-termijn planning in de fysieke distributie. 
Voorbeelden van beslissingen van strategische aard zijn: de be- 
slissing van Volkswagen om een eigen fabriek te bouwen in de 
Verenigde Staten, het vaststellen уап de portfolio van produkt - 
markt - combinaties op lange termijn bij de corporate planning, 
het vaststellen van een structuurplan door een gemeente. 

Tot dusverre zijn de toepassingen op problemen van operationele 
aard het meest talrijk. Toepassingen op strategische problemen 
zijn relatief aanzienlijk minder talrijk; op dit gebied liggen 
nog aanzienlijke mogelijkheden. 


1.3. DOEL VAN HET BOEK 


Het doel van dit boek is de student inzicht te verschaffen in 
de mogelijkheden en de beperkingen van kwantitatieve besluit- 
vormingsmethoden. Daartoe is vanzelfsprekend een minimum aan 
kennis van deze technieken vereist. 


1.4. DE MODELCYCLUS 


1.4.1. DE MODELCYCLUS IN DE SYSTEEMLEER 


In de systeemleer 2) kent men de modelcyclus 3). Deze is sche- 
matisch weergegeven in figuur l.l 


empirisch 
concreet systeem 


- figuur 1.1.: De modelcyclus - 


De modelcyclus is op iedere empirische wetenschap van toepassing. 
In de modelcyclus kan men drie fasen onderscheiden, te weten de 
fase van de abstractie (ook wel aangeduid als de modelbouw) de 
fase van de deductie en de fase van de realisatie (in figuur 

2.1. aangeduid met de letters a, d еп т). Een nadere verfijning 
van deze drie fasen is: 


abstractie | formuleer het probleem 
2. construeer het model 
3 verzamel de gegevens 
4 test het model 
deductie 5 
6 
7 


realisatie 


zoek de oplossing 
test de oplossing 
implementeer de oplossing 


De voorgenoemde zeven fasen geven een indruk van de stadia die 
men bij het toepassen van kwantitatieve besluitvormingstechnie- 
ken in de praktijk zal moeten doorlopen. De gegeven volgorde 
suggereert tevens een tijdsvolgorde. Aan de hier gegeven tijds- 
volgorde zal men echter bepaald niet in alle gevallen de hand 
willen of kunnen houden. Zo zal men bijvoorbeeld de fasen 2 en 
3 in vele gevallen ook gelijktijdig of juist in omgekeerde 
volgorde doorlopen. 

Wij bespreken nu de verschillende fasen aan de hand van een 
voorbeeld. Wij kiezen als voorbeeld een situatie, waarin aan 
een onderzoeker (men kan hier ook lezen: een organisatie-advi- 
seur) wordt gevraagd behulpzaam te zijn bij een voorraadbeheer- 
singsprobleem bij een groothandel. 


1.4.2. HET FORMULEREN VAN HET PROBLEEM 


De eerste stap bij het formuleren van het probleem bestaat uit 
het maken van een inventarisatie van de verschillen tussen de 
gewenste toestand en de werkelijke toestand waarin het concrete 
systeem zich bevindt. Daarbij zal men een aantal in de organi- 
satie werkzame personen interviewen. Het is zeer goed denkbaar, 
dat verschillende personen verschillende visies hebben ten aan- 
zien van de vraag wat als de meest gewenste toestand van het 
systeem wordt beschouwd. 


In ons voorbeeld zou bijvoorbeeld de inkoper kunnen streven naar 
het plaatsen van grote bestellingen, teneinde op deze wijze de 
bestelkosten laag te houden. De verkoper zou kunnen streven naar 
een service-graad van 1002, d.w.z. naar zodanige voorraden dat 
iedere klantenorder steeds onmiddellijk uit voorraad kan worden 
geleverd. Tenslotte zou de financieringsman uit overwegingen 
van solvabiliteit en liquiditeit kunnen streven naar zo laag 
mogelijke investeringen in vaste activa (magazijnruimte) en 
vlottende activa (voorraden). 


De fase van de probleemformulering moet dus in elk geval uit- 

monden in een goede kwalitatieve beschrijving van de doelstel- | 

lingen van alle betrokkenen. In ons voorbeeld zou dit bijvoor- 

beeld kunnen zijn: 

doelstelling 1: streef naar zo laag mogelijke kosten, waarbij 
met kosten wordt bedoeld de som van bestelkosten 
en voorraadkosten. 
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doelstelling 2: Streef naar een zo hoog mogelijke servicegraad. 


In ons voorbeeld zou men de totale kosten (som van bestelkosten 
en voorraadkosten) kunnen aanduiden met uj en de servicegraad 
(dat is het gedeelte van de door klanten gevraagde hoeveelheden, 
dat direct uit voorraad kan worden geleverd) met u2. Men noemt 
dan uj en uy de uitgangsvariabelen (outputvariables). 


De tweede vraag, die men bij de probleemformulering moet beant- 
woorden is: ten aanzien van welke grootheden kan de organisatie 
zelf beslissingen nemen? Anders gezegd: wat zijn de beslissings- 
variabelen (controlled variables)? In ons voorbeeld zouden be- 
slissingsvariabelen kunnen zijn xjj = de bestelhoeveelheid van 
produkt i en хі? = het bestelniveau van produkt i. Tenslotte 
zal men zich moeten afvragen ten aanzien van welke grootheden 
die bij de probleemstelling een rol spelen, de organisatie zelf 
geen beslissingen kan nemen. Anders gezegd wat zijn de ingangs- 
variabelen (input variables, uncontrolled variables). In ons 
voorbeeld zouden uncontrolled variables kunnen zijn: yil = de 
levertijd van produkt i en yi2 = de vraag per periode naar pro- 
dukt 1. 

De fase van de probleemstelling moet dus tevens uitmonden in 
een uitputtende beschrijving van alle beslissings- en ingangs- 
variabelen. 


1.4.3. HET CONSTRUEREN VAN HET MODEL 


1.4.3.1. SOORTEN MODELLEN 


Er bestaan zeer veel verschillende soorten modellen. Wij geven 
in dit boek geen volledige beschrijving van de verschillen- 

de soorten modellen ") maar volstaan met enkele indelingen die 
voor het vervolg van dit boek van belang zijn. 


a. ICONISCHE, ANALOGE OF SYMBOLISCHE MODELLEN 

Iconische modellen zijn schaalmodellen, bijvoorbeeld een schaal- 
model van een te ontwerpen vliegtuig, waarvan bepaalde eigen- 
schappen in een windtunnel worden beproefd. 

Analoge modellen zijn modellen, waarbij een bepaalde verzame- 
ling van eigenschappen van het model wordt gebruikt als repre- 
sentatie van de verzameling van eigenschappen van het concrete 
systeem. Zo kan men bijvoorbeeld hydraulische modellen maken 
van economische systemen, waarbij vloeistofstromen in het model 
de geldstromen in het economische systeem representeren. Een 
ander bekend voorbeeld betreft een model, dat is ontworpen om 

de optimale plaats van een depot te vinden. Men neemt dan een 
horizontale tafel, plakt daarop een kaart waarop de geografische 
plaats van de klanten is aangegeven en boort overal waar een 
klant zich bevindt een gaatje in de tafel. Indien men nu door 
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al deze gaatjes een koordje haalt en aan de koordjes gewichten 
bevestigt, waarvan de grootte evenredig is met de door klanten 
afgenomen hoeveelheden en knoopt men de koordjes boven de tafel 
alle aan dezelfde ring dan zal de plaats van de ring de opti- 
male plaats van het depot aangeven (daarbij zijn dan nog de 
volgende veronderstellingen gemaakt: le er is geen wrijving, 
2e de transportkosten zijn gelijk aan afgelegde weg x hoeveel- 
heid en 3e een rechte lijn op de kaart is een goede benadering 
van de transportafstand). 

Symbolische modellen ook dikwijls aangeduid als mathematische 
modellen vormen het onderwerp van dit boek. Een stelsel van 
mathematische vergelijkingen representeert dan het concrete 
systeem. 


b. OPTIMALISATIEMODELLEN EN BESCHRIJVENDE MODELLEN 

Bij optimalisatiemodellen treft men steeds één of meer beslis- 
singsvariabelen aan. 

In een beschrijvend model zijn geen beslissingsvariabelen opge- 
nomen. Het door ons in de voorgaande paragraaf geformuleerde 
voorraadbeheersingsprobleem leidt tot een mathematisch model 
van de volgende vorm (wij laten van nu af aan gemakshalve de 
index i weg): 


чү - f (x), X>, У], y2) 
ua = Ғ2(хі, Xo, Уі» Y2) 


Dit model is dus een optimalisatiemodel. 


Voorbeelden van beschrijvende modellen zijn: 

- de modellen gebruikt bij het voorspellen van tijdreeksen (bij- 
voorbeeld de exponential smoothing-modellen) 

- een model voor de middellange termijnplanning bij ondernemin- 
gen dat werd opgesteld om de gevolgen van wijzigingen in te 
verkopen hoeveelheden, produktiecapaciteiten, de financiering 
van de onderneming enzovoorts snel te kunnen doorrekenen ? 


с. DETERMINISTISCHE VERSUS STOCHASTISCHE MODELLEN 

In deterministische modellen zijn alle variabelen van determi- 
nistische aard, in een stochastisch model komen één of meer 
stochastische variabelen voor. (Een stochastische variabele is 
een variabele, waarvan de waarde niet met zekerheid bekend is, 
doch onderworpen is aan een kansverdeling.) In ons voorbeeld 
kan men bijvoorbeeld y en у; beide als deterministisch opvat- 
ten (of dit een zinvolle representatie van het concrete systeem 
is moet blijken uit het waarnemen van het concrete systeem). 
Men spreekt dan van een deterministisch voorraadmodel. Dikwijls 
zal men echter situaties aantreffen waarbij men у, (= de vraag 
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per periode) als een stochastische variabele moet beschouwen. 
Hetzelfde geldt voor y) (= de levertijd). Zijn у» en/of у; sto- 
chastische variabelen dan spreekt men van een stochastisch 
voorraadmodel. 

Zijn yj en/of yy stochastisch, dan zullen іп het algemeen ook 
de uitgangsvariabelen u еп uy stochastisch zijn. Nu is men 
gewoonlijk niet zozeer geïnteresseerd in de waarden die uj en 
uy aannemen bij een bepaalde waarde van y; еп уз als wel in de 
kansverdelingen van uj en ua, Vaak zal men ook genoegen nemen 
met het leren kennen van bepaalde parameters van deze kansver- 
delingen, zoals bijvoorbeeld E(uj) en o?{uj} (respectievelijk 
de verwachtingswaarde en de variantie van uj). 


1.4.3.2. HET CONSTRUEREN VAN HET MODEL 


Het gaat er пи om, dat wij de vorm van de functies fj en f2 
specificeren. Vatten wij у) еп у» als deterministische variabe- 
len ор, dan kan men de vorm van de functies fj en Ё, als volgt 
specificeren: 


voorraadhoogte 


D K 
- fig. 2.2. Een deterministisch voorraadmodel - 


OF = xj = bestelniveau 


AD = xj = bestelhoeveelheid 
EC = y т levertijd 

Е , 
tan а = FH 727 de vraag per periode 


Wij beschouwen in totaal n perioden. Zij nu a=de bestelkosten 
per keer еп b - de opslagkosten van één eenheid produkt gedu- 
rende één periode. 


De servicegraad чи; --- 


АС - AD - CD 
Teneinde CD te berekenen merken wij op dat 
CD = BD - ВС. 


Hoek DHB is gelijk aan а; er geldt dus um = tan о 


> BD = HB.tana = ууу» 


то з yy 7 52 


Men kan dit als volgt interpreteren: 

CD = de grootte van de "пеепуегкоор", dat wil zeggen de hoeveel- 
heid gevraagde goederen, die men niet onmiddellijk uit voorraad 
kan leveren. Deze hoeveelheid is gelijk aan de vraag gedurende 
de levertermijn у|Уу2 verminderd met het bestelniveau x2. 


Wij veronderstellen dat degene, die uiteindelijk een beslissing 
moet nemen over de grootte van de bestelhoeveelheid en de hoogte 
van het bestelniveau alleen waarden van xj en za in beschouwing 
neemt die voldoen aan 


хі 2 У1У2 


еп 
0 < хә < У1У2 


Er geldt nu: 


ác. же Op % $ X2 З 3172 
AD AD X] 


Met bovengenoemde veronderstellingen geldt nu 


Usu» 5 1 


Het aantal bestellingen gedurende n perioden is . De be- 


апу, 
X] 


stelkosten gedurende n perioden zijn dus 
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De hoogte van de voorraad bedraagt gemiddeld 40G з l(AD-CD) = 
$ (x1 * x» - y1y2)- 


Zou deze voorraad gedurende n perioden aanwezig zijn, dan zou- 
den de voorraadkosten bedragen: 


(ху + x2 = yiy2)n 


Deze voorraad is echter slechts een bepaald gedeelte van de tijd 
OI A 
п.1. cc aanwezig. 


De voorraadkosten gedurende n perioden bedragen dus 
OI 
jb(x; + X) - yiyo) "яв "0 


Om de verhouding СЕ te berekenen merken wij op 


DK = OC 
en 
X 
AD = tan a > DK = — 
HN У2 
x 
h ÖC Wi 
y2 
OI = OC - ІС» у 22 
yy 
OI y X X Ni 9 oy 
OC X] 72 y2 X] 


De voorraadkosten gedurende n perioden bedragen dus 
Қы A а s 


l Шы e 
jb(xi + хә = у1У2) x 


Het mathematische model voor dit probleem bestaat dus uit een 
stelsel van 2 vergelijkingen, п.1.: 


2 
any (хі + хо - ) 
чі = ME. b dn 
X] X] 
a T Ха А 
LE Se o EE 


X] 
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1.4.4. НЕТ VERZAMELEN VAN GEGEVENS 


De voor het model benodigde gegevens zijn grotendeels van 
kwantitatieve aard. Vaak zal een groot deel van de gegevens in 
verschillende delen van de organisatie in één of andere vorm 
wel beschikbaar zijn. Op de technieken die bij het verzamelen 
van gegevens kunnen worden gebruikt gaan wij hier niet in. 


1.4.5. HET TESTEN VAN HET MODEL 


Men kan het geconstrueerde model op twee wijzen testen. De eer- 
ste manier is een onderzoek naar de logische consistentie van 
het model. Zijn uit het model tegenstrijdige uitspraken af te 
leiden, dan is het logisch inconsistent. Er is dan bij de model- 
constructie een denkfout gemaakt. Dit onderzoek naar de logi- 
sche consistentie vindt geheel plaats binnen het model. 

De tweede manier waarop een model kan worden getest is door te 
onderzoeken of het model met een voldoende mate van nauwkeurig- 
heid de werkelijkheid representeert. Bij het in paragraaf 
1.4.3.2. gegeven voorbeeld kan men bijvoorbeeld de in de prak- 
tijk (in het concrete systeem) gehanteerde waarden voor xj en 
хә invullen. Men verkrijgt dan een voorspelling voor uj en ug. 
De toets bestaat dan hieruit dat men de in de praktijk optre- 
dende waarden van uj en и» vergelijkt met de aldus verkregen 
voorspelde waarden. 

Is het model stochastisch, dan zal men bij het vergelijken van 
de voorspelde en de in het concrete systeem gerealiseerde waar- 
den een statistische toets moeten toepassen. 


1.4.6. HET ZOEKEN VAN DE OPLOSSING 


De volgende fase in de modelcyclus is die van de deductie, d.w.z. 
die van het trekken van conclüsies uit het model. Daarbij kan 
men op twee manieren te werk gaan: door middel van analytische 
methoden dan wel door middel van simulatie. Hier ontmoeten wij 
een vierde indeling van modellen naast de drie in paragraaf 
1.4.3.1. reeds vermelde indelingen, namelijk die van analytische 
modellen versus simulatie-modellen. Bij analytische modellen 
gebruikt men bepaalde wiskundige methoden om de oplossing te 
vinden, bijvoorbeeld de differentiaalrekening. In ons voorbeeld 
kan men zich bijvoorbeeld afvragen onder welke voorwaarden de 
totale kosten (= uj) minimaal worden. 


Bij een gegeven servicegraad u, = s geldt: 


ві З ль YII 
X] 


S = 


S any» + d s? x1? З 
ЕЕЕ 
1 1 
ап 
u) = 22 + lbs? хіп 
1 


шо is minimaal als 


du) _ 
ах pa. 
dus als 
апу? " lbs? ^d. 0 
X] 
gode 6 2ay» 
l 5 b 


en de waarde van и) is dan: 


any 2а 
uj 7 e ب‎ bs? | Se ‘n 
6 


2ay» 
bs? 


> ці = ns v2bay» 


Bij simulatiemodellen wordt de werkelijkheid nagespeeld om een 
oplossing te verkrijgen. Simulatiemodellen kunnen in principe 
zowel deterministisch als stochastisch zijn. 


1.4.7. HET TESTEN VAN DE OPLOSSING 


Bij het testen van de oplossing kan men weer twee wegen bewan- 
delen. Men kan in de eerste plaats de oplossing testen binnen 
het model. In ons voorbeeld zou men waarden voor xj die niet 
ау? 
b 
dergelijke toets gaat men uitsluitend na of de mathematische 
bewerkingen die zijn uitgevoerd om de oplossing te vinden, cor- 
rect zijn. 

In de tweede plaats kan men de oplossing toetsen in de werkelijk- 
heid. In ons voorbeeld zou men de gevonden oplossing bijvoorbeeld 
voor een deel van het assortiment kunnen invoeren. De toets be- 
staat dan hierin, dat men nagaat of dit inderdaad leidt tot la- 
gere kosten dan voor een vergelijkbare groep artikelen, мааг- 
voor men het bestaande systeem van voorraadbeheersing voortzet. 


kd E 3 1 1 A F 
gelijk 211п aan š Киппеп substitueren in uy. В1) een 
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Bij deze toets gaat het om de vraag of het model het concrete 
systeem voldoende nauwkeurig weergeeft, 


1.4.8. IMPLEMENTATIE 


De implementatie zal in het algemeen geleidelijk geschieden. De 
fase van testen zoals in de voorgaande paragraaf beschreven, 
gaat dan geleidelijk over in implementatie. Dikwijls bestaan er 
weerstanden in de organisatie tegen veranderingen. Wij gaan 
daar hier niet op in. 


1.5. MODELLEN ALS BENADERING VAN CONCRETE SYSTEMEN 


Bij het construeren van een model zal men steeds het model ener- 

zijds de werkelijkheid zo goed mogelijk willen laten benaderen 

terwijl het model anderzijds gemakkelijk oplosbaar moet zijn. 

Men zal vrijwel steeds een compromis moeten sluiten, waarbij 

aan beide doeleinden zo goed mogelijk recht wordt gedaan. Men 

kan de werkelijkheid onder andere vereenvoudigen door 

- het weglaten van bepaalde variabelen 

- het aggregeren van variabelen 

- het veranderen van het karakter van variabelen (bijvoorbeeld 
een stochastische variabele vervangen door een deterministi- 
sche, of een discrete variabele door een continue variabele) 

- het veranderen van de relatie tussen twee of meer variabelen 
(bijvoorbeeld het vervangen van een kromlijnig verband door 
een rechtlijnig verband). 


1.6. PROTOTYPE-PROBLEMEN 


Er zijn gedurende het bestaan van de operationele research zeer 
vele problemen opgelost, die men in een aantal prototypen kan 
indelen. Kennis van deze prototypen is buitengewoon nuttig, іп- 
dien men voor een bepaald probleem zelf een model wil constru- 
eren. In zeer vele gevallen kan men zelfs volstaan met het ge- 
bruiken van een bestaand prototype-model. In bepaalde gevallen 
kan men dan van kant en klare computer programma's gebruik па- 
ken bij het zoeken van de oplossing. 

Men kent de volgende prototype-problemen: 

- allocatie-vraagstukken 

- voorraadbeheersingsproblemen 

- wachttijdproblemen 

- netwerkproblemen 

- vervangingsproblemen. 
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In dit boek bespreken wij uitsluitend allocatievraagstukken. 
Voor het oplossen van allocatievraagstukken gebruikt men mathe- 
matische programmeringstechnieken. 


2. LINEAIRE PROGRAMMERING MET 
TWEE VARIABELEN 


2.1. INLEIDING 


Wij geven eerst een voorbeeld van een lineair programmerings- 
probleem. Stel een onderneming heeft bepaalde gegeven capacitei- 
ten van verschillende types machines, een gegeven grootte van 
het personeelsbestand en wellicht eveneens beperkte voorraden 
grondstoffen. De onderneming kan verschillende produkten maken. 
Voor ieder produkt kent men de verkoopprijs, de variabele kosten 
en het beslag per eenheid produkt op de verschillende types 
machines en arbeid evenals de benodigde grondstoffen per eenheid 
produkt. 

Welke produkten (en welke hoeveelheden van die produkten) moet 
men fabriceren indien men streeft naar een zo groot mogelijke 
winst? Dit is een lineair programmeringsprobleem. 

Men kan ook zeggen, dat men uit een groot aantal mogelijke acti- 
viteiten (п.1. alle denkbare produkten in alle denkbare hoeveel- 
heden), een aantal activiteiten moet kiezen op zodanige wijze, 
dat ten eerste aan alle restricties (de beperkte шасһіпе-сарасі- 
teiten, de beperkte hoeveelheden arbeid en grondstoffen) is vol- 
daan en ten tweede een doelstellingsfunctie (de winst) wordt 
gemaximaliseerd. 


2.2. EEN VOORBEELD VAN EEN ONDERNEMING MET TWEE 
PRODUKTEN 


Wij beperken ons aanvankelijk tot een onderneming, die slechts 
twee produkten kan fabriceren. Wij noemen deze produkten A en 
В. Van deze produkten zijn de volgende gegevens bekend: 


produkt A produkt B 
Verkoopprijs 80 59 
Grondstoffen 35 22 


Dekking 45 30 
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Produkt А levert dus een dekking уап f 45,- per eenheid produkt 
en produkt B een dekking van f 30,- per eenheid produkt. Wij 
gaan er in deze opstelling van uit dat alle kosten met uitzon- 
dering van de grondstofkosten als vaste kosten moeten worden 
beschouwd. De onderneming heeft twee machines beide met een ge- 
geven capaciteit, 

Voorts is de hoeveelheid arbeid gegeven; wij gaan er dus van 
uit dat de onderneming gedurende de planperiode geen werknemers 
kan aantrekken of ontslaan. Het capaciteitsbeslag per eenheid 
produkt is gegeven in tabel 2,1, 


Produkt A Produkt B Max. beschikbaar 
Machine І 2 1 120 
Machine 2 1 2 140 
Arbeid 1 1 80 


- tabel 2.1. Capaciteitsbeslag еп beschikbare capaciteiten 
beide uitgedrukt in uren - 


2.3. EEN INTUITIEVE BENADERING 


Bij het zoeken van oplossingen zou men als volgt kunnen redene- 
ren: indien de onderneming alleen produkt A gaat produceren kan 
zij daarvan maximaal 60 eenheden maken (60 is het kleinste van de 
drie getallen 120/2, 140/1, 80/1). De in totaal behaalde dekking 
bedraagt dan 60 х 45 = f 2.700,-. Van produkt B zou men maximaal 
kunnen maken 70 eenheden (70 is het kleinste van de drie getal- 
len 120/1, 140/2, 80/1). De in totaal behaalde dekking bedraagt 
70x 30 г f 2.100,-. Het produceren van 60 eenheden A levert 

dus een hogere dekking dan het produceren van 70 eenheden B. 

Men zou aldus in eerste instantie kunnen besluiten 60 eenheden 
А te gaan produceren. Bij dat produktieprogramma is machine 1 
volbezet; er is echter nog wel ruimte op machine 2, terwijl 
tevens niet alle arbeid is bezet. Van machine 2 zijn nog onbe- 
zet 140 - 60 = 80 uren, terwijl het aantal nog niet benutte аг- 
beidsuren 80 - 60 = 20 bedraagt. Stel пи dat men toch nog één 
eenheid van produkt B zou willen fabriceren. Om na te gaan of 
hierdoor nu de totaal behaalde dekking toeneemt moet men aan 
produkt B de capaciteitskosten van machine 1 toerekenen. Welke 
capaciteitskosten moet men toerekenen? Het antwoord luidt: 

de opportunity-cost, die ontstaat omdat men door B te willen 
maken nu minder van A moet gaan maken 6). Om één eenheid B te 
kunnen maken zal шеп } eenheid A moeten opofferen, dus in plaats 
van 60 eenheden A genoegen moeten nemen met 59,5 eenheden A 
(door } eenheid A minder te maken komt er op machine 1 precies 
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genoeg capaciteit vrij om één eenheid van B te kunnen maken). 
De opportunity-cost als bovenbedoeld is dus f 22,50 per machine- 
uur van machine 1. Dit getal van f 22,50 noemt men de Zmputed- 
cost van machine 1. 

Men kan nu de bijdrage van produkt B tot de totale dekking als 
volgt evalueren: 7 30,- - f 22,50 - 7 7,50. Dit getal van 

f 7,50 per eenheid voor produkt В noemt men de reduced revenue 
van produkt B. 

. Wij hebben dus geconcludeerd, dat door і) eenheid van produkt А 
minder te maken еп in plaats daarvan één eenheid van produkt B 
te produceren de behaalde dekking met f 7,50 toeneemt. Men kan 
dit als volgt eenvoudig controleren: 59,5 х 45 + 1 * 30 = 
270750. 

Men kan nu ook nog eens een halve eenheid А "inruilen!" voor één 
eenheid В: 59 * 45 + 2 ж 30 = 2715,00. 

Ruilt men een halve eenheid A voor één eenheid B, dan betekent 
dit echter wel dat 3/2 uur extra van machine 2 en tevens 5 uur 
extra arbeid moet worden ingeschakeld. Op een gegeven moment 
zal men dus de grens van de maximaal beschikbare hoeveelheid 
arbeid of de maximale capaciteit van machine 2 bereiken. Welke 
van de twee grenzen het eerst wordt bereikt vindt men door te 


vergelijken D en 575 . In de teller van deze breuken staat 


de nog onbezette capaciteit, in de noemer het extra capaciteits- 
beslag veroorzaakt door | eenheid А te ruilen voor één eenheid 


В. Van beide breuken SE. 40 ES? = 53,3 is de eerste de 

4 1/2 3/2 
kleinste. 
Dit houdt in dat men 40 maal een dergelijke ruil kan toepassen 
en dat dan alle arbeid volledig is bezet. Na 40 maal een derge- 
lijke ruil te hebben uitgevoerd hebben wij als produktieprogram- 
ma: 40 eenheden A en 40 eenheden B. De dekking is dan 
40 ж 45 + 40 ж 30 = 3000. 


Dit is tevens het optimale programma, zoals wij in de volgende 
paragraaf zullen laten zien. 


2.4. CONSTRUCTIE VAN HET MATHEMATISCH MODEL 


Vanwege de beperkte produktiecapaciteiten zijn er grenzen aan 
de hoeveelheden, die men van de produkten A en B kan produceren. 
Noemen wij xj = het aantal eenheden van produkt А en х2 = het 
aantal eenheden van produkt B, dan moet gelden: 


2xj + хә 5 120 (machine 1) 
x] + 2x2 < 140 (machine 2) 
X) + X> < 80 (arbeid) 
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Wij zoeken nu zodanige waarden van xj еп xy, dat de waarde van 
de uitdrukking 45xj + 30х2 zo groot mogelijk wordt, terwijl te- 
vens aan bovenstaande ongelijkheden is voldaan. Wij kunnen dus 
voor het in paragraaf 2.2. gestelde probleem het volgende ma- 
thematische model construeren: 


Maximeer 
45хі + 30х2 
Onder de randvoorwaarden 


2X] "N94 5 120 


X1 * 2X2 « 140 
X] + X» < 80 
xy D хә 2 0 


De uitdrukking 45x) + 30x2 noemt шеп de doelstellingsfunctie 
(objective function, goalfunction), de boven vermelde ongelijk- 
heden noemt men de randvoorwaarden of de restricties (restrict- 
ions, constraints). 

Meer algemeen kunnen wij bovenstaand model ook als volgt be- 
schrijven: 


Maximeer 
С]Х] + с2Х2 
onder de randvoorwaarden 
8a11X1 + 8j2X2 S bj 
a21X) + 855X5 < b, 
а31Х1 + аҙ2Х2 < b3 
£i е O, Xo H 
De getallen cj, c5, bj, bo, Ъз, а)), a12, азі, 855, аз), 232 noemt 


men de parameters van het model. De numerieke waarden van deze 
parameters liggen vast, 


2.5. DE GRAFISCHE OPLOSSINGSMETHODE 


Wij kunnen op zeer eenvoudige wijze grafisch aangeven welke pro- 
duktieprogramma's technisch wel en niet mogelijk zijn. In het 
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x], X2-vlak tekenen wij daartoe eerst de rechte lijn met de 
vergelijking 2xj + xj = 120. Zie figuur 2.1. 


x» 


120 


` (40, 60) 


machine 1 


- figuur 2.1. Gebied van mogelijke oplossingen, als alleen 
machine | in beschouwing wordt genomen - 


Het punt P (20, 40) in figuur 2.1. representeert een produktie- 
programma xj = 20 en x; = 40. Daarvoor geldt 2 х 20 + І «40 < 120, 
hetgeen inhoudt dat dit produktieprogramma wat betreft de capa- 
citeit van machine 1 mogelijk is. Het punt Q (40, 60) represen- 
teert een produktieprogramma dat niet mogelijk is, omdat dan de 
capaciteit van machine 1 wordt overschreden: 2 «40 + 1*60 > 120. 
Verder zijn natuurlijk negatieve waarden van xj en x2 niet mo- 
gelijk. Wij concluderen dus, dat alleen die produktieprogram- 
ma's, die kunnen worden weergegeven door een punt binnen of 
juist op de grens van driehoek OAE, gezien de capaciteit van 
machine | mogelijk zijn. 
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In figuur 2.2. hebben we de drie lijnen 


EET Ж хо е 120 
жі» 2х2 = 140 
X) " де * 80 


alle drie getekend. Alleen produktieprogramma's die kunnen wor- 
den weergegeven door een punt binnen of juist op de grens van 
de veelhoek OABCD zijn voor de onderneming technisch gezien mo- 
gelijk. 

Men noemt het gebied OABCD het "area of feasible solutions". 
Een "feasible solution" is een oplossing, die aan alle restric- 
ties voldoet (in ons voorbeeld is een feasible solution een ор- 
lossing die technisch mogelijk is). 


Het probleem is nu dat wij uit de verzameling van feasible so- 
lutions, die oplossing moeten kiezen, die de grootst mogelijke 
dekking geeft. Wij beschouwen daartoe figuur 2.3. 


machine | arbeid machine 2 


- figuur 2.2. Het "area of feasible solutions" - 
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In figuur 2.3. beschouwen wij eerst de lijn Lj. De lijn Lj heeft 
als vergelijking: 


45хі + 30х2 - 1350 


Anders gezegd: de verzameling уап produktieprogramma's met een 
dekking van 1350 wordt weergegeven door de lijn Lj (alleen dat 
gedeelte van lijn Lj, dat in het eerste kwadrant ligt). De ver- 
zameling van produktieprogramma's met een dekking van 1800 
wordt weergegeven door de lijn L5. Lijn Ly is evenwijdig met 
lijn Lj. Trekt men een willekeurige lijn evenwijdig met Lj, dan 
representeert deze lijn een verzameling oplossingen die alle de- 
zelfde dekking geven. Hoe groter de afstand van de oorsprong 
tot de betreffende lijn, des te groter is deze dekking. Lijn Lu 
bevat nog juist één punt van het area of feasible solutions, 
n.l. punt B. Punt B is derhalve van alle feasible solutions 
degene die de hoogste dekking heeft. 


De exacte coordinaten van B vinden wij door het snijpunt te be- 
palen van de lijnen overeenkomend met de restricties voor ma- 
chine 1 en arbeid (zie figuur 2.2.). Punt B is dus het snijpunt 
van 


120 


2X] + X2 
en 

Ny $ x, * BQ 
Oplossen van dit stelsel levert xj = 40, х2 т 40. Het optimale 


produktieprogramma is dus 40 eenheden B. 
De behaalde dekking is dan 45 х 40 + 30 х 40 = 3000. 


get d? `. 2.3; 


2.6. AANNAMES VAN LINEAIRE PROGRAMMERING 


De methode van de lineaire programmering is gebaseerd op een 
aantal veronderstellingen. Deze zijn: 


a. Lineariteit van de doelstellingsfunctie 

De doelstellingsfunctie moet de vorm суху + CoX5 + .... + C 
hebben. 

Doelstellingsfuncties als bijvoorbeeld: 


Cixi + cox]? + CKD TT oor 
of 


C1X] + С2Х1Х2 * ооо о о 


voldoen niet aan deze eis van lineariteit. 

De eis van lineariteit betekent, dat: 

le de bijdrage van activiteit і aan de doelstellingsfunctie 
onafhankelijk is van de niveaus van de andere activiteiten 
en 

2e deze bijdrage recht evenredig is met het niveau van activi- 
teit 1. 
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bijdrage уап xj aan de doelstellingsfunctie 


X. 
l 


- figuur 2.4. Een lineair en een niet-lineair verband tussen 
x. en de bijdrage van x. aan de doelstellingsfunctie - 


Van de in figuur 2.4. gegeven verbanden tussen хі en de bijdrage 
van хі aan de doelstellingsfunctie, leidt dus alleen het door 
lijn L aangegeven verband tot een lineair programmeringsmodel. 


b. Lineariteit van de restricties 

Het beslag van de verschillende activiteiten op de beperkte са- 

paciteiten moet eveneens leiden tot een lineair verband. Dit 

houdt in: 

le het capaciteitsbeslag van activiteit і per eenheid produkt 
is niet afhankelijk van de vraag in hoeverre andere activi- 
teiten eveneens beslag leggen op deze capaciteit en 

2e dit capaciteitsbeslag is niet afhankelijk van het niveau 
van activiteit 1. 


с. Deelbaarheid 

Wij hebben tot dusverre steeds verondersteld, dat de beslissings- 
variabelen xj en хә continue variabelen zijn. Zijn één of meer 
beslissingsvariabelen discreet (bijvoorbeeld: xj moet een niet- 
negatief geheel getal zijn), dan is aan de veronderstelling van 
deelbaarheid van iedere activiteit niet voldaan. Veronderstellen 
wij dat de in dit hoofdstuk besproken onderneming een raffinade- 
rij is, dan is aan de eis van deelbaarheid voldaan. 


Tenslotte merken wij nog op dat het lineaire programmeringsmo- 
del kan worden gekarakteriseerd als: een deterministisch opti- 
malisatiemodel. 
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2-5. Gegeven is het volgende LP-model. 
Maximeer 


10х1 + 12х2 


onder де voorwaarden 


Zei + Ах; 5 60 
хі + хә = 25 
203 KZ, РФ 40 
РО хо о У 


Geef іп een grafiek het gebied van de feasible solutions aan. 
Hoe luidt de optimale oplossing? 

Welke waarde heeft de doelstellingsfunctie in de ant make 
oplossing? 

Een der restricties is in dit model overbodig (redundant). 
Welke restrictie is dat? 


2-6. Gegeven is het volgende LP-model. 
Maximeer 


ERT m R 


onder de voorwaarden 


Xr» $2 


Xj + Xo 10 


х1 e 0, Х2 > 0 


Heeft dit model feasible solutions? 
Heeft dit model een optimale oplossing? 


2-7. Een bouwonderneming moet voor een gebied met een opper- 
vlakte van ЄЄп hectare een bouwplan opstellen. Men overweegt de 
bouw van woningen in acht verdiepingen en de bouw van winkels 
In twee verdiepingen. Daarbij heeft men te maken met de vol- 
gende beperkingen: 


Ж 


Op grond van de bestaande voorschriften mag de verhouding 
vloeroppervlak : grondoppervlak voor woningen niet groter 
zijn dan 2. Ook voor de winkels geldt dat de verhouding tus- 
sen het vloeroppervlak en het grondoppervlak niet groter mag 
zijn dan 2. Daarnaast moet binnen de beschikbare oppervlakte 
van één hectare een —— worden aangelegd en wel 
20, dat voor iedere 50 m? winkelvloeroppervlak één parkeer- 
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plaats met een grondoppervlak van 25 п? wordt gerealiseerd. 

2. Tenminste 30% van de beschikbare oppervlakte moet worden 
besteed aan groenvoorzieningen. Dit betekent dat het grond- 
oppervlak van de te bouwen woningen (dat gelijk is aan 1/8 
van het te bouwen vloeroppervlak, omdat de woningen worden 
gebouwd in acht verdiepingen) plus het grondoppervlak van de 
te bouwen winkels (dat gelijk is aan de helft van het te 
realiseren vloeroppervlak, omdat winkels worden gebouwd in 
twee verdiepingen) plus de parkeerplaatsen tezamen ten hoog- 
ste 70% van het beschikbare grondoppervlak mogen beslaan. 

3. De bouw van woningen en winkels veroorzaakt een additionele 
luchtvervuiling. Uit metingen is gebleken dat één m? win- 
kelruimte 3 eenheden luchtvervuiling veroorzaakt, terwijl 
dit voor één m? woonruimte slechts 1 eenheid bedraagt. De 
maximaal toelaatbare extra luchtvervuiling is gesteld op 
30.000 eenheden. 

4. De bouw van de woningen en de winkels moet door de bouwon- 
derneming worden voorgefinancierd. In het onderhavige pro- 
jekt wenst men ten hoogste f 40 miljoen te investeren. Voor 
de bouw van één m? woonruimte is een voorfinanciering nodig 
уап f 3000,--; voor winkelruimte is dit 7 4000,--. 
(Aanwijzing: formuleer deze restrictie in eenheden van 
f 1000,77.) 

5. De gemeente wenst aan het gebied een duidelijke winkelfunc- 
tie toe te kennen; het aantal m? winkelruimte moet tenmin- 
ste 1000 plus de helft van het aantal te realiseren m? 
woonruimte bedragen. 


De bouwonderneming verwacht met de bouw van woningen en winkels 
een brutowinst van f 200,-- per m? woningruimte en f 400,-- per 
m2 winkelruimte te kunnen behalen. De bouwonderneming wil bin- 
nen de gegeven restricties een bouwplan opstellen, dat een zo 
groot mogelijke brutowinst oplevert. 


a) Formuleer dit probleem als een LP-probleem. 
b) Los het probleem grafisch op. 


2-8. Gegeven is het volgende LP-probleem. 
Maximeer 


2X4 + X» 


onder de voorwaarden 


xj + хо S 8 


ІЛ 


14 
4 


X1 Ф 2х2 


ІМ 


X1 + 2x2 


X15 х2 > 


о 


Geef іп een grafiek het gebied van de feasible solutions еп de 
optimale oplossing aan. ` 
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2-9. Een oliemaatschappij beschikt over drie soorten dieselolie 
(aan te duiden met А, В en С). Het soortelijk gewicht en het zwa- 
velgehalte van deze drie soorten dieselolie is gegeven in onder- 
staande tabel. 


B 
soortelijk gewicht 0,90 0,92 0,86 
zwavelgehalte іп 2 0,28 0,38 0533 


Deze oliemaatschappij wenst nu door menging van deze drie soorten 
dieselolie een dieselolie te verkrijgen met een soortelijk gewicht, 
dat ten minste gelijk is aan 0,89 en met een zwavelgehalte van 

ten hoogste 0,347. | 

De oliemaatschappij kan onbeperkte hoeveelheden van de soorten 

A, B en C inkopen, tegen prijzen van respectievelijk 13 cent, 

9 cent en 10 cent per liter. 

Men wenst het mengsel zodanig samen te stellen dat het zo goed- 
koop mogelijk is, maar toch aan de eisen met betrekking tot soor- 
telijk gewicht en zwavelgehalte voldoet. 


a) Formuleer dit probleem als een LP-probleem met twee variabelen. 
b) Los het probleem op met behulp van de grafische oplossings- 
methode. 


3. SCHADUWPRIJZEN, PARAMETRISCH 
PROGRAMMEREN EN 
GEVOELIGHEIDSANALYSE 


3.1. SCHADUWPRIJZEN 


Wij veronderstellen thans, dat de onderneming uit hoofdstuk 2, 
overweegt de capaciteit van machine 1 uit te breiden zo dat 1 
uur per week extra ter beschikking komt. In plaats van 120 uur 
staat de onderneming dan 121 uur ter beschikking. Alvorens tot 
deze capaciteitsuitbreiding over te gaan wil de onderneming we- 
ten, in hoeverre de behaalde dekking toeneemt door het ter be- 
schikking krijgen van | uur extra van machine 1. In figuur 3.1. 
hebben wij de nieuwe situatie in beeld gebracht. Het optimum 
verschuift nu van B naar Bj. Bj is het snijpunt van 


2X] ту» " 121 
en 


x] * x2 - 80 


De coordinaten van Bj zijn dus xj = 41 en xj = 39. 
De met Bj overeenkomende dekking bedraagt: 


45 * 41 + 30 * 39 = 3015. 


Ten opzichte van de oude situatie is de behaalde dekking dus 
toegenomen met 15. 

Men zegt nu: de schaduwprijs van restrictie l is gelijk aan 15. 
Anders gezegd: de schaduwprijs van een restrictie is gelijk aan 
de marginale toeneming van de doelstellingsfunctie als gevolg 
van een marginale verruiming van die restrictie, onder de ver- 
onderstelling, dat alle andere restricties ongewijzigd blijven. 
De laatste toevoeging (de ceteris paribus clausule) is essen- 
tieel voor de definitie van het begrip schaduwprijs. 

De onderneming zal nu overgaan tot uitbreiding van de capaciteit 
van machine 1, indien de kosten van deze capaciteitsuitbreiding 
kleiner zijn dan de schaduwprijs van restrictie l. 
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machine 1 arbeid machine 2 


- figuur 3.1. De schaduwprijs van restrictie І - 


Men kan de schaduwprijs ook interpreteren als de opportunity 
costs vanhet niet ter beschikking hebben van één extra capaci- 
teitseenheid van machine 1. Zijn deze opportunity costs hoger 
dan de werkelijke kosten, die men voor het verkrijgen van één 
extra capaciteitseenheid zou moeten maken, dan moet men tot uit- 
breiding overgaan. 


3.2. РАКАМЕТКІ5СН PROGRAMMEREN 


In paragraaf 3.1. varieerden wij de constante in het rechterlid 
van restrictie 1 met één eenheid. Wij varieren thans deze con- 
stante, bj, over een groot gebied, п.1. van 0 tot 200. Wij ver- 
onderstellen daarbij dat er overigens in de formulering van het 
model geen wijzigingen optreden. Wij willen thans onderzoeken 
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in hoeverre de waarde уап de doelstellingsfunctie verandert als 
wij bj stapsgewijs varieren уап 0 tot 200. 

Deze veranderingen geven wij aan in tabel 3.1. Wij beschouwen 
nu figuur 3.2. In figuur 3.2. hebben wij getekend de ongewij- 
zigde restricties 2 en 3. Indien wij restrictie 1 buiten be- 
schouwing laten is OFCD het area of feasible solutions. 

Als bj = 0 moet gelden xj = O en хә = 0. Wij noemen de waarde 
van de doelstellingsfunctie 2 = 45xj + 30x2. Als bj = O is dus 
ook z = 0. 


machine 2 


- figuur 3.2. Parametrisch programmeren - 
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bi X] X2 2 | А2 
0 0 0 0 Sege 
1 1 30 30 
2 0 2 60 30 
69 0 69 2070 
70 0 70 2100 30 
71 2/3 69 2/3 2120 20 
72 4/3 69 1/3 2140 20 
99 94/3 60 1/3 2680 
100 20 60 2700 20 
101 21 59 2715 15 
102 22 58 2730 15 
159 79 1 3585 
160 80 0 3600 15 
161 80 0 3600 
200 80 0 3600 


- tabel 3.1. 2 еп Az als functie van bj - 


In figuur 3.2. hebben wij ook getekend de lijn 2x; + x, = 20 еп 
de lijn 45x; + 30х; = 1350. Gemakkelijk is nu in te zien dat in- 
dien bj = 20, P het optimum is. Evenzo zal, als bj = 1, het op- 
timum zijn het snijpunt van de lijn 2xj + xj = 1 met de y-as. 
Dit snijpunt is xj = 0, х2 = 1 еп de waarde van z = 45x] + 30x, 
bedraagt in dat punt 30, Indien bj = 70 wordt punt D het opti- 
male punt. Geeft men bj een nog grotere waarde, bijvoorbeeld 71, 
dan verschuift het optimale punt niet langer langs de х2-ав, maar 


3600 


2700 


2100 


m figuur 2; 3 e 


helling 0 | 


«— helling 15 


€—— helling 20 


€—— helling 30 


70 100 160 


2 als functie van bj - 


200 


39 


- bj 
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langs de restrictie van machine 2. Indien bj = 71 is het snij- 
punt van 2xj + хо = 71 en xj + 2x, = 140 het optimum. De coor- 
dinaten van dit snijpunt zijn 

x| = 2/3, gas 69 2/3 еп z 9 45*2/3 + 30 х 69 2/3 = 2120. 
Indien bj = 100 is C het optimale punt. Bij waarden van bj > 100 
verschuift het optimale punt langs de restrictie van de arbeid. 
Indien bj = 160 is F het optimale punt. Geeft men bj nog grotere 
waarden, dan blijft toch F het optimale punt. De relatie tussen 
z en bj, die wij in tabel 3.1. hebben gevonden, hebben wij in 
figuur 3.3. grafisch weergegeven. In figuur 3.3. is de helling 
van de grafiek niets anders dan de marginale toeneming van z als 
gevolg van een marginale toeneming van bj. Dit is de schaduw- 
prijs van de restrictie. Wij zien dat voor bj = 120 deze scha- 
duwprijs gelijk is aan 15, hetgeen wij ook in paragraaf 3.1. 
reeds hebben geconstateerd. 


3.3. GEVOELIGHEIDSANALYSE 


Vele van de gegevens, die voor een lineair programmeringsmodel 
nodig zijn, zijn niet nauwkeurig van te voren bekend maar be- 
rusten op schattingen. Zo zal bijvoorbeeld in ons voorbeeld de 
verkoopprijs van produkt A een schatting van de verkoopprijs 
in de planperiode kunnen zijn. 

Daarom is het zinvol de vraag te stellen in hoeverre het opti- 
male produktiéprogramma verandert, als de inputgegevens veran- 
deren. Wij beperken ons tot veranderingen in de coefficienten 
van de doelstellingsfunctie, hoewel men natuurlijk dezelfde 
vraag voor veranderingen in de elementen van matrix A kan stel- 
len. 


Wij vragen ons dus af of het optimale programma verandert indien 
cj in plaats van 45 een andere waarde, zeg 44 krijgt. Beschouwen 
wij figuur 2.3. dan zien wij, dat daardoor de lijnen die de 
richting van de doelstellingsfunctie aangeven draaien en wel 
(bij een afneming van cj) tegen de wijzers van de klok in. Bij 
een kleine wijziging van cj blijft punt B het optimale punt. 
Dat wordt anders indien de waarde van cj zodanig wordt verkleind, 
dan lijnen, die de richting van de doelstellingsfunctie aangeven 
evenwijdig gaan lopen met lijnstuk ВС. In dat geval zijn alle 
punten op het lijnstuk BC optimaal. Dit is dus het geval als de 
lijnen 

суху + 30x, = constant 


evenwijdig zijn met de lijn xj + xj = 80. Beide lijnen zijn even- 
wijdig dan en slechts dan indien 


e (1: 90 t 1 


D 


ofwel dan еп slechts dan als cj = 30. Indien с) groter is dan 
45 kunnen wij een soortgelijke redenering toepassen: het gehele 
lijnstuk BA is optimaal, dan en slechts dan als de lijnen 


сухі + 30х2 = constant 


evenwijdig zijn met de lijn 2х) + x = 120, 

Dit is zo d.e.s.d. als су : 2 = 30 : 1 ofwel cj = 60. 

Wij concluderen dus dat het optimale programma ongevoelig is 
voor wijzigingen van с), zolang deze liggen binnen het interval 
30 < cj < 60. 

Bovenstaande analyse noemt men een gevoeligheidsanalyse. 


3.4. VRAAGSTUKKEN 


3-1. De verkoopafdeling van de in vraagstuk 2-1 bedoelde onder- 
neming deelt mee maximaal 400 stuks van produkt a te kunnen 
verkopen. 

a. Voeg deze verkooprestrictie als derde restrictie toe aan het 
іп vraagstuk 3-1 geformuleerde model. Wat is nu de optimale 
oplossing? Hoe groot is nu de waarde van de doelstellings- 
functie in de optimale oplossing? 

b. Hoe groot is de schaduwprijs van de verkooprestrictie? 

c. Indien b4 (dat is de parameter in het rechterlid van de ver- 
kooprestrictie) de waarden 0, 1, 2, ....., 1000 aanneemt, 
welke waarden krijgt de doelstellingsfunctie dan? Maak een 
grafiek, waarin z (dat is de waarde van de doelstellings- 
functie) wordt geschetst als functie van ba, 


3-2. Gegeven is het volgende LP-model: 
Maximeer 


2X] * X5 


onder de voorwaarden 


кү * 548 І 
Xi t 2X5 & 7 II 
мку 2X» < 4 III 
x 2 U, 0: 


a. Los het probleem op met behulp van de grafische oplossings- 
methode. 
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b. Hoe groot is de schaduwprijs уап restrictie II? 

с. Welke waarden krijgt de doelstellingsfunctie als by (dat is 
de parameter in het rechterlid van restrictie II) de waarden 
Oe nee a » 20 aanneemt? Maak een grafiek, waarin 2 (dat 
is de waarde van de doelstellingsfunctie) wordt getekend als 
functie van b». 


3-3. Wij bekijken opnieuw het LP-model van vraagstuk 2-1 (wij 
laten nu dus de in vraagstuk 3-1 geïntroduceerde verkooprestric- 
tie buiten beschouwing). 

Verandert het optimale produktiepakket als de dekkingsbijdrage 
van produkt A verandert van 7 200,- per stuk іп f 160,- per stuk? 
Verandert de waarde van de doelstellingsfunctie door deze wijzi- 
ging? 

Dezelfde vragen voor een verandering van de dekkingsbijdrage van 
produkt A in f 120,- per stuk. 

Binnen welke grenzen Кап с) (dat is de dekkingsbijdrage per stuk 
van produkt A) veranderen zonder dat de samenstelling van het 
optimale produktiepakket verandert? 


3-4. Gegeven is het volgende LP-probleem: 


Maximeer 
2X4 + х2 

onder de voorwaarden 
x4 + 2x2 < 100 1 
X4 * х28 70 Lë 
X1 5 29 III 
хі, X2 2 0 


а. Los dit probleem op met behulp van de grafische oplossings- 
methode. 

b. Hoe groot is de schaduwprijs van restrictie III? 

c. Welke waarden krijgt x (dat is de waarde van de doelstellings- 
functie) als ba (dat is de parameter іп het rechterlid van 
restrictie ІІІ) de waarden 0,1,2,...,200 aanneemt? 

Maak een grafiek, waarin 2 wordt getekend als functie van ba, 


3-5. Beschouw nogmaals vraagstuk 2-7. De bouwonderneming wenst 

in het onderhavige project niet meer dan f 40 mln te investeren, 
omdat de financieringsmiddelen nu eenmaal beperkt zijn. Stel dat 
de bouwonderneming onbeperkt geld kan lenen tegen een interest- 
vergoeding уап 5% gedurende de looptijd van het project (als de 
onderneming f 1000,-- leent, dan moet men daarover dus een inter- 
est ten bedrage van f 50,-- betalen). 


a. Is het voor de bouwonderneming aantrekkelijk tegen deze іпіег- 
estvoet geld te lenen? Met welk bedrag neemt de winst toe 
(с.4. af) als men 7 1000,- leent? 
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b. Als het aantrekkelijk is om geld te lenen, welk bedrag moet 
de bouwonderneming dan lenen als men de winst wil maximeren? 


3-6. Wij bekijken opnieuw het vraagstuk van klokkenfabriek Frisia 
B.V. Verandert het optimale produktiepakket als de brutomarge 
voor de Friese staartklok verandert van f 160,- per klok in 

f 180,- per klok? Verandert de waarde van de doelstellingsfunctie 
door deze verandering? Zo ja, kunt u dan aangeven met welk bedrag 
de waarde van de doelstellingsfunctie verandert? 

Binnen welke grenzen Кап c, (dat is de coëfficiënt van Friese 
staartklokken in de doelstellingsfunctie) variëren, zonder dat 

de samenstelling van het optimale produktiepakket verandert? 


4. DE SIMPLEX-METHODE 


4.1. INLEIDING 


Het is duidelijk, dat de grafische oplossingsmethode alleen kan 
worden toegepast als het aantal variabelen twee bedraagt. Is 

het aantal variabelen groter, dan moet men een andere methode 
toepassen, 

Aangezien het optimum altijd een hoekpunt van het area of fea- 
sible solutions is, zou men voor een probleem met п variabelen еп m 
restricties alle snijpunten van de (m*n) hypervlakken (m hyper- 
vlakken voor de m restricties, n hypervlakken voor de additio- 
nele restricties хі > 0), kunnen uitrekenen. In totaal verkrijgt 


mtn ° e 
men dan ( М ) snijpunten. Een gedeelte daarvan correspondeert 


niet met een feasible solution. Men zou dan vervolgens door 
substitutie van de coordinaten van deze hoekpunten in de restric- 
ties moeten nagaan welke hoekpunten overeenkomen met feasible 
solutions. Tenslotte zou men de coordinaten van de hoekpunten, 
die overeenkomen met feasible solutions moeten substitueren in 

de doelstellingsfunctie. Het hoekpunt waarvan de coordinaten де 
hoogste waarde van de doelstellingsfunctie geeft is dan het op- 
timale punt. 


Bovenstaande methode is zeer omslachtig en geeft als n en m 
groot zijn aanleiding tot zo veel rekenwerk, dat de meest mo- 
derne computers duizenden jaren nodig zouden hebben voor het 
oplossen van één probleem. Daarom is behoefte aan een andere 
methode, Deze methode - de simplexmethode - is reeds іп 1947 
door G.B. Dantzig ontwikkeld. 


4.2. EEN NUMERIEK VOORBEELD 


Wij laten eerst aan de hand van een numeriek voorbeeld stap voor 
stap zien hoe de simplexmethode werkt. Wij kiezen daarvoor het 
probleem: 
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maximeer 
45хі + 30х; 
onder de voorwaarden 
cK) + жз 120 
кр * 2x4 š FAD 
ж v ху 0 
хр «9, %2 2 Q 
De eerste stap is, dat wij dit probleem herschrijven zó, dat 
alle restricties de gedaante van gelijkheden in plaats van оп- 


gelijkheden krijgen. Wij doen dit door aan iedere < - restrictie 
Eén extra variabele toe te voegen: 


ташы XO X3 e 120 
Xj * 2X5 + Xu = 140 
KES MI + xo; = 80 


Aan de variabelen хз, xi en xg stellen wij eveneens de eis, dat 
zij niet negatief kunnen zijn. Daarmee kunnen wij de ongelijk- 
heden vervangen door gelijkheden. De variabelen X3, xy en Ze 
noemt men spelingsvariabelen (slack variables). Men kan x3 in- 
terpreteren als het aantal ongebruikte capaciteitseenheden van 
machine 1. Aan x, en xg kan men een analoge interpretatie geven. 
Ook de doelstellingsfunctie schrijven wij als een gelijkheid, 
namelijk als 


E - 839,  30х2 = 0, 


Wij hebben nu ons probleem herschreven als: 


£ 7 859, * 20% = 0 
oN T $ J. Ni DK e = 120 
ЙІ т 4 * Xi - ]40 
ax де + x; = 80 


Dit is een stelsel van 4 vergelijkingen met 6 onbekenden. 

Dit stelsel is onbepaald, dat wil zeggen het heeft een oneindig 
groot aantal oplossingen. Eén oplossing daarvan is zeer gemak- 
kelijk te vinden. 

Door te stellen xj = O en xy = O is het stelsel wel bepaald. 
Wij vinden dan z = 0, хз = 120, x, т 140, xg = 80. 
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Wij noemen een dergelijke oplossing van een stelsel уап (m+!) 
vergelijkingen met (т+п+1) variabelen, die men vindt door wille- 
keurig n variabelen gelijk aan nul te stellen een bastsoplossing 
(basic solution). 


Wij herschrijven nu het stelsel van vergelijkingen in het on- 
derstaande tableau 


- tabel 4.1. Het eerste simplex-tableau - 


Dit tableau moet men lezen net als het hierboven gegeven stelse- 
van vergelijkingen. De eerste regel moet men dus lezen als: 


a em 45x Т 30x» = 0 
de tweede regel als: 
AK + хә + Xs = 120 enz. 


Іп de eerste kolom van het eerste simplex-tableau hebben wij 
bovendien de zgn. bastsvartabelen aangegeven. Wij zeggen, dat 
Z, X3, Х, en xs de bij dit tableau behorende basisvariabelen 
zijn. Stelt men xj = 0 en xs = 0 dan verkrijgt men als basis- 
oplossing 2 з 0, хз = 120, х, = 140 en хс з 80. Deze waarden 
kan men aflezen in de meest rechtse kolom van het simplexta- 
bleau. Deze kolom duiden wij aan met de letters RL (rechterlid) 
of RHS (right-hand-side). 

Bekijkt men de kolommen van het simplex-tableau afzonderlijk, 
dan bevatten de kolommen уап 2, хз, x, en xg alle slechts één 
cijfer ongelijk nul. Dit cijfer is dan een 1, terwijl de aldus 
verkregen vier cijfers 1 op verschillende regels staan. Omdat 
de kolommen van de basisvariabelen alleen nullen en enen bevat- 
ten kan men de waarden van de basisvariabelen direct aflezen 
in het rechterlid. 


De tweede stap is пи, dat men een andere basisoplossing zoekt, 
die voldoet aan twee eisen: 


4] 


le de nieuwe basisoplossing moet eveneens feasible zijn 

2e de nieuwe basisoplossing moet leiden tot een waarde van z 
die niet lager is dan de waarde van z in de bestaande basis- 
oplossing. 


Wij bekijken daartoe eerst de eerste regel van het simplexta- 
bleau 2 - 45x; - 30x, = 0. 

In de bij dit tableau behorende oplossing is xj = 0 en x» = 0. 
Een oplossing waarbij bijvoorbeeld xj > 0 en x, = 0 is, leidt 
dus tot een grotere waarde van 2 dan de huidige waarde z = 0. 
Wij proberen daarom een basisoplossing te vinden met een waarde 
уап xj > 0. In plaats van xj = 0 en xy = 0 moet men nu naast 

хә = 0 nog een andere variabele gelijk aan nul stellen. Indien 
wij ons niet willen beperken tot basic feaszbie solutions zou- 
den wij willekeurig één bepaalde variabele gelijk aan nul kunnen 
stellen. Wij wensen ons echter uitdrukkelijk te beperken tot 
basic feasible solutions. Wij berekenen daarom de waarden van 
de drie breuken: 


120 140 80 


p" 3 
š і суб АБО ; Ka 
De kleinste van deze drie 15 TS Dit betekent, dat als wij 


de waarde van xj zoveel mogelijk willen vergroten de eerste 

restrictie de bottleneck gaat vormen, niet de restricties 2 of 

3. Dit houdt tevens in dat dan x3, de spelingsvariabele van 

restrictie 1, nul wordt. Wij voeren nu op het simplextableau 

van tabel 4.1. een aantal bewerkingen uit. Na uitvoering van 

deze bewerkingen is xj een basisvariabele geworden, terwijl хз 

niet langer een basisvariabele zal zijn. 

Deze bewerkingen zijn: 

le wij delen de tweede rij door 2, dit is niets anders dan de 
vergelijking 


2X] "INS ". 3 Be 120 
schrijven als 
хі + іх; + ix3 = 60 


Het simplextableau ziet er nu uit als aangegeven іп tabel 
4,4 

2e Wij vermenigvuldigen nu de tweede rij met 45 en tellen deze 
rij bij de eerste op. Dit is niets anders dan de tweede 
vergelijking van het stelsel met 45 te vermenigvuldigen.en 
optellen bij de eerste vergelijking. ; 
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- tabel 4.2. Het tweede simplex-tableau in statu 
nascendi - 


Het tableau wordt dan zoals aangegeven in tabel 4.3.: 


- tabel 4.3. Het tweede simplex-tableau nog steeds іп 
statu nascendi - 


Зе Wij vermenigvuldigen пи de tweede rij met -l en tellen die 
op bij de derde rij. Vervolgens vermenigvuldigen wij de 
tweede rij met -1 en tellen die op bij de vierde rij. Wij 
krijgen dan het tableau van tabel 4.4. 


- tabel 4.4. Het tweede simplex-tableau - 


Wij hebben nu een nieuwe basic feasible solution gevonden, n.l.: 
z = 2700, tjes 60, xš 057 x3 "0, "x е 60: BE 20 


Door het uitvoeren van bovengenoemde stappen hebben wij een 
stelsel van vier vergelijkingen met zes onbekenden verkregen, 
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dat equivalent is met het stelsel weergegeven іп het eerste 
tableau. In het tweede simplextableau is xj een basisvariabele. 
De bij xj behorende kolom luidt: 


ооо 


Wij hebben van het oorspronkelijke element іп de tweede rij van 
de kolom behorende bij xj (het getal 2) een 1 gemaakt door te 
delen door 2. Vervolgens hebben wij de rest van de kolom beho- 
rende bij xj "schoongeveegd", d.w.z. wij hebben zodanige bewer- 
kingen toegepast, dat in de rest van de kolom behorende bij xj 
nu nog uitsluitend nullen staan. Het element in de tweede rij 
іп de kolom behorende bij xj in het eerste tableau noemt men 
het pZvot-element. Wij hebben in tabel 4.1. het pivot-element 
omcirkeld. 


De derde stap is nu, dat wij ons afvragen of wij de optimale 
oplossing reeds hebben gevonden. De vraag luidt dus: is de basic 
feasible solution van het tweede tableau (tabel 4.4.) de opti- 
male oplossing of is er een andere basic feasible solution met 
een waarde van z 2 2700? Wij beschouwen hiertoe de eerste rij: 


z - 7іхә + 221хз = 2700 
ofwel 


z з 2700 + 7ixo we 223x3 


In de basisoplossing van het tweede tableau is х; = 0 еп x3 = 0. 
Indien wij nu x een waarde > O geven neemt de waarde van 2 

toe en wel met 7} voor iedere eenheid х;. In het eerste tableau 
luidde de eerste rij z = 45x + 30х;. De bijdrage van x, aan de 
doelstellingsfunctie was toen 30. In het tweede tableau is de 
bijdrage van xy aan de doelstellingsfunctie nog slechts 71. Dit 
komt doordat men om een éénheid x, te kunnen produceren Í een- 
heid xj moet opgeven. Daardoor neemt de waarde van de doelstel- 
lingsfunctie niet toe met 30 maar slechts met 30 - 1 ж 45 = 7l. 
(Men herleze nu ook de intuïtieve benadering van paragraaf 
2535). 

Dat men om 1 eenheid x te kunnen produceren Í eenheid xj moet 
opgeven leest men uit de tweede rij: 


хі + іх) + іхз = 60 


In de basisoplossing is xj = 60, х; = 0 en хҙ = 0. Wil men 
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хә = | maken, dan moet men de waarde van xj wel verkleinen en 
wel tot 59}, aangezien x3 niet negatief mag worden. Het element 
van kolom x, in de rij corresponderend met basisvariabele х] 
geeft dus aan hoeveel eenheden xj men moet opgeven om 1 eenheid 
van X te kunnen produceren. 

Evenzo geven de elementen van х; in de rijen corresponderende 
met de basisvariabelen x, en xg aan hoeveel eenheden van x, 
respectievelijk xs men moet opgeven om één eenheid van x, te 
kunnen produceren: 


3/2х» т. 1/2х3 + m." BO 


In де basisoplossing behorende bij het tweede tableau is 

x, = 80. Wil men x2 = 1 maken, dan zal ofwel x, < 80 en хз = 0 
ofwel xi = 80 en хз > 0 moeten zijn. De laatste mogelijkheid 
(x, = 80 en x3 > O is niet aantrekkelijk: men zou dan namelijk 
via de vergelijking van de tweede rij ook xj verder moeten la- 
ten dalen dan hierboven reeds besproken; dit zou inhouden dat 
wij er niet meer op moeten rekenen dat de waarde van de doel- 
stellingsfunctie toeneemt met 30 - | ж 45 = 7} als x) toeneemt 
met één eenheid. Wij houden dus x3 = 0 en zien dat wij dan om 
хә met één eenheid te kunnen laten toenemen xi moeten laten 
dalen met 3/2. Меп kan dit ook inzien als volgt: ruilt шеп één 
eenheid x voor 1 eenheid xj, dan is daarvoor van machine 2 
extra benodigd 2 - | = 3/2 uur. De waarde van de spelingsvaria- 
bele xi zal dus afnemen met 3/2 als wij de waarde van x5 met 1 
laten toenemen. Dit resultaat hadden wij ook in paragraaf 2.3. 
reeds gevonden. 

Hoever kunnen wij nu х. maximaal laten toenemen? In plaats van 
1 eenheid x moeten wij opofferen } eenheid xj. Wij hebben in 
onze basisoplossing 60 eenheden xj. Wij kunnen dus maximaal 60 
eenheden xj opofferen. Het aantal eenheden x, dat wij maximaal 
in de oplossing kunnen introduceren is dus - voor wat bétreft 
de restrictie weergegeven іп de tweede rij - 


ор 120 


i 
Voor wat betreft de derde rij geldt het volgende: wil men x 
met 1 laten toenemen, dan zal x, moeten afnemen met 3/2. Men 
heeft in de basisoplossing thans x, = 80 (d.w.z. er zijn nog 
80 uur onbezet op machine 2). Men kan dus wat dat betreft x2 
maximaal laten toenemen tot 


89. 
375" 39 473 


Voor wat betreft de vierde rij vindt men op analoge wijze dat 
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men х2 maximaal Кап laten toenemen tot 
0 


T 40. 


Deze drie voorwaarden tesamen nemend concluderen wij dat wij 
х2 maximaal kunnen laten toenemen tot 40. De rekenregel is hier 
dus weer: bepaal de kleinste van de drie breuken 


T, m en Р . De kleinste van deze drie breuken is T Dit 


komt overeen met de rij van variabele xs. Variabele x5 gaat nu 

uit de basis om plaats te maken voor variabele x3. 

Wij transformeren nu het tableau van tabel 4.4. met als pivot- 

element het element іп de rij overeenkomend met хс en in de 

kolom overeenkomend met x5. Dit element hebben wij іп tabel 4.4. 

weer omcirkeld. De transformatie bestaat uit de volgende stap- 

pen: 

1. vermenigvuldig de vierde rij met 2 

2. tel de met 7} vermenigvuldigde nieuwe vierde rij op bij de 
eerste rij 

3. tel de met -і vermenigvuldigde nieuwe vierde rij op bij de 
tweede rij 

4. tel de met -3/2 vermenigvuldigde nieuwe vierde rij op bij de 
derde rij. 


Wij krijgen dan het tableau van tabel 4.5. 


15 3000 
40 


- tabel 4.5. Het derde simplex-tableau - 


Wij vragen ons thans weer af of wij de nu gevonden basisoplos- 

sing z = 3000, xj = 40, x> = 40, хз = 0, x, = 20, х5 = 0 nog 

kunnen verbeteren. Wij beschouwen daartoe weer de eerste rij: 
SÉ 15х3 + 15х5 = 3000 


ofwel 


š = 3000 = IIE; > 15х6 
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In de huidige basisoplossing is x3 = 0 еп xs = 0. Zouden wij of- 
wel x3 ofwel x5 een waarde # 0 geven, dan neemt de waarde van 

z dus niet toe maar áf, Het is nu niet meer mogelijk één van 

de niet-basis variabelen (x3 of x5) in de basis te introduceren 
teneinde daardoor z te doen toenemen. Dit betekent dat z = 3000, 
xj = 40, хә = 40, хз = 0, хц = 20 en xs = 0 de optimale oplos- 
sing is 7). 


De gehele Simplex-procedure kan met het opschrijven van de drie 
tableau's worden samengevat. Zie tabel 4.6. 


- tabel 4.6. Samenvatting simplex tableaus - 


4.3. SAMENVATTING VAN DE SIMPLEX-METHODE 


In de voorgaande paragraaf hebben wij de simplex-methode aan de 
hand van een voorbeeld toegelicht. Een lineair programmerings- 
model, dat uitsluitend < -restricties bevat en waarvan de doel- 


stellingsfunctie gemaximaliseerd moet worden, kan men als volgt 
schrijven: 


Maximeer 


CX] + C2 X> + ооа о т EnXn 
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onder de voorwaarden 


а ү 1х + ар2хо Р чом + а|аХи < b, 


I^ 
с 
N 


а21Х1 + а22Х2 "asa о ann 


IA 
су 


P P sen Т 
àm1*1 3n2*2 Зоо т 


X]; X»; ....»%9 ха Z 0 


Wij veronderstellen voorlopig, dat de getallen bj, bo, ...., bg 
alle niet negatief zijn. De simplex-methode kan dan als volgt 
worden samengevat: 

Stap 1: Voeg aan alle restricties een spelingsvariabele toe en 
schrijf de doelstellingsfunctie als 


во бук; e 53D * 23s T Cain = 0 


Schrijf nu het eerste simplextableau op. In dit tableau is 
(Z; Xotif Xp42» o; Kn+m) de basisoplossing (xq4], ....., 
Xn+m Zijn de spelingsvariabelen). 


Stap 2: Optimaliteitstest. Wij onderzoeken de optimaliteit van 
de basisoplossing door te kijken naar de elementen in de z-rij. 
Zijn deze alle > 0, dan is de basisoplossing behorende bij dit 
tableau de optimale oplossing. Indien één of meer elementen іп 
de z-rij < 0 zijn, gaan wij naar stap 3. 

Stap 3: Selectie van de variabele, die wij in de basis willen 
opnemen. Wij kiezen daarvoor die variabele waarvan het element 
in de z-rij het kleinst із 8). 

Stap 4: Selectie van de variabele, die wij uit de basis willen 
verwijderen. Wij bekijken daartoe de elementen in de kolom van 
de variabele, die wij in de basis willen opnemen. De elementen 
in deze kolom die > 0 zijn, delen wij op de overeenkomstige 
elementen in het rechterlid. Van de aldus gevonden breuken se- 
lecteren wij de kleinste. Het getal dat de noemer vormt van 
deze kleinste breuk is nu het pivot-element. De rij waarin het 
pivot-element voorkomt geeft de variabele aan die wij uit de 
basis gaan verwijderen. 

Stap 5: Veeg met het pivot-element de kolom van de variabele, 
die wij in de basis willen opnemen, schoon. Keer nu terug tot 
stap 2. 


Wij herhalen de stappen 2, 3, 4 en 5 totdat wij bij stap 2 con- 
cluderen, dat wij de optimale oplossing hebben gevonden. 


Een basisoplossing, waarvan de waarden van de basisvariabelen 
alle > 0 zijn, noemt men een nzet-gedegenereerde oplossing. 
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Een basisoplossing, waarvan de waarden van één of meer basis- 
variabelen - 0 zijn noemt men een gedegenereerde oplossing. In- 
dien geen gedegenereerde oplossingen voorkomen, leidt de simplex- 
methode in een eindig aantal stappen naar de optimale oplossing. 
Indien wel gedegenereerde oplossingen voorkomen is het mogelijk 
dat men na het uitvoeren van een bepaald aantal iteraties te- 
recht komt bij een basisoplossing, die men al eerder gevonden 
had. Men zegt dan, dat er sprake is van cycling. Er zijn tot 
dusverre geen gevallen bekend, waarin in praktijkproblemen 
cycling optrad. Overigens kan men door het aanbrengen van enkele 
simpele wijzigingen in de simplex-methode ervoor zorgen, dat 

men ook indien gedegenereerde oplossingen voorkomen de optimale 
oplossing in een eindig aantal stappen vindt: 


4.4. INTERPRETATIE VAN DE SIMPLEX-TABLEAUS 


Bij de bespreking van de simplex-methode in paragraaf 4.2. zijn 
wij reeds uitgebreid ingegaan op de betekenis van de elementen 
in de simplextableaus. Wij willen het daar besprokene thans kort 
samenvatten en beschouwen daartoe tabel 4.6. 

Kijken wij naar het tweede tableau in tabel 4.6. dan vinden wij: 

- in het rechterlid de waarden van de basisvariabelen; de oplos- 
sing behorende bij tableau 2 is dus: 2 = 2700, x; = 60, 

X. = 80, те = га 

- in де z-rij де reduced revenues (voor de echte variabelen) 
en de imputed costs (voor de spelingsvariabelen). De reduced 
revenue van produkt B bij deze oplossing is dus 7} (de reduced 
revenues hebben altijd een -teken іп de z-rij). De imputed 
cost van machine 1 is bij deze oplossing 22} (slack variable 
хз hoort bij machine 1); 

- in de kolommen van de niet-basisvariabelen, het aantal een- 
heden dat men van de basisvariabelen moet opofferen om één 
eenheid van de betreffende niet-basisvariabele in de oplossing 
te kunnen introduceren. Wil men bijvoorbeeld één eenheid хз 
in de oplossing introduceren (d.w.z. wil men bijvoorbeeld 1 
uur machine 1 vrijhouden voor bijvoorbeeld onderhoudswerkzaam- 
heden), dan moet men l eenheid van produkt A opgeven en krijgt 
men tevens } uur op machine 2 en 5 uur arbeid vrij. 


Het laatste simplex-tableau kan men volkomen analoog interpre- 
teren. De reduced revenues van alle echte variabelen zijn nu 
gelijk aan oul, De imputed cost van restrictie І is gelijk aan 
15 evenals de imputed cost van restrictie 3. De imputed cost 
van restrictie 1 is niets anders dan de schaduwprijs уап re- 
strictie 1. In paragraaf 3.1. concludeerden wij ook reeds, dat 
де schaduwprijs van restrictie 1 gelijk aan 15 was. 
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Het laatste simplex-tableau geeft dus behalve de optimale ор- 
lossing ook de schaduwprijzen van de restricties. Пе schaduw- 
prijs van restrictie 2 is gelijk aan nul, die van restrictie 3 
is gelijk aan 15. 


4.5. VRAAGSTUKKEN 


4-1. Gegeven is een onderneming, die in principe drie produkten 
(А, В еп С) zou kunnen maken. Van deze drie produkten zijn de 
volgende gegevens bekend: 


А В С 
Verkoopprijs 110 100 105 
Grondstofkosten 60 70 45 
Dekking I 50 30 60 
Loonkosten 60 20 40 
Dekking II -10 10 20 


Men wenst nu het produktiepakket te bepalen, dat leidt tot een 

zo hoog mogelijke winst. Daarbij zijn de volgende gegevens van 

belang: 

- de onderneming heeft thans 250 man direct produktief perso- 
neel; uit sociale overwegingen wenst men de werkgelegenheid 
op dit niveau van 250 man te handhaven; uitbreiding of inkrim- 
ping van het aantal personeelsleden is dus niet mogelijk; 

- er zijn 200 machines van type І; na aftrek van tijd nodig voor 
onderhoud zijn per dag beschikbaar 1400 uren voor machines 
van type І; | 

- ег zijn 200 machines уап type ІІ; па aftrek van tijd nodig 
voor onderhoud zijn per dag beschikbaar 1500 uren voor machi- 
nes van type II; 

- per dag zijn beschikbaar 2000 direct produktieve manuren; 

- het capaciteitsbeslag van de produkten is als volgt (uren per 
eenheid produkt): 


C 
manuren 3 1 2 
machines І 1 2 2 


machines II 2 2 1 
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- de loonkosten bedragen f 20,- рег manuur ofwel 7 160,- per 
mandag; 

| - de kosten van salarissen van niet-direct produktief personeel, 

de afschrijvingen, de rente en de overige vaste kosten bedra- 

gen in totaal f 10.000,- per dag. 


Vragen: 


a. Formuleer het LP-model en los het op. 


b. Hoe groot is het verlies van de onderneming per dag? Hoe 
groot is de schaduwprijs van arbeid? Kan men de winstgevend- 
heid van de onderneming verbeteren door 1 man extra in dienst 
te nemen of kan men, als men de winstgevendheid van de onder- 
neming zou willen verbeteren, beter overwegen | man te ont- 
slaan? Met welk bedrag neemt het verlies af bij deze beide 
potentiele maatregelen? 


4-2, a. Los het probleem van vraagstuk 2-1 op met behulp van de 
 simplex-methode (opmerking: de in vraagstuk 3-1 geïntrodu- 
ceerde verkooprestrictie blijft buiten beschouwing). 


b. Stel, dat een nieuw produkt C wordt ontwikkeld, waarvan de 
gegevens zijn: dekkingsbijdrage per eenheid f 250,-, сарасі- 
teitsbeslag per eenheid produkt 4 uren уап afdeling І en 2 
uren van afdeling II. Hoe luidt nu het optimale produktie- 
pakket? Beantwoord deze vraag door het probleem opnieuw - nu 
voor 3 produkten op te lossen. 


c. Met welk bedrag neemt пи de behaalde dekking af, als men ten- 
minste één eenheid van produkt B zou willen blijven produceren. 


4-3. a. Los het probleem van klokkenfabriek Frisia B.V. op met 
behulp van de simplex-methode (zie vraagstuk 27234 


b. Hoeveel uren van de houtafdeling zijn in het optimale pro- 
gramma niet benut? 
Hoe groot is de schaduwprijs van arbeid in de spuiterij? 


c. De stylist van Frisia B.V. ontwerpt een nieuw soort stoel- 

 tjesklok. Deze is technisch geheel gelijk aan de bestaande 
stoeltjesklok en vraagt precies evenveel arbeid als de be- 
staande stoeltjesklok in de spuiterij en in de montage. Qua 
houtbewerking is de stoeltjesklok nieuwe stijl echter aanzien- 
lijk bewerkelijker: de nieuwe stoeltjesklok vraagt 2 uur ar- 

 beid van de houtafdeling. De brutomarge van de nieuwe klok 
bedraagt 240,- per klok. 
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Wordt deze nieuwe klok іп het optimale produktieprogramma 
opgenomen? Beantwoord deze vraag zonder het probleem opnieuw 
(met 3 produkten) op te lossen (zie de opmerking bij de op- 
lossing van vraagstuk 4-2). 


d. Los het probleem van Frisia B.V. nu op voor 3 produkten met 
behulp van de simplex-methode, 


4-4. а. Frisia B.V. heeft nu naast de 3 bestaande produkten ( 


(friese staartklok, friese stoeltjesklok en friese stoeltjes- 


klok nieuwe stijl) de mogelijkheid de originele uurwerken 
die Frisia B.V. zelf maakt, te verkopen aan andere fabrikan- 
ten van stijlklokken. In dat geval bedraagt de brutomarge 
f 60,- per verkocht uurwerk. 

Los het probleem van Frisia B.V. nu met 4 produkten op met 
behulp van de simplex-methode. 


b. Wat is nu de schaduwprijs van arbeid in de spuiterij? 


4-5. Gegeven is een boer, die in principe drie soorten gewas- 
sen kan verbouwen; deze gewassen worden aangeduid met de let- 
ters А, В еп С. De boer bezit 60 hectare bebouwbare grond. 

Eén ha grond ingezaaid met gewas A levert 2 ton verkoopbaar 
produkt van gewas A ор; voor de gewassen B еп С is dit respec- 
tievelijk 1,8 en 1,6 ton. De verkoopprijzen bedragen (in gul- 
dens per ton): 


A 800,- 
4715, - 
б, 180; 


De variabele kosten (zaaigoed, brandstof voor landbouwwerktui- 
gen, etc.) bedragen 7 100,- per ha ingezaaid met gewas A, 

f 55,- per ha ingezaaid met gewas B en f 35,- per ha ingezaaid 
met gewas C. 

Gedurende de maanden april tot en met september moet de boer 
20 uur arbeid verrichten voor iedere hectare ingezaaid met 
gewas A, 10 uur arbeid voor iedere hectare ingezaaid met 

gewas В en 15 uur arbeid voor iedere hectare ingezaaid met 
gewas C. 

De boer moet thans beslissen welke hoeveelheden hij van de 
drie gewassen in het komende jaar gaat verbouwen. Hij wenst in 
de maanden april tot en met september ten hoogste 1000 uur te 
werken en binnen die restrictie zijn bruto-winst te maximeren. 


a. Formuleer bovenstaand probleem als een LP-probleem. 
b. Los het probleem op met behulp van de simplex-methode. 
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4-6. Gegeven is een chemische fabriek: bij het produktieproces 
van deze fabriek komen twee afvalstoffen (aan te duiden met A en 
B) vrij. Deze afvalstoffen worden geloosd op het oppervlaktewater. 
Aan het lozen van deze afvalstoffen zijn tot dusverre geen kosten 
verbonden. 

Het lozen van de afvalstoffen brengt een aanzienlijke verontrei- 
niging van het oppervlaktewater teweeg. De engineering-afdeling 
heeft nu een ontwerp gemaakt voor een installatie, waarmee per 
dag 100 kg van stof A en 200 kg van stof B aan het afvalwater 

kan worden onttrokken; met behulp van deze installatie kunnen de 
stoffen A en B worden verwerkt tot halffabrikaten voor andere 
chemische industrieën. In principe kunnen vier soorten halffabri- 
katen worden geproduceerd (aan te duiden met 1, 2, 3 en 4). Enke- 
le gegevens met betrekking tot deze halffabrikaten zijn: 


1 2 3 4 
Verkoopprijs per eenheid ід зе з 32,90 AO TR 
Grondstofkosten (excl. А еп В) мж. «Р жо ñ y> 5,77 
Grondstofkosten A p.m. p.m. p.m. p.m. 
Grondstofkosten B p.m. p.m. p.m. p.m. 


Voor het produceren van één eenheid van de halffabrikaten 1, 2, 
3 en 4 zijn nodig (kg): 


1 4 
van stof А 2 2 0 1 
van stof B 1 5 2 1 


De vaste kosten van de nieuwe installatie worden geraamd op 
f 300.000,-- per jaar ofwel (250 werkdagen per 18481) 7. 1.200,*— 
per dag. 


Als men besluit deze nieuwe installatie in gebruik te nemen, wel- 
ke halffabrikaten moet men dan gaan produceren als men streeft 
naar een zo groot mogelijke winst (c.q. naar een zo klein moge- 
lijk verlies)? 

a. Formuleer dit probleem als een LP-probleem. 

b. Los het probleem op met behulp van de Simplex-methode. 

c. Is het - als men puur bedrijfseconomisch redeneert - aantrekke- 
lijk om deze nieuwe installatie in gebruik te nemen? 

d. Op de engineering-afdeling krijgt iemand nog een nieuw idee, 
namelijk om een halffabrikaat 5 te gaan produceren met behulp 
van de nieuwe installatie. De verkoopprijs van halffabrikaat 5 
bedraagt f 20,-- per eenheid, de grondstofkosten (excl. de 
grondstofkosten voor А en B) f 13,-- per eenheid. Voor het 
produceren van één eenheid van halffabrikaat 5 is nodig 1 kg 
van stof A en 2 kg van stof B. Verandert door de aanwezigheid 
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van deze nieuwe mogelijkheid de samenstelling van het optimale 
produktiepakket, zoals berekend bij vraag b? Beantwoord deze 
vraag zonder het probleem opnieuw op te lossen (zie de opmer- 
king bij de oplossing van vraagstuk 4-2). 


4-7. De in vraagstuk 4-6 bedoelde chemische fabriek moet nu voor 

het lozen van de afvalstoffen een verontreinigingsheffing betalen. 

Deze heffing bedraagt f 5,-- per kg voor stof A en f 3,-- рег kg 

voor stof В. : 

a. Formuleer nu opnieuw het probleem van de chemische fabriek als 
LP-probleem. 

b. Los het probleem op met behulp van de Simplex-methode. 


S. HET GEBRUIK VAN KUNSTMATIGE 


EN SURPLUS-VARIABELEN 


5.1. INLEIDING 


Tot dusverre hebben wij uitsluitend problemen beschouwd, die 
konden worden geformuleerd als: 
maximeer 


C1X]| + с2Х2 ғ оооо о + Coin 
onder de voorwaarden: 


8а11Х1 + а12Х2 deli ate. d + ах 


IA 
с 
mer 


a5j]X] + а22Х2 + seses + ах £ b» 


amx] + 8m2X2 %..... + AmnXn 5 bm 
waarbij de getallen 

bis: bost soris By SH 
zijn. 


Wij beschouwen nu het volgende probleem. Een veehouder heeft de 
beschikking over twee soorten veekoeken die wij A en B noemen. 
Een kilogram veekoeken van type А bevat 1 eenheid mineralen en 
1 eenheid vitaminen. Een kilogram veekoeken van type B bevat | 
eenheid mineralen en 2 eenheden vitaminen. Een kilogram A kost 
f 3,-, Eén kilogram B kost f 5,-. De veehouder wenst voor zijn 
veestapel een menu te kiezen, dat per dag minstens 6 eenheden 
mineralen en 8 eenheden vitaminen bevat en dat tevens zo goed- 
koop mogelijk is. Wat is de optimale samenstelling van dit menu? 
Wij herkennen hierin onmiddellijk een lineair programmerings- 
probleem, namelijk: 

minimeer 


Зх: * 3X9 


onder de voorwaarden 
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IV 
о 


Ket: хо (mineralen) 


V 
оо 


ki +. 2х, 2 (vitaminen) 


Indien wij nu proberen, dit probleem in bovengenoemde vorm te 
schrijven krijgen wij: 
maximeer 


"Эк: * 5X» 

onder de voorwaarden: 
"oo S 79 
ехе 2X5. 5 76 


Xj, 4$, E 9 


Wij zijn er nu wel in geslaagd de » -restricties om te zetten 
in < -restricties, maar wij hebben nu negatieve getallen in het 
rechterlid. Dit houdt in dat wij, als wij de in paragraaf 4.3. 
beschreven oplossingsmethode toepassen, geconfronteerd worden 
met een simplex-tableau waarvan de basisoplossing niet feasible 
is (wij moeten dan immers een spelingsvariabele toevoegen aan 
iedere restrictie; de basisoplossing van het eerste tableau be- 
staat dan uit deze spelingsvariabelen; deze spelingsvariabelen 
hebben dan negatieve waarden, hetgeen overeenkomst met een in- 
feasible solution). 

In het algemeen kan men stellen, dat men iedere > -restrictie 
kan transformeren in een < -restrictie door te vermenigvuldigen 
met -І. Мі) spreken nu af, dat wij bij het toepassen van de 
simplex-methode iedere restrictie steeds zo zullen formuleren, 
dat de constante in het rechterlid > 0 is. Het bovengenoemde 
probleem (het zgn. dieetprobleem), formuleren wij dus als: 
minimeer ? 


Зх] » 5X2 
onder de voorwaarden 
Kit х2 2 6 


8 


IV 


X] + 2х2 


ху, хә 2 0 
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5.2. DE INTRODUCTIE VAN SURPLUS-VARIABELEN 
Wij kunnen naar analogie van de spelingsvariabelen nu in de 
restricties surplus-variabelen introduceren. Wij krijgen dan 


het volgende probleem: 


бт 3x] "S 5х2 = 0 


XT T ao Ra 


хі + 2x2 T XH 


waar x3 voorstelt het aantal eenheden mineralen dat men méér 
geeft dan strikt nodig is. Aangezien wij weer eisen хз > 0 be- 
tekent xj + x) - хз = 6 hetzelfde als xj + x; 2 6. Aan xy kan 
men een analoge interpretatie geven. Het eerste simplex-tableau 
wordt nu (zie tabel 5.1): 


2 X] X2 X3 XL RL 
1 "3 zë 0 0 0 
0 1 1 wr 0 

0 l 2 0 Si 8 


- tabel 5.1. simplex-tableau voor dieetprobleem met surplus 
variabelen - 


Het probleem is nu, dat dit tableau geen basis-oplossing bevat. 
De kolommen van x3 en х, bevatten namelijk het element -1. Wil 
men dit getal -1 іп de tweede rij transformeren in +1 door de 
hele rij te vermenigvuldigen met -1, dan wordt het rechterlid 
van +6 ook omgezet in -6. Ор deze wijze vinden wij geen basic 
feasible solution, die als startoplossing voor de simplex-me- 
thode kan fungeren. 


5.3. DE INTRODUCTIE VAN KUNSTMATIGE VARIABELEN 

Om een basic feasible solution te vinden die als startoplossing 
voor de simplex-methode kan fungeren introduceren wij naast de 
surplus-variabelen tevens kunstmatige variabelen (artificial 


variables): 


Xi: 722” X4 +х5 = 6 


II 
oo 


хі + 2х2 TX * X6 
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Deze kunstmatige variabelen worden alleen geïntroduceerd om een 
basic feasible solution te vinden. Wij willen de kunstmatige 
variabelen zo snel mogelijk elimineren en schrijven daarom de 
doelstellingsfunctie als: 

minimeer 


3x1 * 5х2 + Мх 5 + MXg 


waar M een zeer groot positief getal is (zeg bijvoorbeeld 
M = 1.000.000). 
Wij hebben dus nu het volgende probleem: 


ECCO) 5. - Mxg - Mxg = 0 
Ei * 2х2 "Же Ф х5 = 6 
Жү: * 2х2 7 хц жа” 8 


Het eerste simplex-tableau wordt nu: 


- tabel 5.2. Simplex-tableau voor dieetproble 
variabelen en kunstmatige variabelen - 


Dit simplex-tableau bevat geen basic feasible solution. Door nu 


de met M vermenigvuldigde tweede rij en de met M vermenigvuldig- 
de derde rij bij de eerste rij op te tellen krijgen wij: 


se Ak AN -5%3М | 


- tabel 5.3. Simplex-tableau voor dieetprobleem met basic 
feasible solution - 


Een iets eenvoudiger notatie verkrijgen wij door het tableau 
van tabel 5.3. te schrijven als in tabel 5.4. 


- tabel 5.4. Simplex-tableau voor dieetprobleem met basic 
feasible solution - 


Men moet het tableau van tabel 5.4. lezen net als dat van tabel 
5.3. In de x,-kolom іп de z-rij staat dus (-3%2М) en niet -3. 
De notatie van tabel 5.4. is gekozen omdat deze overzichtelijker 
is dan die van tabel 5.3. 


5.4. TOEPASSING VAN DE SIMPLEX METHODE OP HET DIEET- 
PROBLEEM | 


Wij kunnen nu op het tableau van tabel 5.4. de simplex-methode 
toepassen. Wij bepalen eerst welke variabele wij in de basis 
willen opnemen. Wij wensen nu de doelstellingsfunctie te mini- 
meren. Daarom kiezen wij die variabele, waarvan de coefficient 
in de z-rij het grootst is. Voor xj is deze coefficient -3+ 2M, 
voor x, is deze coefficient -5 + ЗМ. Aangezien M zeer groot is, 
is -5 + ЗМ > -3%2М. Wij kiezen dus xj . Nu bepalen wij welke 
variabele wij uit de basis gaan verwijderen. Wij berekenen daar- 
toe 6/1 en 8/2. De kleinste van deze twee breuken is 8/2. Daar- 
mee hebben wij het pivot-element gevonden (wij hebben het pivot- 
element in tabel 5.4. omcirkeld). Variabele xg gaan wij uit de 
basis verwijderen. Wij gaan nu de kolom van x2 schoonvegen met 
het pivot-element. Hiertoe delen wij eerst de xg-rij door 2. 


- tabel 5.5. Tweede simplex-tableau voor dieetprobleem - 
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Zo vinden wij de laatste rij van het tweede tableau (zie tabel 
5.5.). Om nu іп de z-rij een 0 te krijgen moeten wij de met 

(5- 3M) vermenigvuldigde laatste rij van tabel 5.5. optellen bij 
de 2- rij van tabel 5.4, Wij vinden dan bijvoorbeeld voor de co- 
efficient van xj in de z-rij: 


73 42M LAM) * 4 = -1 + 1M 


en dit hebben wij opgeschreven in de z-rij van tabel 5.5. Veel 
eenvoudiger en volstrekt gelijkwaardig is de volgende methode: 
tel de met 5 vermenigvuldigde laatste rij op bij de eerste rij 
en tel dan de met -3M vermenigvuldigde laatste rij op bij de 
tweede rij. Aangezien wij in de notatie van tabel 5.4. de let- 
ter M weglaten, kunnen wij het tweede gedeelte van bovenstaande 
zin ook vervangen door: tel de met -3 vermenigvuldigde laatste 
rij op bij de tweede rij. Wij kunnen dus bij het schoonvegen 
van de kolommen de z-rij behandelen als twee afzonderlijke rijen 
die wij ieder voor zich schoonvegen. Men noemt de tweede rij 
ook wel de hulpdoelstellings functie. 

In tabel 5.5. komt пор één kunstmatige variabele (nl. xs) voor 
in de basisoplossing. De waarde van de doelstellingsfunctie is 
20 + 2М, dus een zéér groot getal. Pas als de waarde van de 
doelstellingsfunctie een getal is, waarin M niet meer voorkomt, 
hebben wij een feasible solution van het oorspronkelijke pro- 
bleem zoals geformuleerd іп paragraaf 5.1. Dit is het geval dan 
en slechts dan als in de rij van de hulpdoelstellingsfunctie 

in het rechterlid een 0 voorkomt. 

Wij voeren nu nog een iteratie uit, Wij selecteren xj als vari- 
abele die wij in de basis willen introduceren (wij kiezen uit 
xj en xi met als coefficienten -} + ІМ en -5/2 + ІМ, Aangezien 
-4 + JM > - 5/2 + ІМ kiezen wij xj). De variabele s verlaat nu 
de basis en na schoonvegen krijgen wij tabel 5.6. 


- tabel 5.6. Derde simplex-tableau voor dieetprobleem - 


In dit tableau bestaat de gehele tweede rij uit nullen, behalve 
de elementen in de kolommen van de kunstmatige variabelen. Dit 
is geen toeval; laat men de kolommen van de kunstmatige varia- 
belen buiten beschouwing, dan is immers in tabel 5.4. de tweede 
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rij al een lineaire combinatie уап de derde en vierde rij (па- 
melijk eenvoudig de som van de derde en vierde rij). Indien men 
erin slaagt in het rechterlid іп de tweede rij een 0 te verkrij- 
gen, zal dus de gehele tweede rij met uitzondering van de ko- 
lommen der kunstmatige variabelen uit nullen bestaan. 

Wij hebben nu een basis-oplossing waarin de kunstmatige variabe- 
len niet meer voorkomen. Wij passen nu de optimaliteitstest 

toe. 


Aangezien wij thans de waarde van de doelstellingsfunctie mini- 
meren, hebben wij de optimale oplossing gevonden indien géén 
der elementen in de z-rij > 0 is. Alle elementen in de z-rij 
met uitzondering van de elementen in de kolommen der kunstmatige 
variabelen zijn < 0. Aangezien wij de kunstmatige variabelen 
buiten beschouwing laten, hebben wij in tabel 5.6. de optimale 
oplossing gevonden. Deze luidt: 2 = 22, xj = á, X? = 2, хз = 0. 
en x, = 0. 

Wij vatten de gehele procedure nog eens samen in tabel 5.7. 

In tabel 5.7. hebben wij als tableau 0 het tableau van tabel 
5.2. geschreven. Dit tableau bevat nog geen basisoplossing. Wel 
weten wij dat wij van de kunstmatige variabelen gemakkelijk 
basisvariabelen kunnen maken, hetgeen wij aangeven door deze 
variabelen tussen haakjes te plaatsen in de kolom "Basis". 


Tableau | Basis 


- tabel 5.7. Simplex-procedure voor dieetprobleem - 


5.5. VOORBEELD VAN EEN PROBLEEM MET EEN - -RESTRICTIE 


Tot besluit van dit hoofdstuk geven wij nog een voorbeeld van 
een probleem met ёёп < -restrictie, één > -restrictie еп één 
= =restrictie, namelijk: 

maximeer 


lOx; + 12x9 

onder de voorwaarden 
2x] + 4x2 < 60 
Xi t X 25 
2хі + 2х2 > 40 


X], X2 2 0 


Wij schrijven dit probleem als: 


maximeer 
10x] + 12х; -Мх, -Mxg 
onder de voorwaarden 
2х] + 4х2 + хз - 60 
кі» Xo + Xi = 25 
2ху + 2x» - xs + xç = 40 


Aan de = -restrictie voegen wij dus uitsluitend een kunstmatige 
variabele toe, aan de > -restrictie een surplus variabele èn 
een kunstmatige variabele. Wij willen de doelstellingsfunctie 
maximeren; om nu zo snel mogelijk de kunstmatige variabelen 

te elimineren hebben wij de kunstmatige variabelen in de doel- 
stellingsfunctie de coefficient -M gegeven, waarbij M een groot 
positief getal is, De oplossing van dit probleem is gegeven in 
tabel 5.8. 
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Basis 


N 


Tableau 


—- tabel 5,8 e 
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5.6. VRAAGSTUKKEN 


5-1. Een fabrikant van haarwaters wenst een nieuw soort haar- 
water aan het bestaande assortiment toe te voegen, Het produk- 
tieproces bestaat uit twee fasen. In de eerste fase worden een 
drietal standaardmengsels А, В еп С gemaakt. In de tweede fase 
worden deze standaardmengsels in een bepaalde verhouding gemengd 
tot het eindprodukt. 

De nieuwe soort haarwater moet zowel anti-roos als anti-haaruit- 
val bestanddelen bevatten en wel tenminste 40 mg anti-roos be- | 
standdeel en tenminste 60 mg anti-haaruitval bestanddeel per 

100 ml. 

De standaardmengsels bevatten per ml de volgende hoeveelheden 
anti-roos en anti-haaruitval bestanddelen: 


А В 6 
anti roos 1 1 2 
anti haaruitval 1 0 1 


De kostprijzen van de standaardmengsels А, B еп С bedragen res- 
pectievelijk 30, 20 en 50 cent per 100 ml. 


De fabrikant wenst het nieuwe type haarwater op zodanige wijze 
samen te stellen uit de drie standaardmengsels, dat de grond- 
stofkosten van het nieuwe mengsel per 100 ml zo gering mogelijk 
zijn. 


а. Formuleer dit probleem als een LP-probleem en los het op met 
behulp van de simplex-methode. 


b. Als men als eis stelt, dat het nieuwe type haarwater tenminste 
61 mg anti-haárui tval-middel moet bevatten, stijgen dan de 
grondstofkosten van het nieuwe produkt? Zo ja, met welk be- 
drag рег 100 пі? 


с. Los het probleem opnieuw op met behulp van de simplex-methode 
nu met een kostprijs voor standaardmengsel C van 10 cent in 
plaats van 50 cent per liter. 


5-2. Een veehouder beschikt over 4 soorten veekoeken А, В, С еп 
D. De vitaminegehalten van de й soorten veekoeken zijn vermeld 
in onderstaande tabel (eenheden per stuk): 


Veekoek: А В С D 
Vitamine VI 4 0 2 1 
V2 2 1 1 0 
УЗ 0 0 2 2 
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Deze veehouder wenst een zodanig menu samen te stellen, dat 
iedere koe per dag tenminste 50 eenheden vitamine Vl, 40 een- 
heden vitamine V2 en 10 eenheden vitamine V3 krijgt. De inkoop- 
prijs van de veekoeken is: 


А f 4,- per stuk 
B 7.1,--рек stuk 
С FI,“ DOC stuk 
D f 4,- per stuk 


De veehouder wenst tevens een zo goedkoop mogelijk menu samen 
te stellen. 


Formuleer bovenstaand probleem als een LP-probleem en los het 
op met behulp van de simplex-methode. 


5-3. Los het model van vraagstuk 2-7 op met behulp van de 
simplex-methode. (Aanwijzing: laat ter beperking van het reken- 
werk de overbodige restricties weg.) 


5-4. Los het probleem van vraagstuk 2-8 op met behulp van de 
simplex-methode. 


5-5. Los het probleem van de oliemaatschappij van vraagstuk 2-9 
op met behulp van de simplex-mehtode. 


6. HET DUALE PROBLEEM 


6.1. DE ZIENSWIJZE VAN EEN BUITENSTAANDER 


Wij bekijken thans opnieuw het probleem van de optimale samen- 
stelling van het produktiepakket uit paragraaf 2.2.: 
maximeer 


45хі + 30х; 
onder de voorwaarden 
Aki 22 Ro 5 840 Масһіпе 1 
хі + 2x2 < 140 Масһіпе 2 
кі Є у» $ 90 Arbeid 
X], хә 2 0 


Wij beschouwen dit probleem nu echter vanuit het gezichtspunt 
van een buitenstaander, die overweegt alle produktiecapacitei- 
ten waarover de onderneming beschikt te huren. Deze buitenstaan- 
der vraagt zich nu af welke huurprijzen hij zal moeten betalen. 
Hij streeft er enerzijds naar in totaal een zo klein mogelijke 
som te betalen voor het huren van alle produktiecapaciteiten, 
terwijl anderzijds de door hem geboden huurprijzen voor de on- 
dernemer acceptabel moeten zijn. Stel deze buitenstaander is 
bereid u, gulden per uur te betalen voor machine l, ua gulden 
рег uur voor machine 2 еп из gulden per uur voor de arbeid. In 
totaal zal deze buitenstaander voor alle produktiecapaciteiten 
tesamen dus moeten betalen: 


12041 + 1400 + 80из 


De ondernemer, die een dergelijk bod ontvangt, zal nu als volgt 
redeneren: de dekking van produkt A bedraagt 45 gulden per een- 
heid, Produkt А vraagt 2 uren van machine 1, 1 uur van machine 

2 еп І uur arbeid, die bij verhuur 241 + u, + из gulden opbren- 
gen. Het in produktie nemen van produkt A blijft dus aantrekke- 


72 


lijk zolang 45 > 20) + uj + из. Zolang deze restrictie geldt 
zal de ondernemer niet tot verhuur overgaan. De buitenstaander 
moet dus zodanige prijzen uj, uy en оиз bieden, dat geldt: 


201 + uz + ua 2 45 
Evenzeer moet gelden: 
üi + 2U) fus 2 39 


(immers indien dit niet zou gelden is het voor de ondernemer 
aantrekkelijker om produkt B te maken, dan om capaciteiten te 
gaan verhuren). De buitenstaander wordt dus geconfronteerd met 
het volgende probleem: 

minimeer 


120u, + 140ч2 + 80u3 

onder de voorwaarden 
2u; + йо * Mg e 43 ці» ورلا‎ uz 2 0 
uj + 24274-93 € 30 


6.2. HET PRIMALE EN HET DUALE PROBLEEM 


Men kan ieder lineair programmeringsprobleem schrijven in de 
vorm: 
maximeer 


C1X| + CX Фо аж З Coin 


onder де voorwaarden 


IA 


a]]X] + 812X2 + ..... + аху < bI 


a2]X] + 822X2 + evene + Ann < b> 


àm|X] T Xm2X2 + ۰۰۰۰۰ + amnXn < bm 


Deze vorm noemt men de primale vorm van het lineair programme- 
ringsprobleem. (Het is mogelijk dat in de primale vorm negatieve 
constanten in het rechterlid voorkomen.) Bij ieder lineair pro- 
grammeringsprobleem behoort het zogeheten duale probleem. Dit 
duale probleem kunnen wij als volgt formuleren. Bij iedere re- 
strictie hoort een duale variabele; er zijn dus m duale variabe- 
len uj, l2, ....., ug. De doelstellingsfunctie van het duale 
probleem is bu) + b2u2 + ..... + bgum. 

Deze doelstellingsfunctie moet worden geminimeerd. Bij iedere 
variabele van het primale probleem behoort een restrictie, er 
zijn dus n restricties en deze hebben de volgende vorm: 


IV 


"1 


IV 


"2 


ajnu] + a2nu2 + ..... + AmnUm È Cn 
Tenslotte moet ook gelden: 


ці, UZ, ....., Um 2 0 


Het probleem van de buitenstaander in paragraaf 6.1. is dus het 
duale probleem van de ondernemer uit paragraaf 2.2, 


6.3. EEN NUMERIEK VOORBEELD 


Ter illustratie van het in de voorgaande paragraaf besprokene 
geven wij hier een numeriek voorbeeld. In paragraaf 5,5. ont- 
moetten wij het volgende probleem: 

maximeer 


10x] + 12x9 
onder de voorwaarden 

2x] + 4x, < 60 

хі ко. 25 

2x; + 245 > 40 

Xi któ 2 U 
Om van dit probleem het duale probleem te kunnen vinden moeten 
we het eerst in de primale vorm schrijven. Van de > -restrictie 
maken wij gemakkelijk een < -restrictie, n.l. door te vermenig- 


vuldigen met -l. De = -restrictie schrijven wij als twee restric- 
EIER 5.1. 


Xr * xw 5.25 
en 
Zi + хо 2 £5 
Hiervan vermenigvuldigen wij de laatste restrictie met -1 ten- 


einde de primale vorm te verkrijgen. Wij vinden aldus 
maximeer 


10x] + 12x5 


onder de voorwaarden 
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2x| + Axa < 60 
4*9 4» S. 29 
“Қа 5 4.5 "42 
"2*4 «фо S -40 


Het duale probleem is: 
minimeer 
60uj + 25u, = 25u3 = 40u, 


onder de voorwaarden 


IV 


10 
12 


дит $ ú, * ца > Ruh 


IV 


чи; ^ Чо ^ ug * Жы 


Has HA EU‏ ورلا وولا 


6.4. HET OPLOSSEN VAN HET DUALE PROBLEEM 


Wij lossen thans het in paragraaf 6.1. geformuleerde duale pro- 
bleem op. Wij schrijven dit probleem als: 


minimeer 
12041 + 1400, + 80из + Mug + Muz 
onder de voorwaarden 
2uj + цо * ^ $3 Ui * ME - 45 
312057 2ы, Фо Mà SM + u7 є 30 


De oplossing van het duale probleem lezen wij af uit het laatste 
tableau van tabel 6.1.: uj = 15, ug = 0, ug = 15. De buiten- 
staander moet dus 15 gulden per capaciteitseenheid van machine 
1, 0 gulden per capaciteitseenheid voor machine 2 en 15 gulden 
per uur arbeid bieden. Deze prijzen zijn de schaduwprijzen van 
de restricties van het oorspronkelijke probleem. Het minimale 
bod, dat de buitenstaander met succes kan uitbrengen is 3000, 
hetgeen precies gelijk is aan de maximale winst van de onder- 
nemer. 

De z-rij van het laatste tableau kan men schrijven als: 


z = 3000 + 2002 + 40ui, + 4005 


NI— NI— NI— 
NI— Nol NI— 


20/3-100/3 -160/3 
0 0 0 

1/3 -2/3 1/3 
1/3 -2/3 


- tabel 6.1. Het optimale produktiepakket. 
Oplossing van het duale probleem - 


Wil men in de oplossing uy = 1 maken dan zal dus de waarde van 
Z toenemen met 20. Men kan dit als volgt interpreteren: indien 
de buitenstaander 1 gulden per capaciteitseenheid van machine 

2 gaat bieden, dan moet hij zijn totale bod met 20 gulden ver- 
hogen om het nog juist acceptabel te doen zijn. Dit komt omdat 
de ondernemer nog 20 ongebruikte eenheden van machine 2 heeft. 
Wil men in de oplossing u, = 1 maken, dan zal de waarde van 

2 toenemen met 40. Men kan dit als volgt interpreteren: indien 
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de buitenstaander aan surplusvariabele u, de waarde 1 wil geven, 
betekent dit dat in feite de eerste restrictie wordt 


241 + us + из 2 46 


De buitenstaander moet dan zijn totale bod met 40 verhogen, om- 
dat de opbrengst voor de ondernemer van produkt А met 1 gulden 
toeneemt en de ondernemer in zijn optimale programma 40 eenheden 
А produceert. 

Het is dus geen toeval, dat in de z-rij van het laatste tableau 
van het duale probleem de waarden van de basis-variabelen van 

de optimale oplossing van het primale probleem staan. 

Zonder bewijs geven wij de volgende stelling: 

Als het primale probleem een eindige optimale oplossing heeft, 
heeft het duale probleem eveneens een eindige optimale oplossing. 
De optimale waarde van de doelstellingsfunctie van het primale 
probleem is dan gelijk aan de optimale waarde van de doelstel- 
lingsfunctie van het duale probleem. 


6.5. VRAAGSTUKKEN 


6-1. a. Formuleer het duale probleem van het in vraagstuk 2-1 
geformuleerde probleem. 
Welke betekenis kan men aan de duale variabelen toekennen? 
Op welke wijze kan men de restricties van het duale probleem 
interpreteren? Hoe kan men de doelstellingsfunctie van het 
duale probleem interpreteren? 


b. Los het duale probleem op met behulp van de і об-е. 


6-2. а. Formuleer het duale probleem van het in paragraaf 5.1. 
gegeven dieet probleem, 


b. Los het duale probleem op met behulp van de simplex-methode. 
6-3. a. Formuleer het duale probleem van het in vraagstuk 5-2 


gegeven dieetprobleem. 
b. Los het duale probleem op met behulp van de simplex-methode. 


6-4. Gegeven is het volgende probleem: | 
minimeer 


10x] + 14x» * 5X3 
onder de voorwaarden 
хет” Mit 2x3 2 16 
dE X4 * Зхз 2 12 
Xi о x. * ах; - 20 
X4, Хә, X3 20 


a. Schrijf dit probleem in de primale vorm. 
b. Formuleer het duale probleem. 


6-5. Beschouw opnieuw het probleem van klokkenfabriek Frisia, 
zoals geformuleerd in vraagstuk 2-2 en schrijf het bij dit vraag- 
stuk behorende LP-model op. 

a. Voeg aan iedere restrictie een duale variabele toe. Geef bij 
elke duale variabele aan welke (economische) betekenis men 
aan deze variabele kan toekennen. 

b. Formuleer het duale probleem. 

C. Los het duale probleem op met behulp van de simplex-methode. 


6-6. Beschouw opnieuw de bouwonderneming van vraagstuk 2-7 en 

schrijf het bij dit vraagstuk behorende LP-model op. 

a. Voeg aan iedere restrictie een duale variabele toe. Geef bij 
elke duale variabele afzonderlijk aan welke (economische) 
betekenis men aan deze variabele kan toekennen. 

b. Schrijf het probleem in de primale vorm. 

с. Formuleer het duale probleem. 


7. GEVOELIGHEIDSANALYSE 


71. GEVOELIGHEIDSANALYSE VAN DE COEFFICIENTEN VAN DE 
DOELSTELLINGSFUNCTIE 


In paragraaf 3.3. hebben wij aan de hand van de grafische oplos- 
singsmethode reeds laten zien hoe men de gevoeligheid van de 
gevonden oplossing voor een verandering in één der coefficienten 
van de doelstellingsfunctie kan onderzoeken. Wij willen thans 
laten zien hoe men bij een groter probleem (d.w.z. een probleem 
met meer dan twee beslissingsvariabelen) deze gevoeligheidsana- 
lyse uitvoert. 


Wij beschouwen het volgende probleem: 
maximeer 


8хі + 6x2 + 5хз + 4x, 
onder de voorwaarden 
хі + 2x2 + хз * хь 5 60 І X4,X5,X4,X, 2 0 
2x; + 2Х2 + хз t 3x, 5 100 kt 


De oplossing van dit probleem is gegeven in tabel 7.1. Voor zéér 
kleine wijzigingen in cj (bijvoorbeeld c; - 7,99 іп plaats van 
cj = 8) zal de optimale oplossing vermoedelijk niet gevoelig 
zijn. Preciezer gezegd: de bij cj = 8 behorende optimale oplos- 
sing xj = 40, x; = 0, хз = 20 en x, = 0 is vermoedelijk ook de 
optimale oplossing als cj - 7,99. 


Wij gaan nu de coëfficiënten van de doelstellingsfunctie één 
voor één varieren en vragen ons af binnen welke grenzen deze 
coëfficiënten kunnen varieren zonder dat de optimale oplossing 
verandert. Als wij dus cj variëren, nemen wij aan dat сә, сз 

en с, ongewijzigd blijven en zoeken dan de getallen а en Ь мааг- 
voor geldt: voor iedere cj die voldoet aan а < cj < Б, geldt 

dat de bij c = 8 behorende optimale oplossing (xj = 40, x2 = 0, 
хз = 20, х, = 0) de optimale oplossing 15. 


Voor c5, сз en сі zoeken wij soortgelijke intervallen. 


^ CADRE: 74157 


Wij merken op dat wij kunnen schrijven voor de waarde уап 2 in 
de optimale oplossing: 


420 = 8 * 40 + 5 ж 20 


Men zegt пи: де z-rij in het laatste tableau is een lineaire 
combinatie van de z-rij in het eerste simplex-tableau en de 
tweede en derde rij in het laatste simplex-tableau. Wij kunnen 


schrijven: 
20 *iB-*-0 +. S w 2 
яті *8*72*45 3 (*]) 
A O + бе (-1) + 5 « 2 
З. 0*8*]^*5 x» (-1) 
Ofwel, als wij іп plaats van de getallen 8, 6, 5, 4 schrijven 
С|» C2, C3 en сі: 
CR MERE я Ü + c. «2 
eten + сүзге есе (SE) (А) 
2 з 0 + с) * (21) + ca * 2 
3= 0 + Ci * l + C3 * (7-1) 
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In paragraaf 4.4. zagen wij reeds dat in de z-rij van het laats- 
te tableau voor de kolommen van de echte variabelen de reduced 
revenues voorzien van een -teken verschijnen. De reduced revenue 
van produkt 2 is hier dus -4 en dit betekent, dat met de huidige 
waarden van сі, сэ, сз en сц introductie van één eenheid van 
produkt 2 de waarde van de doelstellingsfunctie met 4 zou doen 
dalen. Men kan dit ook inzien als volgt. Voor de eerste verge- 
lijking van stelsel A kan men schrijven: 


-4 = es * C3 * Ü е C3 #¿ 2 


ofwel: de reduced revenue van produkt 2 is gelijk aan de bruto- 
marge van produkt 2 verminderd met het aantal еепһедеп van pro- 
dukt 1 dat шеп zou moeten opgeven om één eenheid van produkt 2 
te kunnen maken x de bruto-marge van produkt І en eveneens ver- 
minderd met het aantal eenheden van produkt 3, dat men zou moe- 
ten opgeven om één eenheid van produkt 1 te kunnen maken x de 
bruto-marge van produkt 3. 

Een soortgelijke interpretatie kan men geven aan de andere drie 
vergelijkingen van stelsel А, met dit verschil, dat men bij de 
slack-variabelen spreekt over imputed costs іп plaats van over 
reduced revenue. De derde vergelijking bijvoorbeeld kan men als 
volgt interpreteren: de imputed costs van de eerste restrictie 
is gelijk aan het aantal eenheden van produkt 1 dat men moet 
opofferen om één capaciteitseenheid van restrictie І ongebruikt 
te laten x de bruto-marge van produkt 1 plus het aantal eenhe- 
den van produkt 3 dat men moet opofferen om Eén capaciteitseen- 
heid van restrictie I ongebruikt te laten x de bruto-marge van 
produkt 3. | 


Uit bovenstaande interpretatie Кап men tevens concluderen, dat 
de optimale oplossing bij verandering van de waarden van cj, C2, 
сз en с, de optimale oplossing blijft, zolang geldt dat de мааг- 
den van de uitdrukkingen: 


“ба т бі $. D 9 са * 2 
«січові 2 4766 7009 
дз gy « Cel) жж Z 
0 +c) Г +. cs k teal) 
alle > 0 zijn. 
Vragen wij ons af binnen welke grenzen nu с; kan varieren zonder 


dat de optimale oplossing verandert, dan moet blijkbaar gelden 
voor een nieuwe waarde 62: 


- сі%8ж%0%5х2>0 
ofwel 


63 € 10 


Het blijkt in de praktijk gemakkelijker te zijn om te rekenen 
met de verandering in с». Deze noemen wij Лс (er geldt dus 


с: = сә + Aes). Wij vinden dan: 


- (6 + Ас) + 8 «0 + 5 ж 2 > 0 
ofwel 

Aes < 4 
Deze grens kunnen wij onmiddellijk іп het laatste simplex-ta- 
bleau aflezen. Wij constateren dus: 
Bij toeneming van de bruto-marge van een niet-basisvariabele 
blijft de optimale oplossing ongewijzigd, zolang deze toeneming 
kleiner is, dan de coefficient van die variabele in de £-ri] 
van het laatste tableau. Anders gezegd: bij toeneming van de 
bruto-marge van een niet-basisvariabele blijft de optimale op- 


lossing ongewijzigd, zolang deze toeneming kleiner is dan de 
absolute waarde van de reduced revenue van die variabele. 


Voor basisvariabelen ligt de situatie iets gecompliceerder. Kie- 
zen wij als voorbeeld variabele 1, dan weten wij dat de optimale 
oplossing niet verandert zolang aan de volgende 4 ongelijkheden 
is voldaan: 

ыт Еее ROA 5*2» 0 

“RO (ER A03 є 2-45 w (el), > 0 

Се 08 де) (41) + 5» 2 + Q 

ü * (& * Bc) * 12.754 (2-1) » 0 


Voor deze vier vergelijkingen kan men ook schrijven: 


A *0* ACY > 0 
Рт 2 %* Де; > 0 
a” iE AGT > 0 
з * до: > 0 
Aan deze vier vergelijkingen moet tegelijkertijd zijn voldaan. 


Er moet dus gelden 7 +2Acj > 0 en 3 + Ac > 0 ofwel Ac, > - 31 
en Ас) > - 3. De grootste beperking wordt gevormd door 
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Ac, > - 3, zodat de benedengrens voor Ас) gelijk is aan -3. 
De bovengrens wordt gegeven door Ac; < 2. Wij hebben nu dus ge- 
vonden: 


es A < AC] ec» 


De bovenstaande vier eenvoudige vergelijkingen kan men ook di- 
rect uit het laatste simplex-tableau afleiden. De getallen 4, 
7, 2 en 3 zijn de coéfficiénten in de z-rij van de niet-basis- 
variabelen en de getallen 0, 2, -І en І zijn de coefficienten 
іп de хі-гі| van de niet-basisvariabelen. 
Men kan derhalve de grenzen, waarbinnen Ac; mag varieren als 
volgt bepalen. Vorm de breuken met als teller de coefficienten 
уап de niet-basisvariabelen in de z-rij en als noemer de coeffi- 
ciénten van de niet-basisvariabelen in de xj-rij en voeg aan 
deze breuken een -teken toe. Wij krijgen dan de breuken: 

4 7 2 3 


SH po L 


Kies nu van de aldus gevonden getallen van alle negatieve getal- 
len de grootste en van de positieve getallen de kleinste (niet- 
gedefinieerde breuken, zoals 4/0 in dit voorbeeld, blijven bui- 
ten beschouwing). Wij vinden dan -3 en +2: 


өз 261. < 2 
ofwel (с? = 8 + Aci): 
e Е? € TO 
Bovenstaande gevoeligheidsanalyse berust geheel op de veronder- 


stelling, dat alle overige gegevens (dus ook de andere coeffi- 
ciënten van de doelstellingsfunctie) ongewijzigd blijven. 


7.0. GEVOELIGHEIDSANALYSE VAN HET RECHTER LID 


In hoofdstuk 3, 4 en 5 bespraken wij het probleem van het opti- 
male produktiepakket: 


maximeer 
45x + 30x2 

onder de voorwaarden 
2хі ж ко 5 129 І 
хі + 2х; < 140 II 


IA 


44.045 X4 80 ^ III 
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Іп paragraaf 3.1. hebben wij reeds laten zien dat voor een ver- 
ruiming van de eerste restrictie met één eenheid (d.w.z. voor 

bj = 121 in plaats van bj = 120) de waarde van de doelstellings- 
functie toeneemt met de schaduwprijs van restrictie | (in dit 
voorbeeld behoort bij restrictie 1 een schaduwprijs uj = 15). 

De vraag is nu: binnen welke grenzen kan bj varieren zonder dat 
de waarde van de schaduwprijs van restrictie 1 verandert? 

Om deze vraag te beantwoorden beginnen wij met de volgende ор- 
merking. Een verandering van bj kan de optimale oplossing op 
twee manieren beïnvloeden. In de eerste plaats kan een veran- 
dering van bj leiden tot andere waarden van de basisvariabelen, 
terwijl de verzameling basisvariabelen ongewijzigd blijft. Dit 
betekent, dat in het laatste simplex-tableau alleen in de laats- 
te kolom (die van het rechterlid, waarin de waarden van de 
basisvariabelen staan) wijzigingen optreden, maar dat alle an- 
dere kolommen ongewijzigd blijven. Ook de schaduwprijzen blijven 
dan ongewijzigd. 

In de tweede plaats kan de verandering in b, zo groot zijn, dat 
andere variabelen in de basis worden opgenomen. Dit leidt tot 
een andere verzameling basisvariabelen en tevens tot andere 
schaduwpri jzen. 

Wij concluderen dus, dat de schaduwprijs van restrictie | (nl. 
uj; = 15) ongewijzigd blijft zolang de verandering van bj zo 
klein is, dat de huidige verzameling basisvariabelen D Kaj 
Xu) ongewijzigd blijft (zie tabel 7.2.). 


- tabel 7.2. Laatste simplex-tableau voor het probleem van 
het optimale produktiepakket - 


Stel nu dat de waarde van bj wordt gewijzigd in bl = 121. Wij 
zouden dan het probleem opnieuw kunnen oplossen en zouden dan 
vinden het simplex-tableau van tabel 7.3.: 


- tabel 7.3. Laatste simplex-tableau als ы = pe] ~” 


In feite hebben wij het tableau van tabel 7.3. niet gevonden 
door het probleem opnieuw met behulp van de simplex-methode op 
te lossen, maar door als volgt te redeneren. In de kolom van 

de bij restrictie 1 behorende spelingsvariabele staat hoeveel 
eenheden wij van de basisvariabelen moeten opofferen om één 
eenheid van die spelingsvariabele in de oplossing te kunnen іп- 
troduceren. Verruiming van restrictie 1 met één eenheid (d.w.z. 
wijziging van bj = 120 in bi = 121) is precies het tegenoverge- 
stelde. Wij kunnen dus door de verruiming met één eenheid van 
restrictie 1 één eenheid xj en één eenheid x, extra in de oplos- 
sing brengen, terwijl wij één eenheid x? moeten opgeven (Opm. : 
de schaduwprijs van restrictie | is gelijk aan 

(La 45 ASTI k 0 = 1 * 305» 15, d.w.2. gelijk aan dé extra eenbe- 
den van de basisvariabelen die wij in de oplossing kunnen intro- 
duceren c.q. moeten opgeven vermenigvuldigd met de dekkingsbij- 
dragen per eenheid). Wij vinden dus: 


3000 15 3015 
Nieuwe oplossing - K + E E 41 
20 1 21 
40 "i 39 
Veranderen wij de waarde van bj in Ы = 120 + ЛЬ}, dan is de 
nieuwe oplossing: 
3000 15 3000 + 154b) 
Nieuwe oplossing = 29 + Ab} b L. 40% Д; 
20 1 20 з» Ab) 
40 e, 40-9 . Ab 
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Zolang de waarden van alle variabelen in deze nieuwe oplossing 
niet negatief zijn, is deze nieuwe oplossing feasible en dus 
optimaal, Er moet dus gelden: 

40 + Ab) > 0 

20 + Ab; 

40 Ж Abi 


IV 


0 


IV 


Aan deze drie ongelijkheden is voldaan als geldt 
"20. ЛЬ; < 40 


Wij concluderen dus dat de schaduwprijs u] = 15 ongevoelig is 
voor wijzigingen in bj, zolang geldt 


100 < bi < 160 


Om deze "range of feasibility" voor Ab; te vinden, kan men dus 
de volgende procedure toepassen: Vorm de breuken met als teller 
de waarden van de basisvariabelen in het rechterlid en als noe- 
mer de coefficiënten in de kolom van de bij restrictie 1 beho- 
rende spelingsvariabele en voeg een -teken toe. Wij krijgen dan: 


и g 
ke? P < 


Kies nu van de negatieve getallen de grootste еп van de рові- 
tieve getallen de kleinste, Dit zijn de grenzen voor Ab): 


-20 < hb, < 40 


De boven beschreven procedure geldt alleen voor een verandering 
van < -restricties. Voor een verandering van een > -restrictie 
geldt een soortgelijke procedure. In plaats van de getallen in 
de kolom van de spelingsvariabele neemt men dan de getallen in 
de kolom van de surplus-variabele, terwijl men bovendien aan 

de zo gevormde breuken geen -teken behoeft toe te voegen. Het 
bewijs van de juistheid van deze procedure laten wij aan de 
lezer over. 
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7.3. VRAAGSTUKKEN 
7-1. Wij herhalen de probleemstelling van vraagstuk 4-1: 
Produkt А В C Capaciteit 


Dekking per 7 50;- У З» 7 60,- 


eenheid 

Capaciteitsbeslag 3 1 2 2000 iat 

in manuren 

Capaciteitsbeslag 1 2 2 1400 (machine-uren 
machines type І type 1) 
Capaciteitsbeslag 2 2 1 1500 (machine-uren 
machines type ІІ type 11) 


Het laatste simplex-tableau luidt: 


Hierin zijn x4, хэ en хз echte variabelen (produkten A, B en C) 
епоху, X5 еп хе spelingsvariabelen behorende bij de restricties 
manuren, machine-uren type I en machine-uren type II. 


a. Zou men het verlies van deze onderneming tot nul kunnen reduce- 
ren, als men 25 man zou kunnen ontslaan? 


b. Stel dat de grondstofkosten van produkt A toenemen van f 60,- 
tot f 70,- per eenheid produkt. Heeft deze wijziging invloed 
op de optimale samenstelling van het produktiepakket? 

Heeft deze wijziging invloed op het verlies van de onderne- 
ming per dag? Zo ja, hoe groot is dit verlies nu? 


c. Stel dat de onderneming een mogelijkheid ziet om de verkoop- 
prijs van produkt B te verhogen. Hoe hoog zou deze verkoop- 
prijs moeten zijn teneinde opneming van produkt B in het | 
produktiepakket voor de onderneming aantrekkelijk te maken? 
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7-2. Wij herhalen de probleemstelling van vraagstuk 4-2. b: 


Produkt A B eure Epara 
capaciteit 

Dekking per 

eenheid 7.200,--.7 300,- 7 250;- 

Capaciteitsbeslag 

afdeling I 2 4 4 1800 

Capaciteitsbeslag 

afdeling II 4 5 2 3000 


Het laatste simplex-tableau luidt: 


1/3 gS 
*1/6 1/3 


Hierin zijn X], X en хз echte variabelen (produkten A, B en C) 
en x, en xs spelingsvariabelen behorende bij afdeling I en II. 


a. Stel, dat de dekkingsbijdrage van produkt B toeneemt tot 
f 330,- per eenheid, wordt dit produkt dan in het optimale 
produktiepakket opgenomen? | 


b. Stel dat de dekkingsbijdrage van het nieuwe produkt C daalt 
tot f 225,- per eenheid, verandert dan de samenstelling van 
het optimale produktiepakket? 

Heeft deze wijziging invloed op de behaalde dekking? 
Zo ja, hoe groot is dan nu de behaalde dekking? 


c. Kunt u zonder het probleem opnieuw met behulp van de simplex- 
methode op te lossen berekenen met welk bedrag de behaalde 
dekking toeneemt als de beschikbare capaciteit van afdeling 
II toeneemt van 3000 tot 3500 uur? Zelfde vraag voor een 
toeneming van de capaciteit van afdeling II van 3000 tot 
4000 uur. 
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7-3. Wij gaan uit van het іп vraagstuk 4-6 gegeven probleem. 
Het laatste simplex-tableau behorende bij dit probleem luidt: 


a. Stel dat de verkoopprijs van halffabrikaat 1 daalt van f 17,-- 
tot f 14,50, verandert dan de samenstelling van het optimale 
produktiepakket? Heeft deze daling van de verkoopprijs van 
halffabrikaat 1 invloed op de in totaal behaalde dekking? Zo 
ja, hoe groot is dan nu de in totaal behaalde dekking? 


b. Als de verkoopprijs van halffabrikaat 2 toeneemt van f 32,-- 
tot f 35,--, wordt halffabrikaat 2 dan opgenomen in het opti- 
male produktiepakket? 


c. Op de engineering-afdeling krijgt iemand anders nog een idee. 

` Door een kleine wijziging in het ontwerp van de nieuwe instal- 
latie lijkt het mogelijk om per dag 200 kg van stof A aan het 
afvalwater te onttrekken in plaats van 100 kg. De vaste kosten 
van de installatie nemen door deze wijziging toe van f 1200,-- 
per dag tot 7:1230;%- per dag. 
Kunt u uitgaande van de gegevens in bovenstaand simplex-tableau 
berekenen of het in gebruik nemen van de gewijzigde installatie 
nu - bedrijfseconomisch gezien - wel aantrekkelijk is? 


8. ENKELE MET LINEAIRE 
PROGRAMMERING VERBAND 
HOUDENDE ONDERWERPEN 


8.1. MEERVOUDIGE DOELSTELLINGEN 


Het komt dikwijls voor, dat men in een praktijksituatie te ma- 
ken heeft met meervoudige doelstellingen. Bij wijze van voor- 
beeld beschouwen wij het volgende probleem. Een onderneming 
wenst de samenstelling van het produktenassortiment te kiezen 

op zodanige wijze, dat enerzijds de dekkingsbijdrage zo groot 
mogelijk is, terwijl men anderzijds wenst te streven naar een 

zo groot mogelijke werkgelegenheid. De onderneming kan twee pro- 
dukten maken, waarvan de gegevens zijn vermeld in tabel ы” 


Produkt 1 Produkt 2 
Verkooppijs 49 64 
Grondstofkosten 15 20 
Dekking I 34 44 
Loonkosten 20 40 
Dekking ІІ 14 4 
«сарв 8,1, ~ 


Er zijn twee machines; het capaciteitsbeslag per machine per 
eenheid produkt en de maximaal beschikbare capaciteiten zijn 
gegeven in tabel 8.2.: 


Produkt 1 Produkt 2 Maximale capaciteit 
Machine 1 4 2 280 
Масһіпе 2 2 2 200 
Arbeid 2 4 


“-Хазі6,2:.- 
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Wij kunnen dit probleem formuleren als een lineair-programme- 
ringsprobleem als wij één van de twee doelstellingen opnemen 
als een restrictie en de andere doelstelling beschouwen als de 
doelstellingsfunctie. Wij kiezen voorlopig het streven naar een 
zo groot mogelijke dekkingsbijdrage als doelstellingsfunctie en 
formuleren het streven naar een zo groot mogelijke werkgelegen- 
heid als restrictie. 


Het aantal produktief bestede arbeidsuren is 2x; + 4x; als мі) 
xj eenheden van produkt 1 еп x eenheden van produkt 2 maken. 
Daarnaast kunnen wij besluiten een bepaald aantal arbeidsuren 
niet te gebruiken. Om deze mogelijkheid in het model tot uik” 
drukking te brengen voeren wij in een fictief produkt 3, waarvan 
de gegevens zijn vermeld in tabel 8.3. 


Kostprijsgegevens: Capaciteitsbeslag: 
Verkoopprijs 0 Machine 1 0 
Grondstofkosten 0 Machine 2 0 
Dekking І 0 Arbeid 1 
Loonkosten 10 
бакқа id -10 


- tabel 8.3. Gegevens fictief produkt 3 - 


Het model kan nu als volgt worden geformuleerd: 
maximeer 


z з l4xj + Ах; - 10x3 


onder de voorwaarden 


аху + 2х2 < 280 machine | 
2x] + 2х; < 200 machine 2 
2x| + йхо + хз 2 У werkgelegenheid 


Wij zijn nu met name geïnteresseerd іп de relatie tussen z en 
у, met andere woorden wij zoeken de vorm van de functie 

z = f(y). 

Wij lossen daartoe het bovenstaande LP-probleem enkele malen 
op voor verschillende waarden van y. Voor y = 100 krijgen wij 
het volgende probleem: 
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2 - 14х) - 4х) + 10x3 ою 
áx) + 2х, + XL = 280 
які +. 28o * xe - 200 
2х] + 4X2 + хз 7 Xg * 100 


De simplex-tableaus zijn gegeven in tabel 8.4. In tabel 8.4. 
hebben wij geen kunstmatige variabele behoeven te introduceren 
omdat bij de laatste rij zeer gemakkelijk хз als basisvariabele 
kan worden geintroduceerd. Wij hebben tableau 1 laten ontstaan 
door de met -10 vermenigvuldigde laatste rij van tableau 0 op 
te tellen bij de eerste rij van tableau 0. 


De optimale oplossing is z з 980, х] = 70, x = 0, хз = 0, 
X5 = 60, xg - 40. De surplus variabele behorende bij 

2x] + 4x2 + x3 2 100 is dus gelijk aan 40. Dit houdt in, dat 
voor waarden van y < 140 geldt: z = 980. 


| 2 
4 X6 
X5 
XI 
Wij lossen nu het probleem op soortgelijke wijze op voor у = 200. 
Het laatste simplex-tableau is gegeven in tabel 8.5. 


De optimale oplossing is nu: z = 920. Wij zijn nu speciaal gein- 
teresseerd in de schaduwprijs van de laatste restrictie. Deze 
is gelijk aan 1. Voor у = 199 geldt nu z = 921 en voor y = 201 
geldt 2 - 919. De grenzen waarbinnen deze schaduwprijs geldig 


blijft vinden wij door op te schrijven: 


60 40 20 
1/9011. 1/3 * "A3 
360, 120 

-60 


De grenzen voor Ay zijn dus 
-60 < Ay « 120 

ofwel 
140« y « 320 


Binnen deze grenzen is dus de schaduwprijs van de laatste ге- 
strictie gelijk aan l. 
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Wij lossen nu het probleem nogmaals op voor y = 360. Het laatste 
simplex-tableau is gegeven in tabel 8.6. 


TEDA 8.0; 


Voor у = 360 geldt dus 2 = 600. De schaduwprijs van de werkge- 
legenheidsrestrictie is 5 en de grenzen waarbinnen deze schaduw- 
prijs constant blijft zijn: 
40 80 CG 

тах(---, St < Ау < вітро 
ofwel 

-40 < Ay < 40 
ofwel 

320 < y < 400 


Tenslotte lossen wij het probleem nog eenmaal op voor y = 500. 
Het laatste tableau is dan gegeven in tabel 8.7, 


^ Tabel 8.7, ~ 


Voor у - 500 geldt dus 2 - -600. De schaduwprijs van de werkge- 
legenheidsrestrictie is 10, Deze schaduwprijs blijft constant 

100 
"I 


zolang < Ay ofwel zolang y > 400 is. 
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Wij kennen пи de waarde уап 2 voor iedere waarde van y. In fi- 
guur 8.1. hebben wij z = f(y) in beeld gebracht. 


Wij voeren nu in het begrip effictênte oplossing. Wij gaan uit 
van een probleem met meer dan één doelstelling, waarbij deze 
doelstellingen alle gemaximeerd moeten worden. Wij noemen 

dan een oplossing efficient, indien er geen andere oplossing 
bestaat die ten aanzien van tenminste één doelstelling een һо- 
gere waarde heeft terwijl de waarden van de andere doelstellin- 
gen niet lager zijn. In figuur 8.l. is P, geen efficiente oplos- 
sing, immers oplossing P2 is qua werkgelegenheidsdoelstelling 
beter dan Pi en qua winstdoelstelling niet slechter. Men zegt: 
oplossing Ру wordt gedomineerd door P5. Evenzo wordt Ру gedomi- 
neerd door Рз. P wordt ook gedomineerd door Ps, А, Pg, В еп 
Pi, Ре ig ook geen efficiënte oplossing, want Ps wordt gedomi- 
neerd door A. Wel efficient zijn alle oplossingen die overeen- 
komen met punten op de gebroken halve rechte ABCD ....... 
Het in figuur 8.1. gearceerde gebied kan men naar analogie van 
de bij lineaire programmering gebruikelijke terminologie de 
verzameling van de feasible solutions noemen. Deze verzameling 
is altijd convex, d.w.z. dat als twee punten P en Q beide tot 
de verzameling behoren alle punten op het lijnstuk PQ ook tot 
de verzameling behoren. 


Uit de verzameling van efficiente oplossingen (deze verzameling 
bestaat dus uit de punten op de gebroken halve rechte 
АВСР......) zal men nu één oplossing moeten kiezen. De tech- 
niek van de lineaire programmering biedt bij deze keuze verder 
geen hulp. 
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Ps (140, 980) 
ў 
d Pe 


320, 800 


od y 


D 


- figuur 8.1. De verzameling van efficiënte oplossingen 
voor het winst-werkgelegenheidsprobleem - 
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8.2. GEHEELTALLIGE PROGRAMMERING 


Tot dusverre hebben wij steeds aangenomen dat de beslissings- 
variabelen x; continue variabelen zijn, die alle waarden > 0 
kunnen aannemen. In de praktijk doet zich echter meermalen de 
situatie voor dat één of meer van de beslissingsvariabelen dis- 
crete variabelen zijn. Een van de mogelijkheden om dit probleem 
te omzeilen is de discrete variabelen te behandelen als conti- 
nue variabelen. De gevonden oplossing wordt dan vervolgens voor 
de discrete variabelen afgerond. Wij beschouwen het volgende 
voorbeeld: 

maximeer 


25X] + 1 5X» 


onder de voorwaarden 


5X| + 4X» < 20 I 
8x, + 4x2 < 24 St 
Х1» Х2 ED 


waar хі, хә gehele niet-negatieve getallen zijn. Indien wij dit 
probleem oplossen door eerst xj en x? als continue variabelen 


te beschouwen vinden wij (zie figuur 8.2.): xj = 1l РУ 3 
Wij willen nu deze oplossing afronden еп vinden dan x; = |, 


хо = 3, z = 70 (bij dit afronden moeten wij er natuurlijk voor 
zorgen dat de oplossing binnen het gebied van de feasible 
solutions blijft; in dit voorbeeld is gemakkelijk in te zien, 
dat wij door naar beneden af te ronden inderdaad een feasible 
solution verkrijgen, bij grotere problemen is dit niet altijd 
even gemakkelijk in te zien). 


In figuur 8.2. hebben wij alle toegelaten geheeltallige oplos- 
singen aangegeven. Tevens zien wij door de figuur nauwkeurig 
te bekijken, dat (2,2) een betere oplossing is dan (1,3). Voor 
xı ® 2, хо 2 geldt 2 е 90. 

Voor problemen als het bovenstaande zijn aparte oplossingsme- 
thoden ontwikkeld. Men duidt de techniek van het oplossen van 
dergelijke problemen aan als geheeltallige рговгамаегіпв (іп- 
teger programming). Zijn sommige beslissingsvariabelen continu 
en andere beslissingsvariabelen geheeltallig dan spreekt men 
van mixed integer programming. 
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^ figudr 8.2. ~ 


8.3. VRAAGSTUKKEN 


8-1. De boer van vraagstuk 4-5 wordt in feite geconfronteerd 
met een keuzeprobleem met twee doelstellingsfuncties. Hij wenst 
enerzijds een zo hoog mogelijke bruto-winst en anderzijds 
zoveel mogelijk vrije tijd. Gedurende de maanden april tot en 
met september kan de boer zelfs als hij dat zou willen niet 
meer dan 1500 uur werken; werkt hij gedurende die periode 
1500 - у uur, dan heeft hij y uur vrije tijd. Zij z de maximale 
bruto-winst, die de boer kan behalen, indien hij tenminste 

y uur vrije tijd wenst. 


а. Maak een grafiek van de functie z = £05. 
b. Wat is de verzameling van efficiënte oplossingen? 


8-2. Een onderneming beschikt over een leegstaande bedrijfshal 
met een vloeroppervalkte van 1250 m?. Men overweegt deze hal te 
gaan benutten, door er nieuw aan te schaffen machines іп te 
installeren. Daarbij kan men kiezen uit twee typen machines, 


98 


die wij aanduiden met machines van type I en machines van type II. 
Voor één machine van type I is een vloeroppervlak van 200 m? 
nodig, voor één machine van type II een vloeroppervlak van 300 m°. 
De onderneming kan ten hoogste f 1,5 mln besteden voor de aanschaf 
van de nieuwe machines. De aanschafprijzen van de machines zijn: 

f 400.000,-- voor een machine van type I еп f 300.000,-- voor 

een machine van type II. 

Door het aanschaffen van machines van type I neemt de jaarlijkse 
bruto-winst toe met f 10.000,-- per aangeschafte machine. Voor 
machines van type II is dit 7 12.000,--. 


a. Formuleer dit probleem als een LP-probleem. 

b. Los het probleem op met behulp van de grafische oplossings- 
methode alsof de beslissingsvariabelen ook gebroken (niet- 
geheeltallige) waarden konden aannemen. Rond de gevonden uit- 
komst naar beneden af op gehele getallen. 

c. Is de zo gevonden oplossing de optimale oplossing voor dit 
probleem? Zo niet, wat is dan wel de optimale oplossing? 


8-3. Gegeven is een onderneming, die in principe drie produkten 
(A, B en C) zou kunnen maken. Van deze drie produkten zijn de 
volgende gegevens bekend: 


A B C 
Verkoopprijs 110 100 105 
Grondstofkosten 60 70 45 
Dekking I 50 30 60 
Loonkosten 60 20 40 
Dekking II -10 10 20 


Verder zijn de volgende gegevens bekend: 

- er zijn 200 machines van type I; na aftrek van tijd nodig voor 
onderhoud zijn per dag beschikbaar 1400 uren voor machines van 
type I 

- er zijn 200 machines van type II; na aftrek van tijd nodig voor 
onderhoud zijn per dag beschikbaar 1500 uren voor machines van 
type II 

- het capaciteitsbeslag van de produkten is als volgt: 


А В С 
manuren 3 1 2 
machines І 1 2 2 
machines II 2 2 1 


- de directe loonkosten bedragen f 20,-- per manuur 

- de kosten van salarissen van niet-direct produktief personeel, 
de afschrijvingen en de overige vaste kosten bedragen in totaal 
f 10.000,-- per dag. 
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De leiding van deze onderneming wenst twee doelstellingen па te 
streven: enerzijds een zo hoog mogelijke winst, anderzijds een 
zo groot mogelijke werkgelegenheid. 


a. Noem de behaalde winst per dag z en het aantal betaalde man- 
uren per dag y. Bepaal dan z als functie van y en teken de 
grafiek. | 

b. Wat is de verzameling van efficiënte oplossingen? 


T PRODUKTIEPLANNING 
MET BEHULP VAN 


LINEAIRE PROGRAMMERING 


9.1. PROBLEEMSTELLING 


Wij zullen in dit hoofdstuk laten zien hoe men vraagstukken met 
betrekking tot produktieplanning kan oplossen met behulp van 

lineaire programmering. Wij beschouwen daartoe een sterk ver- 

eenvoudigd voorbeeld. | 

Een onderneming produceert twee produkten A en B. De vraag naar 

die produkten in de komende drie perioden is bekend (zie tabel 


9.1. E 
periode 
1 2 3 
Produkt A 1000 2500 6000 
Produkt B 1500 1000 5000 


- Tabel 9.1. Vraag per periode - 


Bij deze onderneming kunnen drie capaciteitssoorten worden оп- 
derscheiden, п.1. machines, arbeid en opslagcapaciteit. De be- 
schikbare capaciteiten zijn gegeven in tabel 9.2.: 


periode 
1 2 3 
Масһіпев (чгеп) 560 560 560 
Arbeid (uren) 350 300 350 
Opslagcapaciteit 


(m? magazijnruimte) 8000 8000 8000 


- Tabel 9.2. Beschikbare capaciteiten per periode - 
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Het capaciteitsbeslag рег eenheid produkt voor de produkten А 
еп B is gegeven in tabel 9.3.: 


Produkt A Produkt B 


Machines (uren) 0,12 0,1 
Arbeid (uren) 0,05 0,08 
Opstagcapaciteit (m?) 1.3 2.0 


- Tabel 9.3. Capaciteitsbeslag per eenheid produkt - 


In tabel 9.2. zijn het aantal beschikbare machine-uren en ar- 
beidsuren gegeven. Dit aantal uren betreft het aantal uren in 
de normale werktijd, Daarnaast is overwerk mogelijk en wel tot 
maximaal 504 van de in een bepaalde periode binnen de normale 
werktijd beschikbare uren. 

De werkelijke voorraden aan het begin van periode | en de mini- 
maal gewenste voorraden aan het einde van periode 3 zijn gege- 
ven in tabel 9.4. 


Voorraad aan het begin Gewenste voorraad aan het 


van periode 1 eind van periode 3 
А 900 800 
В 400 600 
- Tabel 9.4 > 


De kosten van overwerk bedragen f 15,- per overgewerkt uur. De 
voorraadkosten bedragen 1,47 van de kostprijs van het opgesla- 
gen produkt per periode. De kostprijs van produkt A is f 20,- 
van produkt B 7 15,- per eenheid, Bij welk produktieprogramma 
is de som van voorraadkosten en overwerkkosten minimaal? 


9.2. FORMULERING VAN HET MODEL 


Allereerst geven wij aan welke variabelen wij als beslissings- 
variabelen in het model wensen op te nemen. Noemen wij A1, А», 
Аз, Ві, B2, B3 de te produceren hoeveelheden gedurende de door 
de index aangegeven periode еп SA], 5А;, SA3, SB], 5В2, SB; de 
voorraden van het betreffende produkt aan het eind van de door 
de index aangegeven periode, dan geldt: 


ЗА; = 900 + A, - 1000 

SA, = 900 + Aj - 1000 + A 

SA3 = 900 + Aj - 1000 + А 

SB) = 400 + Bj - 1500 

SB, = 400 + Bj - 1500 + В 

SB3 = 400 + Bj - 1500 + B 
ofwel 

SA) = Aj - 100 

SA, = Aj + А; - 2600 

SA3 = Aj + А; + Аз - 8600 

58; = Bj - 1100 

SB, = Bj + В2- 2100 

SB3 = Bj + В, + Вз - 7100 


Wij kunnen nu bij het constr 
belen А; еп В; ofwel de vari 


2 E 2500 
фо 2500 + Аз - 6000 
2 - 1000 
> - 1000 + Вз - 5000 


ueren van het model ofwel de varia- 
abelen SA; en SB; (i = 1, ..., 3) 


als de beslissingsvariabelen gaan beschouwen. Kiest men Aj en 


В; (і 1, ..., 3) als de be 
het model construeren als vo 
Allereerst specificeren wij 


slissingsvariabelen, dan kan men 
lgt: 
de restricties. Wij beginnen met de 


9 capaciteitsrestricties (één restrictie per capaciteitssoort 


en per periode). 


De capaciteitsrestricties in 


luiden als volgt: 


MAC LT OFT * ЧІВ; S 
MAC 2t :0,12А2 * 0,1В9:5 
MAC 3: 0,1243 + 0,153 S 
LAB 1: 70,05А1 * 9,088; 5 
ҺАВ 2: 0,05А2 + 0,08В2 < 
LAB 3: 0,05Аҙ + 0,08B3 < 


Bij bovenstaande restricties 
nomen het aantal beschikbare 
De restricties inzake de ops 


zake de machine-uren en arbeidsuren 


840 

840 
840 

529 

450 

525 

hebben wij in het rechterlid opge- 


uren inclusief 507 overwerk. 
lagcapaciteit luiden als volgt: 
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STO}: 1,5(Aj -100) + 2(Ві - 1100) < 8000 

STO»: 1,5(А) +A, - 2600) + 2(B + BJ -2100) < 8000 

5Т03: 1,5(Aj * A5 +A3 - 8600) + 2(Ву +В, + Вз - 7100) < 8000 
In iedere periode moet tenminste zoveel worden geproduceerd dat 
aan de vraag kan worden voldaan, Wij kunnen dit tot uitdrukking 
brengen door te eisen SA; > 0 en SDi- S D voor 3 wi, sii, 3, In 
de laatste periode moet zelfs gelden 5Аҙ > 800 en SB3 > 600. De 
restricties SA3 > 0 en SB3 > 0 worden daardoor overbodig. Wij 
krijgen derhalve zes vraag-restricties: 

DEA 1: А) ~ 1002 0 

DEA 2: А) + А; - 2600 > 0 
DEA 3: Aj + А; + Аз - 8600 > 800 

DEB 1: Bj - 1100 > 0 

DEB 2: В) + B; - 2100 > 0 

DEB 3: Bj + В, + Вз - 7100 > 600 
Wij gaan nu over tot het specificeren van de te minimeren doel- 
stellingsfunctie. Deze luidt: 

0,28(A1-100) + 0,28(Aj+A5-2600) + 0,28(A1*A5*A4-8600) + 

+ 0,21(В1-1100) + 0,21(B1*B5-2100) + 0,21(Bj+B5+B3-7100) Ф 

- 15 01 * 15 02 Ф 15 Оз 
waar 0), 0» en 03 voorstellen het aantal uren overwerk іп de 


perioden 1, 2 еп 3. Het aantal uren overwerk in de eerste perio- 
de wordt gegeven door: 


01 - 0,05A, + 0,08В) - 350 
als de uitdrukking іп het rechterlid > 0 is en 0) = 0 als 


0,05A, + 0,08В) - 350 < 0 

is, 

Wij nemen nu in het model op de gelijkheidsrestrictie 
01 = 0) = 0,05A1 + 0,088; * 330 


waarbij wij dus tevens 0) еп 0; als beslissingsvariabelen intro- 
duceren. Wij stellen vanzelfsprekend voor iedere beslissingsva- 

riabele, dus ook voor 0) en 0) de eis van non-negativiteit, dus 

01 > 0 en u, 2 0, 
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In de optimale oplossing zal nu altijd gelden ofwel 0; > 0 еп 
01 = 0 ofwel 0; = 0 en Оу > 0 doch nimmer O, > 0 en U > 0. Im- 
mers indien 0; > 0 еп Ui > 0, dan zou men zeer gemakkelijk een 
betere oplossing kunnen vinden door 0) еп 0) beide zodanig te 
verlagen dat juist een van beide juist gelijk aan nul wordt. 


Wij moeten derhalve de volgende drie restricties nog aan het 
model toevoegen: 
OVE 1: 0) - Uy = 0,05Aj ж 0,088; - 350 
OVE 2: 02 - 02 Р 0,05А2 - 0,08В: - 350 
OVE 3: 04 ^ Us. 0,0543 + 0,08B3 - 1350 
Het volledige model luidt nu als volgt 10), 
minimeer 
0,84A1 + 0,56А2 + 0,28Аҙ + 0,63Bj + 0,42В2 + 0,2183 + 
+ 15011115 Dg + 15 98 


onder de voorwaarden 


MAC 1: 0,12А) + 0,1B < 840 
MAC 2: 0,12А; + 0,1В2 < 840 
MAC 3: 0,12A4 + 0,1B3 < 840 


LAB La 0,О5А; Ж 0,088; 5 525 
0,05А; + 0,088, < 450 
0,05A4 + 0,08B3 < 525 
STO 1: LB" 4-28 $o 10350 
2 1,5Aj + 1,5А2 + 2B, + 2В2 < 16100 
STO 3: ‘ISA “3,540 Є I, SAGE 280% 48р + 28; < 35100 
OVE 1: 0,05Aj + 0,088; - 0] + Up = 350 

2: 0,05А; + 0,08В; - 02 + 0; = 300 
OVE 3: 0,05Аҙ + 0,08B3 ~ Оз + Оз = 350 
DEA 1: Aj 2 100 

2: А, - Ав 2 ВИН 
DEA 3: А) + А2 + Аз > 9400 
DEB 1: Bj > 1100 
DER 2: -Bj + 82% %100 
DEB 3: Pi * Во 7/53 2 7700 


t^ 
> 
оз 
(о N 
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9.3. OPLOSSEN VAN НЕТ MODEL 


Het model kan nu worden opgelost met behulp van de simplex- 
methode. Hiervoor kan men gebruik maken van IMOPTI of van OPTI- 
MIZER. Voor de wijze van invoer van de gegevens verwijzen wij 
naar de beschrijvingen van deze beide Programma's. IMOPTI draait 
alleen op grote computers; een beschrijving van IMOPTI is ver- 
krijgbaar bij het Instituut voor Ontwikkeling van het Wiskunde 
Onderwijs te Utrecht. OPTIMIZER draait op verschillende microcom- 
puters (voor de IBM-Personal Computer is de minimale configuratie 
48K met disk drive en printer). De manual en het programma zijn 
verkrijgbaar bij SOFTKEY B.V. te Deventer. 


De inputgegevens en de output van IMOPTI drukken wij hierna af ор 
de pagina's 106 еп 107. De inputgegevens en de output van OPTI- 
MIZER zijn afgedrukt op de pagina's 108, 109 en 110. 


De output van IMOPTI bestaat alleen uit het laatste simplex- 
tableau. De output van OPTIMIZER bestaat uit drie onderdelen. 
Eerst geeft OPTIMIZER van het laatste simplex-tableau het rechter- 
lid (de waarden van de doelstellingsfunctie, van de beslissings- 
variabelen en van de slack en surplus variabelen) en de z-rij (de 
schaduwprijzen en de reduced costs). Vervolgens geeft OPTIMIZER 
de resultaten van de gevoeligheidsanalyse van de coëfficiënten 
van de doelstellingsfunctie en tenslotte de resultaten van de 
gevoeligheidsanalyse van het rechterlid. 


Wij bespreken nu eerst de optimale oplossing. Deze luidt: 


Aj = 100 
A> з 5522,22 
A3 = 3777,78 
Bj = 3534,72 
B2 = 298,61 
B3 = 3866,67 


Van de overige beslissingsvariabelen komen alleen 0) en Оз in 
de basis voor. Uj - 62,22 hetgeen inhoudt dat 62,22 van de in 
periode І іп de normale werktijd beschikbare uren niet worden 
gebruikt. 03 - 148,22 hetgeen inhoudt dat in periode 3 geduren- 
de 148,22 uur overwerk zal moeten plaatsvinden. 0, еп U, komen 
niet in de basis voor. In periode 2 kunnen alle іп de normale 
werktijd beschikbare uren nuttig worden gebruikt, terwijl er in 
periode 2 geen overwerk behoeft plaats te vinden. Deze zelfde 
informatie kunnen wij ook verkrijgen door de waarden van de 


L INEAIR-FROGRAMMERINGSFROGRAM IMOFTI 
WILT U DE GEGEVENS VIA ТАРЕ INLEZEN?TYPE JA OF NEE ? 
CM LWA+1 = 25077В, LOADER USED 41400Вгее 


GEEF LE VOLGENDE SPECIFICATIE VAN HET РКОВІ.ЕЕМ: 
AANTAL BESLISSINGSVARIABELEN ? 12 


TOTAAL AANTAL RESTRICTIES 718 
AANTAL =“ RESTRICTIES 7 9 
AANTAL = RESTRICTIES ? 3 
AANTAL >= RESTRICTIES ? 6 


MOET ПЕ DOELSTELLINGSFUNKTIE GEMAXIMEERD OF GEMINIMEERD WORDEN 
ТҮРЕ МАХ OF MIN ? min 


GEEF DE 12 NAMEN VAN DE BESLISSINGSUARIABELEN?T 
81 "82 э,аз rbi »b2 »b3 01 702 „o3 rul ru2. sud г 


GEEF DE BIJBEHORENDE COEFFICIENTEN VAN DE DOELSTELLINGSFUNKTIE IN 
DEZELFDE VOLGORDE? 
.84».56».28».63».42».21,15,15»,15,0,0»0 
GEEF ACHTEREENVOLGENS VOOR ELK VAN DE 18 RESTRICTIES EEN REGEL MET DE 
NAAM VAN ПЕ RESTRICTIE GEVOLGD DOOR DE TECHNISCHE COEFFICIENTEN EN HET 
RECHTERL ID? 
DE 9 =< RESTRICTIES 
1?maci/.12»0»0».1»0»0»0»0»0»0»50»0»840 
2?mac2/0» .12»0»0» .1»0»0»0»0»07050»840 
3?mac3/0»0».12»0»0» .1»0»0»0»0»0»507840 
4?1abi/.05»0»0».08»0»0»0»0»0»0»0»0» 525 
5?1ab2/0».05»0»0».08»0»0»0»0»070»05 450 
6?1ab3/0»0».05»0»0».,08»0»0»0»0»0»0»523 
7?65101/1.5,0,0, 2.0, 0» 0, 0, 0,0,0, 010350 
8?5102/1.5,1.5, 0, 2,2,0. 0, 0,0,0,0,0,16100 
9?sto3/1.5»1.5»1*5»2»2»2»0»0»0»0»0»505 35100 


DE 3 - RESTRICTIES 
1?0vel/.,05»0»0».08»0»0»-1»0»0»1»0»0» 350 


2?0ve2/0».05»0»0».08»0»0»-1»0»0»1»0» 300 
3?0ve3/0»0».05»0»0».08»050»-1»0»0»1» 350 


DE 6 »- RESTRICTIES 
1?deal1/1»0»0»0»0»0»0»0»050»0»0*100 


2?dea2/1»1»0»0»0»0»0»0»0»0»0»0»2600 
3?dea3/1»1»1»0»0»0»0»0»0»0»0»0»7400 
A?debi1/0s0»0»1»0»0»0»0»0»070»0»1100 
5?deb2/0s0»0»1»1»0»0»0»5050»0»0»2100 


&?deb3/0s0»s0»1»1»1»0»0»0»0»0s0»7700 


ECHTE VARIABELEN ^1 А2 АЗ кі B2 B3 01 02 03 
U1 92 уз 

ПЕА1 DEA? ПЕАЗ DERI DEB2 DERZ 

МАС1 масо MAC3 LAB1 LAE? LAB3 65701 STO? STO3 


SURFLUS VARIABELEN 
SLACK VARIABELEN 


ARTIFICIAL VARIABELEN OVE1 OVE? OVE3 
ART 1 ART 2 ART 3 ART 4 ART 5 ART 6 


EINBOFLOSSING МА 14  ITERATIES 
TABLEAU ZONDER DE KOLOMMEN ПЕК ARTIFICIAL VARIABELEN 


KASIS 
VARIA NIET -RASISVARIABELEN 
BELEN 
01 02 u2 03 81702 DEA3 
MAC 1 0.00 -1.25 1.25 0.00 -.06 -.03 
МАС2 0.00 1.25 -1.25 0.00 .06 .15 
DER2 0.00 0.00 0.00 0.00 1.33 2.00 
LAB1 0.00 -1.00 1.00 0.00 -.05 -.03 
LAR? 0.00 1.00 -1.00 0.00 0.00 0,00 
LAR3 0.00 0.00 0.00 0,00 .05 .08 
STO1 0.00 -25.00 25.00 0.00 -1.28 wl KE 
DER1 0.00 12.50 -12.50 0.00 .64 .33 
STO3 0.00 0.00 0.00 0.00 0.00 1.50 
01 -1.00 -1.00 1.00 0.00 -.05 - 03 
DEA2 0.00 0.00 0.00 0.00 -1,11 -2.67 
В3 0.00 0.00 0.00 0.00 -1.33 -2.00 
^1 0.00 0.00 0.00 0.00 0.00 0.00 
A2 0.00 0.00 0.00 0.00 -1.11 -2.67 
АЗ 0.00 0.00 0.00 0.00 1.11 1.67 
Ві 0.00 12.50 -12.50 0.00 .64 33 
B2 0.00 -12.50 12.50 0.00 e 69 1.67 
03 0.00 0.00 0.00 -1.00 -.05 -.08 
DOELF 15.00 12.38 2.63 15.00 66 1.69 


VERVOLG TAELEAU 


BASIS 
'ARIA NIET-BASISUVARIABELEN 
BELEN 
DEA1 MAC3 DERZ WAARDE BASIS 
МАС1 .06 -.28 -.03 474 .53 
МАС2 -.06 1.28. .13 147.47 
DER2 0.00 16.67 1.67 1733.33 
LAR1 0.00 -.22 -.02 237.22 
LAR2 0.00 0.00 0.00 150.00 
LAEZ 0.00 .22 .10 26.78 
5701 .25 -5.56 -.56 3130.56 
DER1 ‚63 2.78 .28 2434.72 
STO3 0.00 0.00 2.00 5600.00 
01 0.00 -,22 -.02 62.22 
ПЕА2 0.00 -22.22 -2.22 3022.22 
B3 0.00 -16.67 -2.67 3866.67 
А1 -1.00 0.00 0.00 100.00 
^2 1.00 -22.22 -2.22 5522.22 
АЗ 0.00 22.22 2.22 3777.78 
В1 63 2.78 28 3534.72 
B2 7-63 13.89 1.39 298.61 
03 0,00 -.22 -.10 148.22 
SAL. поз -9621,847222 


108 


PROJECT NAME 


L 


INPRO 


5 
N 
Q Q Q ә 


MAC3 


LAB] 


LAB2 


LAB3 


.08 


5Т01 


5Т02 


STO3 


ОУЕ1 


ОУЕ2 


Q е e e өм e б 


ОУЕЗ 
.08 


. Hes Hoo a. 


Q 
ui 
Q 


BH ef ee 


Q 


> 
N 
Q 


. 
ui 
ө 


% 
Ut 
P> 
° 
ui 


. 
UY 
= 
$ 
Ut 


Q 
ui 
Q 


A3 
03 


Q е 


ві B2 
01 02 
„63 „42 
0 0 
Ві B2 
01 02 
el H 
0 0 
0 21 
0 0 
0 0 
0 0 
.08 0 
0 0 
0 .08 
0 0 
0 0 
0 0 
2 0 
0 0 
2 2 
0 0 
2 2 
0 0 
.08 0 
1 0 
0 .08 
0 1 
0 0 
0 0 
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РЕА] 1 0 0 0 0 
0 0 0 0 0 0 
H GE 100 

DEA2 1 1 0 0 0 
H 0 H H H H 
0 СЕ 2600 

DEA3 1 1 1 0 H 
0 0 0 0 0 0 
0 СЕ 9400 | 

DEB] 0 0 0 1 0 
0 0 0 0 0 H 
H GE 1100 

DEB2 0 0 0 1 1 
0 0 0 0 0 0 
0 СЕ 2100 

DEB3 0 0 0 1 1 
1 0 0 0 0 0 
0 СЕ 7700 


PROBLEM SOLUTION 


All items not listed below have value 0. 


Type of optimization: MINIMIZE 
OPTIMAL SOLUTION: GOAL- 9621.848 
DECISION VARIABLE OPTIMAL VALUE 
Al 99.99988 
A2 5522.223 
A3 3777.777 
ві 3534.722 
B2 298.6109 
B3 3866.667 
03 148.2223 
01 62.22226 


SLACK AND SURPLUS IN CONSTRAINTS 


TYPE RESOURCE VALUE 
SLACK МАСІ 474.5278 
SLACK MAC2 147.4721 
SLACK LABl 237.2223 
SLACK LAB2 150 
SLACK LAB3 26. 77776 
SLACK STO1 3130.557 
SLACK STO3 5600 


SHADOW PRICES FOR CONSTRAINTS 


RESOURCE VALUE 
_MAC3 -5.472223 

STO2 -.6636111 

ОУЕ2 -2.625001 

ОУЕЗ -15 

DEA1 -.14875 

DEA3 -1.686667 


DEB3 -1.957222 
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REDUCED COSTS FOR DECISON VARIABLES 


VARIABLE COST 

01 15 

02 12:375 
02 2.625001 
03 15 


RANGES ОҒ OBJECTIVE FUNCTION COEFFICIENTS 


VARIABLE LOWER BOUND PRESENT VALUE UPPER BOUND 
. 69125 .84 UNBOUNDED 
„31375 ‚56 .70875 

-880803.9 .28 .5262501 
.42 „63 „8679999 
„1820001 | .42 „6300001 

-.1183334 .21 1761609 
9.536743E-07 15 UNBOUNDED 
2.625001 15 UNBOUNDED 
2.016303 15 5.637148Е%07 

-12.375 0 2.625001 

-2.625001 0 UNBOUNDED 

-15 0 UNBOUNDED 


RANGE OF RESOURCES 


RESOURCE LOWER BOUND PRESENT VALUE UPPER BOUND 
МАСІ 365.4725 840 UNBOUNDED 
MAC2 692.5275 840 UNBOUNDED 
MAC 3 818.5003 840 976.0001 
LAB] 287.778 525 UNBOUNDED 
LAB2 300 450 UNBOUNDED 
LAB3 498.2222 525 UNBOUNDED 
5Т01 7219.446 10350 UNBOUNDED 
STO2 15670.01 16100 17317.39 
STO3 29500 35100 UNBOUNDED 
ОУЕ1 287.778 350 UNBOUNDED 
OVE2 276.1115 300 417.978 
ОУЕЗ UNBOUNDED 350 498.2222 
ОБАЇ 0 100 3995.557 
DEA2 UNBOUNDED 2600 5622.223 
DEA3 8266.666 9400 9579.164 
DEB] UNBOUNDED 1100 3534.724 
DEB2 UNBOUNDED 2100 3833.334 


DEB3 6340 7700 7914.997 
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spelingsvariabelen LAB 1, LAB 2 en LAB 3 te beschouwen. Deze 
komen alle in de basis voor, hetgeen inhoudt dat deze restricties 
niet bindend zijn. LAB ] = 237,22; in het rechterlid van ге- 
strictie LAB 1 stond 525, zodat het aantal nuttig gebruikte uren 
in periode | gelijk is aan 525 - 237,22 - 287,78. In de normale 
werktijd zijn dus nog 350 - 287,78 - 62,22 uren niet gebruikt, 
hetgeen wij ook reeds weten uit Ui = 62,22. 


Eveneens niet bindend zijn de restricties MAC 1, MAC 23, STO 1 
en STO 3 (de bijbehorende spelingsvariabelen komen immers in de 
basis voor). In periode | zijn nog 474,53 machine-uren beschik- 
baar, terwijl aan het eind van periode 1 nog 3130,56 m? opslag- 
ruimte beschikbaar is. 


Wel bindend zijn de restricties STO 2 en MAC 3. Aan het eind 
van periode 2 is de opslagcapaciteit volledig benut. Zou er één 


m^ extra opslagruimte beschikbaar zijn, dan zouden de kosten 
met f 0,66 dalen van f 9621,85 tot f 9621,19. 


9.4. GEVOELIGHEIDSANALYSE VAN HET RECHTERLID 


Stel nu, dat de onderneming de mogelijkheden om voor periode 2 
magazijnruimte bij te huren wil verkennen en dat de volgende 
magazijnen te huur worden aangeboden: 


Magazijn Oppervlakte in m? Huurprijs voor één periode 
A 600 f. 300,.—- 
B 900 f 540,-- 
C 1800 7. 990,-- 


Is het verstandig één of meer van deze magazijnen te huren? Zo 
ja, welke? 

Om deze vraag te kunnen beantwoorden moeten wij weten binnen wel- 
ke grenzen de schaduwprijs van f 0,66 per m? geldig blijft. In de 
output van OPTIMIZER kunnen wij deze grenzen direct aflezen onder 
het kopje "RANGES OF RESOURCES". Wij vinden daar bij STO2 dat de 
schaduwprijs geldig blijft zolang 15670,01 « bg « 17317,39, waar- 
bij bg de coëfficiënt in het rechterlid van restrictie STO2 voor- 
stelt. De huidige waarde van bg is 16100, zodat de schaduwprijs 
geldig blijft zolang -429,99 « Abg < 1217,39. 


Als wij alleen beschikken over de output van IMOPTI, dan moeten 
wij zelf een gevoeligheidsanalyse van het rechterlid van restric- 
tie STO2 uitvoeren met behulp van de gegevens van het laatste 
simplex-tableau. Wij vinden dan, dat de schaduwprijs van 66 cent 
per m^ geldig blijft zolang geldt: 
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-147,47 -1733,33 -26,78 -2434,72 -3777,78 


mex( Gee 9679 6005. 10» уві 57 ATE ee 

-3534,72 -298,61 ЖАТТАУ ° 07 
DE. . t р 0 ké My «c mal erp rape 
-3130,56 -62,22  -3022,22  -3866,67 -5522,22 
z1798 ГО" ТТИ TR о РА , 
-148,22 , 
20,05 
ofwel 


-432,77 < Abg < 1244,4 


Dit interval wijkt enigszins af van het door OPTIMIZER gegeven 
interval. De afwijking wordt veroorzaakt door afrondingsfouten 
(de met behulp van IMOPTI berekende bovengrens van 1244,4 is 


met name doordat het getal in de noemer‏ ; و 
.0‚ 


is afgerond op slechts één significante decimaal kunnen grote 
afrondingsfouten optreden). 


berekend uit 


Wij concluderen dus, dat indien men magazijn A of B huurt men 
kan rekenen op een kostenverlaging van 600 respectievelijk 900 
maal 66 cent. Over de kostenverlaging die kan worden bereikt 

als men magazijn C of een combinatie van meer dan één magazijn 
huurt kunnen wij geen uitspraak doen. De netto besparingen bij 
het huren van magazijn A of B kunnen wij nu als volgt berekenen: 


Magazijn A Magazijn B 
Kostenvermindering 600 ж 0,66 - 396 900 ж 0,66 - 594 
Meerkosten door magazijnhuur 300 540 
Netto besparing 96 54 


Wij concluderen dus, dat men indien men gedwongen zou zijn te 
kiezen tussen A en B moet kiezen voor het huren van magazijn А. 
Mogelijk levert het huren van А Әп B of van magazijn С of van 
een andere combinatie van meer dan één magazijn een grotere 
netto besparing ор dan f 96,-. Als wij dit willen nagaan zullen 
wij het probleem nog een of meer keer met de simplex-methode 
moeten oplossen voor andere waarden van bg. 
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9.5. | GEVOELIGHEIDSANALYSE VAN DE COEFFICIENTEN VAN DE 
DOELSTELLINGSFUNCTIE 


In de optimale oplossing wordt de piek in de vraag naar beide 
produkten voor een deel opgevangen door overwerk (in periode 3) 
en voor een ander deel door in het voren (in periode 2 en in 
periode 1) te produceren. Blijkbaar is de som van overwerkkosten 
en voorraadkosten op die manier het kleinst. 


Stel nu eens, dat de overwerkkosten in periode 3 niet f 15,* maat 
slechts f 5,- per uur bedragen, terwijl alle andere gegevens 
gelijk blijven. Moet men dan in periode 3 meer gaan overwerken? 

Om deze vraag te beantwoorden moet men weten binnen welke grenzen 
со (dat is de coëfficiënt van 03 in de doelstellingsfunctie) kan 
variëren zonder dat de optimale oplossing verandert. In de output 
van OPTIMIZER vinden wij direct het antwoord. De gezochte grenzen 
zijn: 2.016303 < Cg < 5.637148 х 10°. De benedengrens is dus afge- 
rond f 2,02 en de bovengrens ruim f 56 miljoen per uur. Deze 
bovengrens is dus praktisch oneindig groot. 


Beschikken wij alleen over de output van IMOPTI, dan moeten wij 
zelf de gevoeligheidsanalyse uitvoeren met behulp van de gegevens 
uit het laatste simplex-tableau. Wij vinden dan**): 


С ИБ 1859. NA 1.96 
max( — , ——. ت‎ 


min с» 
0 9 


ofwel TIe ж Bsp < = 


ofwel Ke cg < = 


Door afrondingsfouten verschilt dit interval enigszins van het 
door OPTIMIZER gegeven interval. 


De conclusie is echter duidelijk: als de overwerkkosten in de 
derde periode f 5,- per uur bedragen in plaats van f 15,- per 
 uur blijft de optimale oplossing ongewijzigd. 


9.6. VRAAGSTUKKEN 


9-1. Caravanfabriek Van der Woude B.V. heeft op een gegeven 
moment 200 werknemers in de direct productieve afdelingen in 
dienst. De vraag naar caravans vertoont een sterk seizoenspatroon. 
Voor de komende drie planningsperioden verwacht men de volgende 
vraag: 
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Planningsperiode Vraag (aantal caravans) 
1 600 
2 800 
3 2000 


Voor het makan van één caravan zijn 80 manuren nodig. Na aftrek 
van snipperdagen, ziekte, vakantie en dergelijke is het aantal 
per man, per periode beschikbare uren: 


Periode Aantal uren beschikbaar per man 
1 580 
2 580 
3 500 


De directe loonkosten per man bedragen f 10.000,- per periode. 
Het is mogelijk het aantal direct-produktieve personeelsleden 

te veranderen aan het begin van iedere periode. Een verminde- 
ring van het aantal direct produktieve personeelsleden veroor- 
zaakt afvloeiingskosten van f 15.000,- per af te vloeien perso- 
neelslid. Een vergroting van het aantal personeelsleden veroor- 
zaakt wervingskosten van f 1000,- per te werven personeelslid, 
Van der Woude B.V. kan een ongelimiteerd aantal caravans op- 
slaan. De voorraadkosten bedragen f 180,- per caravan per pe- 
riode. De beginvoorraad is 200 caravans. Aan het eind van de 
derde periode wenst men tenminste een voorraad van 200 caravans 
te hebben. 

De produktieplanner van van der Woude B.V. wenst een zodanig 
produktieplan op te stellen, dat de som van loonkosten, voorraad- 
kosten, afvloeiingskosten en wervingskosten zo klein mogelijk is. 


a. Formuleer dit probleem als een LP-probleem. 
b. Los het probleem op met gebruikmaking van IMOPTI of van 
OPTIMIZER. 


9-2. Het probleem van Van der Woude B.V. is in vraagstuk 9-1. 
sterk vereenvoudigd weergegeven. In werkelijkheid produceert 
van der Woude B.V. naast caravans namelijk ook directieketen. 
Voor de komende 3 perioden is de vraag naar caravans en direc- 


tieketen: 
periode caravans directieketen 
1 600 725 
2 800 750 


3 2000 800 


115 


Voor het maken van een directiekeet zijn 60 manuren nodig (voor 
een caravan 80 uren). 

Aan het begin van iedere periode kunnen personeelsleden worden 
ontslagen of aangenomen. Ontslag van één personeelslid kost 

f 15.000,-, het werven van een personeelslid kost 7 1000,-. De 
directe loonkosten bedragen 7 10.000,- per personeelslid per 
periode. 

Het is mogelijk de directe personeelsleden te laten overwerken. 
Het maximum daarbij is 202 van het in een bepaalde periode in de 
normale werktijd beschikbare aantal manuren. De kosten van 
overwerk bedragen 7 30,- per uur. 

Het aantal werkuren per periode gedurende normale werktijd is 
gegeven in onderstaande tabel: 


periode aantal uren beschikbaar per man 
1 580 
2 580 
3 500 


Het aantal personeelsleden is thans 267. Gezien de in de fabriek 
beschikbare ruimte en de aanwezige gereedschappen en machines 
kunnen maximaal 300 personeelsleden tegelijkertijd in de fabriek 
werken, Teneinde in de zomermaanden een zo groot mogelijke pro- 
duktie te kunnen bereiken is in voorgaande jaren gedurende en- 
kele maanden in 2-ploegendienst gewerkt. De kosten van de ploe- 
gentoeslag bedragen 254 van de normale loonkosten. Uit ervaring 
weet men, dat 2/3 deel van het aantal personeelsleden niet be- 
reid is in ploegendienst te werken. Personeelsleden die gedu- 
rende een bepaalde maand in ploegendienst werken mogen geen 
overwerk verrichten. 

De voorraadkosten bedragen f 180,- per periode voor een caravan 
en f 160,- per periode voor een directiekeet. De beginvoorraad 
is 200 caravans en 100 directieketen. Men wenst aan het eind 

van periode 3 tenminste 200 caravans en 100 directieketen in 
voorraad te hebben. | 


De produktieplanner van van der Woude B.V. staat nu voor de 

taak een produktieplan te ontwerpen dat zodanig is dat 1€ steeds 
aan de vraag kan worden voldaan en 29 de som van normale loon- 
kosten, overwerkkosten, ploegentoeslag, afvloeiings- en wervings- 
kosten en voorraadkosten zo klein mogelijk is. 


a. Geef aan welke beslissingsvariabelen u zou willen gebruiken 
in het voor dit probleem op te stellen LP-model. Geef iedere 
beslissingsvariabele een naam. 


b. Schrijf de doelstellingsfunctie van het model volledig uit. 
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с. Geef aan welke restricties u bij dit LP-model zou willen 
onderscheiden. Geef iedere restrictie een naam. Omschrijf 
bij iedere restrictie of groep van restricties de betekenis 
van de restrictie(s) met enkele woorden. 


d. Schrijf de restricties volledig uit. 


9-3. Neduco B.V. is een producent van luierbroekjes. Deze luier- 
broekjes worden afgezet in vier geografische gebieden (Benelux, 
Duitsland, Frankrijk en Italië). 

De vraag per periode in deze vier gebieden is voor de komende 
drie perioden (uitgedrukt in 1000 dozen): 


periode 
1 2 3 
Benelux 200 240 280 
Duitsland 300 320 340 
Frankrijk 120 140 180 
Italië 80 80 90 


Neduco beschikt over fabrieken in Nederland en in Duitsland. In 
beide fabrieken wordt per periode gedurende 160 uur gewerkt in 
normale werktijd; in periode 3 is in de Duitse fabriek echter 
slechts 100 uur beschikbaar in verband met het begin van de 
bedrijfsvakantie. In de Nederlandse fabriek produceert men 2000 
dozen per uur; in de Duitse fabriek is dit 3000 dozen per uur. 
In beide fabrieken is overwerk mogelijk tot ten hoogste 25 рго- 
cent van de normale werktijd. De kosten van overwerk bedragen 
voor de Nederlandse fabriek f 2000,-- per uur, voor de Duitse 
fabriek 7 2800,-- рег uur. 

In Frankrijk en Italië beschikt Neduco over een depot, іп Neder- 
land en in Duitsland is er een magazijn voor gereed produkt op 
het fabrieksterrein. De in een bepaalde periode geproduceerde 
hoeveelheden worden nog gedurende dezelfde periode naar één van 
de depots òf naar één van de fabrieksmagazijnen gebracht. De 
opslagruimte in de depots en in de fabrieksmagazijnen bedraagt 
(in 1000 dozen): 


Nederland 500 
Duitsland 600 
Frankrijk 200 
Italië 100 


De opslagkosten per periode zijn gegeven in onderstaande tabel 
(bedragen per 1000 dozen): 


Nederland 7 Мы” 
Duitsland 1 4D 
Frankrijk f 24677 
Italië 7 Taa 
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De transportkosten (in guldens per 1000 dozen) zijn gegeven іп 


onderstaande tabel: 


Van: Nederland Duitsland 
Nederland 0 200, 
"ed Duitsland 14709. о 0 
aar: Frankrijk 1600 ,-- 1500,-- 
Italië 4000,-- 3400. == 


Пе beginvoorraad aan het begin van periode 1 en de minimumvoor- 
raad (dat is de voorraad die aan het einde van iedere periode in 
ieder depot c.q. magazijn aanwezig moet zijn in verband met 
onvoorspelbare vraagfluctuaties) zijn (in 1000 dozen): 


beginvoorraad minimum voorraad 
Nederland 140 120 
Duitsland 150 150 
Frankrijk 100 80 
Italië 60 50 


De produktieplanner van Neduco B.V. heeft tot taak een produktie- 
en transportplan op te stellen zó, dat de voorraad in geen van 

de depots of magazijnen ooit barada de minimum voorraad daalt, 

en zó dat de som van voorraadkosten, transportkosten en overwe 
kosten zo klein mogelijk is. 


a. Geef aan welke beslissingsvariabelen u zou willen gebruiken in 
het voor dit probleem op te stellen LP-model. Geef iedere 
beslissingsvariabele een naam. 

b. Schrijf de doelstellingsfunctie van het model volledig uit. 

c. Geef aan welke soorten restricties u in dit model wilt opnemen. 
Geef iedere restrictie een naam. 

d. Schrijf de restricties volledig uit. 


10. НЕТ OPSTELLEN VAN EEN 
STRUCTUURPLAN VOOR EEN 
GEMEENTE MET BEHULP VAN 
LINEAIRE PROGRAMMERING 


10.1. INLEIDING 


In het voorgaande hoofdstuk beschreven wij een toepassing van 
lineaire programmering op een probleem van operationele aard. 

In dit hoofdstuk beschrijven wij een toepassing van lineaire 
programmering op een probleem van strategische aard, nl. het 
opstellen van een structuurplan door een gemeente. Wij beschrij- 
ven deze toepassing aan de hand van een voorbeeld 5-2; 


10.2. PROBLEEMSTELLING 


Gegeven is onderstaande plattegrond van een plattelandsgemeente: 


opp. in ha.: 80 120 110 60 


Zone 1 en zone 4 zijn weide-gebieden; in deze zones is het ge- 
hele grondoppervlak in principe te gebruiken voor woningbouw. 
Zone 2 bestaat uit een zeer aantrekkelijk bos- en heidegebied; 
woningbouw in zone 2 is wel mogelijk, maar zal dan ten koste 
gaan van een waardevol natuurgebied. Men wenst dit natuurgebied 
zo veel mogelijk te sparen. Zone 3 is de bestaande dorpskern. 
Van zone 3 is thans 70 ha. bebouwd met woningen en centrumvoor- 
zieningen. 


Men wil in deze gemeente voor 5000 gezinnen woningen bouwen. 

Er worden twee typen woningen onderscheiden, nl. gesubsidieerde 
woningen en zgn. vrije-sector woningen. De vraag is nu het aan- 
tal te bouwen woningen per type en per zone te bepalen. De vol- 
gende gegevens zijn verder nog bekend: 
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- Een gesubsidieerde woning vraagt een grond-oppervlak van 400m?, 
een vrije sector woning 800m?. 

- in zone 3 bevinden zich 400 enigszins vervallen woningen, die 
eventueel voor sloop in aanmerking kunnen komen. Per te slopen 
woning komt dan eventueel 1000m? grondoppervlak vrij. Voor 
iedere te slopen woning moet uiteraard wel een nieuwe woning 
worden gebouwd, eventueel in een andere zone. 

- In zone 3 bevinden zich thans de centrumvoorzieningen die zijn 
afgestemd op het bestaande aantal woningen. Ieder nieuw gezin, 
dat zich in de gemeente komt vestigen heeft centrumvoorzienin- 
gen nodig. Deze centrumvoorzieningen zullen alle in zone 3 
moeten worden gerealiseerd. Per gezin is daarvoor 50m? grond- 
oppervlak nodig. 

- Voor het totaal aantal nieuw te bouwen woningen moet 2/3 deel 
in de gesubsidieerde sector en 1/3 deel in de vrije sector 
worden gebouwd. 

- Zowel om redenen van architectonische als van sociologische 
aard wenst men in de zones 1, 2 en 4 tenminste 252 gesubsidi- 
eerde woningen en tevens tenminste 257 vrije-sector woningen 
te bouwen. 

- De "wellfare-coefficients" zijn als volgt !?); 


zone gesub. woning vrije-sector woning 
l 1,0 4,9 

2 2,4 6,6 

3 4,5 8,5 

4 3:4 7,0 


Afbraak van woningen: -80,0 per woning. 

- Men wenst een zodanig plan te ontwerpen, dat enerzijds de 
waarde van de wellfare-function zo groot mogelijk is en ander- 
zijds een zo groot mogelijk deel van het natuurgebied in zone 
2 wordt gespaard. 


10.3. DE FORMULERING VAN HET MODEL 
Wij onderscheiden 9 beslissingsvariabelen, nl. 


Dij = het aantal te bouwen woningen van type і in zone і. 
j = l slaat op een gesubsidieerde woning, j = 2 op een 
vrije sector woning. Door beide woningtypen in elk van de 
vier zones te bouwen, hebben wij hiermee 8 beslissingsva- 
riabelen gedefinieerd. 


120 


DM = het aantal іп zone 3 af te breken woningen. 


Wij formuleren nu eerst de doelstellingsfunctie. Wij willen ma- 
ximeren 

В11 » 4,9D1» jd 2,4021 т 6,6022 M 4,5031 + 8,5032 + 3,1953 

+ TOD: 80DM 


Wij formuleren nu de restricties. 


Voor iedere zone geldt een restrictie ten aanzien van het maxi- 
male grondgebruik. Voor de zone 1 is dat: 


4D11 + 8012 < 8000 


waarbij als eenheid van oppervlakte de are is gekozen. 
De vier restricties ten aanzien van het grondgebruik luiden 


ZON. 12:76011 7 8012 < 8000 
ZON 2: 40 21 + 8022 « 12000 


A 


ZON 4: АВ + 902 « 6000 


In zone 3 is іп totaal 11000 аге beschikbaar. Dit beschikbare 
oppervlak moet worden verminderd met het reeds bebouwde gedeelte 
(7000 are) en met de grond, die voor voorzieningen nodig is 
(5000 * 0,5 = 2500 are). In het rechterlid van restrictie ZON 3 
staat daarom 11000 - 7000 - 2500 = 1500 are). 


In zone 3 mogen maximaal 400 woningen worden gesloopt. 
Dit geven wij aan in de volgende restrictie: 


MDEM: DM < 400 


Er moeten in totaal 5000 woningen worden gebouwd. Dit aantal 
moet worden vermeerderd met het aantal te slopen woningen: 


ТИМЕ: D+ Dj2* Dat Daat Dar D32* Duit Du27 DM = 5000 
Een derde deel van het aantal te bouwen woningen moet in de 
vrije sector worden gerealiseerd en tweederde deel in de gesub- 


sidieerde sector: 


DWTY: Diit Dar 031+ Dui = 2(012+ 022+ 032+ Duz) 
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In zone 1 moet van het totaal aantal te bouwen мопіпреп ten- 
minste 257 in de gesubsidieerde sector worden gebouwd: 


АВС 1: Din E 0,25(D1 71+ 012) 
ре vijf andere restricties van dit type zijn: 


ARC 2: Dj 2 0,25(021% D33) 


ARC 3: Dur 20,25(р,1% Dus) 
ARC 4: Di2 2 0,25(011% Dj») 
ARC 5: D22 > 0,25(021% Daa 
ARC 6: Dya > 0,25(Duj+ Du) 


Het volledige model luidt als volgt: 
^ | maximeer 


5» 4,912 + 2,4021 + 6,6022 + 4,5D31 * 8,5032 
+ 3,10, + 7,0,2 - 80DM 


onder de voorwaarden 


ZON 1: 4D11 * 8D1» « 8000 
ZON 2: 4D» IX 802 2 « 12000 
ZON 3: 4D31 + 8032 - JODM < 1500 


ZON 4: 4D, + 80» < 6000 
МРЕМ : DM < 400 

ARC 1: -0,75011 + 0,25012 < 0 
ARC 2: -0,75021 + 0,250, < 0 
ARC 3: -0,750,1 + 0,250,2 < Q 
ARC 4: 0,25011 - 0,75012 < 0 
ARC 5: 0,2501 - 0,7502 < 0 
ARC 6: 0,250, - 0,750,2 < 0 


ТИМЕ : рут+ Пі2% П21% Daat D31* Daat Різ Dua DM - 5000 
DWTY : раче 2D1 2+ Danz 2022% Ра” 20323 Ріг 20,2 0 
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10.4. ПЕ OPLOSSING VAN НЕТ MODEL 


Volgens de іп paragraaf 10.2. gegeven probleemstelling wenst 
men enerzijds de welfare-function te maximeren en anderzijds 
een zo groot mogelijk natuurgebied in zone 2 te sparen. Dit is 
dus een probleem met twee doelstellingen. Noemen wij de waarde 
van de welfare-function z en het aantal te sparen are natuur” 
gebied in zone 2 у, dan zijn wij dus met name geïnteresseerd 
in z als functie van y. 

Wij lossen het probleem daarom een aantal malen op voor уег- 
schillende waarden van у, namelijk 


y = 0, 500, 1000, 2000 еп 2650 are 


Іп het rechterlid van restrictie 2 staan dan achtereenvolgens 
de getallen 12000, 11500, 11000, 10000 en 9350. | 
De input en de output van ІМОРТІ voor deze 5 alternatieven is 
gegeven op de pagina's 109-118. 


In tabel 10.1. hebben wij de waarden van у, 2 еп up Voor de vijf 
alternatieven samengevat. Met u, is bedoeld de schaduwprijs van 
het natuurgebied in zone 2 uitgedrukt in eenheden van de wel- 
fare-function per are. Deze schaduwprijs vinden wij in de rij 
van de doelstellingsfunctie in de kolom van spelingsvariabele 


ZON 2. 
Alternatief У 2 02 
А 0 18872,50 0,30 
В 500 18725,00 0,30 
С 1000 16115,00 15,07 
D 2000 942,86 15,20 2 
E 2650 -8939,46 15,20 


- Tabel 10.1. Aantal аге te sparen natuurgebied in zone 
2 (у), waarde van de welfare-function (2) en schaduwprijs 
van één are natuurgebied in eenheden van de welfare- 
function (u>) - 


In tabel 10.2. geven wij een overzicht van de basisvariabelen 
behorend bij de vijf alternatieven. 


123 


L INEATR-FROÜGRANMERINGSEROGRAM ІМОРТІ 
WILT U TE GEGEVENS VIA TAPE ІМІЕ?ЕМ?ТҮРЕ JA OF NEE? 
CM LUA*1 = 25022В» LOADER USED  41400BNEE 


GEEF ПЕ VOLGENDE SPECIFICATIE VAN HET PROBLEEM? 
AANTAL BESLISSINGSVARIARELEN ? 9 


TOTAAL. AANTAL RESTRICTIES #13 
AANTAL. =< RESTRICTIES ? 11 
AANTAL. = RESTRICTIES ? 2 
AANTAL. >= RESTRICTIES ? o 


MOET ПЕ LOELSTELLINGSFUNKTIE GEMAXIMEERD OF GEMINIMEERI WORDEN 
ТҮРЕ МАХ OF MIN ? МАХ 


GEEF ПЕ 9 NAMEN VAN ПЕ BESLISSINGSUARIABELEN? 
011 »D12 »D21 ,022 „031 +032 »I41 ,п42 ,рм 


GEEF DE БІЛБЕНОКЕМПЕ COEFFICIENTEN VAN DE DOELSTELLINGSFUNKTIE IN 
DEZELFDE VOLGORDE? 

1,4,9,2.4,6,6,4,5,8,5,3,1,7,-80 
GEEF ACHTEREENVOLGENS VOOR ELK VAN DE 13 RESTRICTIES EEN REGEL MET DE 
NAAM VAN ПЕ RESTRICTIE GEVOLGE DOOR ПЕ TECHNISCHE COEFFICIENTEN EN HET 
КЕСНТЕКІ INS 
ПЕ 11 == RESTRICTIES 

1720М1/4»,8,0»0»0»0»0»0,0»8000 

2?220N2/0»0»4»8»0»0»0»0»0»12000 

JTZO0N3/0»0»0»0»4»8»0»0»-10»1500 

4?020N4/0»0»0»0»0»0»4»8»0» 6000 

О?МПЕМ/0,0,0,0,0,0,0,0,1,400 

Ф"АКСІі/-,75,.25,0,0,0,0,0,0,0,0 

2?АКС2/0,0»-.25,.25,0,0»0»0,0»0 

ВТАКСЗ/0,0»,0»0»0»0»-.275»,.,25,0,0 

??РАКС4/,25,-.,75,0,0,0,0,0»0,0»0 

ІОРАКС5/0,0,,25,-,75,0,0,0,0,0,0 

11?7ARC6/0»0»0»0»0»0».25»-,75»0»0 


DE 2 - RESTRICTIES 
1?TTDIWE/1»1»1»1»1»1»1»1»-1» 5000 


2?DUTY/1»-2»1»-2»1»-2»1»—-2»050 
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ECHTE VARIABELEN ^ ^" 011 bi12 02111022: D31 2132 DAP CDAS DM 

SLACK VARIABELEN 70мі ZON2 ZON3 ZON4 МПЕМ ARCI ARC2 ARCI ARCA 
ARCS ARC 

ARTIFICIAL VARIABELEN ТОМЕ DUTY 


EINDOFLOSSING МА 13 ITERATIES 
TABLEAU ZONDER DE KOLOMMEN ПЕК ARTIFICIAL VARIABELEN 


BASIS 


VARIA NIET-RASISVARIARBELEN 
BELEN 
20м4 ARCO 20м2 n32 20М3 ARC 5 
п21 0.00 0.00 .15 0.00 0.00 1.60 
ARC2 0.00 0.00 .10 0.00 0.00 1.40 
П31 0.00 0.00 0.00 2.00 .25 0.00 
20М1 1.00 0.00 1.00 0.00 1.00 0.00 
МПЕМ 0.00 0.00 0.00 0.00 0.00 0.00 
ARC1 -.10 -1.40 10 -1.75 ^. 19 -1.40 
п22 0.00 0.00 .05 0.00 0.00 -.80 
п42 „05 -.80 0.00 0.00 , 0.00 0.00 
П11 -.15 -1.60 -.15 -2.00 -.25 -1.60 
141 .15 1.60 0.00 0.00 0.00 0.00 
АКСА 0.00 1.00 0.00 1.25 .06 1.00 
ARCI .10 1.40 0.00 0.00 0.00 0.00 
пі2 -. 05 80 "05 1.00 0.00 .80 
HOELF .42 1.68 . 30 3.40 . 87 :88 


VERVOLG TABLEAU 


BASIS 
VARTA NIET--RASISVARTABELEN 
BELEN 
ПМ WAARDE BASIS 

n21 0.00 1800.00 
ARC? 0.00 1200.00 

31 -2,50 7: 0 

ZONI -4.67 833.33 

MDEM 1.00 400.00 

ARCI 1.46 2.08 

n22 0.00 600.00 

142 0.00 300.00 

D11 1.83 2580.33 

141 0.00 900.00 

ARCA "6071 510,42 

ARC 3 0.00 600.00 

п12 -33 766.67 
DOELF 68.95 18872.50 


SAL ПО 18872,500000 
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SIU 
COMMAND- ІМӘРІІ 


LINEAIK--FROGRAMMERINGSFROGRAM IMOFTI 


WILT U ПЕ GEGEVENS VIA TAFE ІМІ.Е?ЕМ?ТҮРЕ JA OF NEE ? 
CM LWUA*1 = 250776» LOALER USED 41 400 KNEE: 


TOTAAL. AANTAL RESTRICTIES ?13 
AANTAL. «є RESTRICTIES ? 11 


AANTAL. RESTRICTIES ? 2 


AANTAL “>= RESTRICTIES ? 0 


MOET ПЕ DOELSTELLINGSFUNKTIE GEMAXIMEERD OF GEMINIMEERI! WORDEN 
TYFE MAX OF MIN ? MAX 


GEEF ПЕ 9 NAMEN VAN ПЕ RESLISSINGSUARIABELEN'^? 
П11 »D12 »D21 ,П22 s031 ,п32 ,П41 ,П42 "ТІМ 


GEEF ПЕ BIJREHORENIE COEFFICIENTEN VAN DE NOEL STELLINGSFUNKTIE IN 
DEZELFDE VOLGORDE? 

1,4.9,2.4,6.6,4.5,8.,5,3.1,7»,-80 
GEEF ACHTEREENVOLGENS VOOR ELK VAN ПЕ 13 RESTRICTIES EEN REGEL MET ПЕ 
NAAM VAN ШЕ RESTRICTIE GEVOLGD LOOR DE TECHNISCHE COEFFICIENTEN EN HET 
RECHTERL LL 
ПЕ 11 =< RESTRICTIES 

1720М1/4,8,0,0,0,0,0,0,0,8000 

2?20м2/0,0,4,8,0,0,0,0,0»11500 

3?20N3/0»0»0»0»4»8»0»0»-10»1500 

4?20M4/0»0»0»0»0»0»4»8»0» 6000 

S?MDIEM/0O»Os0O»0»0»0»0»0»1»400 

6TARC1/-.75».25»0»0»0»0»0»0»050 

7?”АКС2/0»0»-.25».25,0»0»0»0»0»0 

В?АКСЗ/0»0»0»0»0»0»-.25».25,0»0 

9?АКС4/.25»-.2%,0»0»0»0»0»90»0»90 
10тАакс5/0,0,,25,-,75,0,0,0,0»0»0 

11?АКСЬ/0,0,0,0,0,0,,25у-,75,0,0 


ПЕ 2 = RESTRICTIES 
1TTDUE/1»1»1»1»1»1»1»1»—-1» 5000 


2?ПМТҮ/1»-2,1»-2,1»,-2,1»-2,0»0 
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ECHTE VARIABELEN Dii peas рої Dee: паз 022 -041 [42 ІМ 


SLACK VARIABELEN ZON1 ZON2 ZON3 ZOMA MEM ARCI ARC2 ARC3 ARCA 
ARCS ARCE 
ARTIFICIAL VARIABELEN TOWE DUTY 


EINDOFLOSSING NA 13  ITERATIES 
TARLEAU ZONDER ШЕ KOLOMMEN ПЕК ARTIFICIAL VARIABELEN 


RASIS 
VARTA NIET--RASISVARTABELEN 
BELEN 
20М4 ARCÓ ZON? 132 ZÜNS ARCI 
IER 0.00 0.00 .15 0.00 0.00 1.60 
акс2 0.00 0.00 .10 0.00 0.00 1.40 
DA) 0.00 0.00 0.00 2.00 «29 0.00 
ZON: 1.00 0.00 1.00 0.00 1.00 0.00 
МПЕМ 0.00 0.00 0.00 0.00 0.00 0.00 
ARC1 -.10 -1.40 -.10 "1.75 see ER -1.40 
n22 0.00 0.00 .О5 0.00 0.00 -.80 
п42 .05 -80 0.00 0.00 0.00 0.00 
111 -.15 "1.60 -.15 -2.009 -.25 "1.60 
П41 .15 1.60 0.00 0.00 0.00 0.00 
ARCA 0.00 1.00 0.00 1.25 .06 1.00 
ARC3 .10 1.40 0.00 0.00 0.00 0.00 
112 -.05 .80 -.05 1.00 0.00 e HO 
DOELF .42 1.68. e 30 5.40 :87 e 88 


VERVOLG TABLEAU 


BASIS 

VARIA NIET-RASISUARIARELEN 
BELEN 

DM WAARLE BASIS 

р21 0.00 1725.00 

ARCA 0.00 1150.00 

I31 -2.50 575.00 

ZON1 -4,67 553.33 

MDEM 1.00 400.00 

АКС1 1.46 52.08 

122 0.00. 575.00 

п42 0.00 500.00 

рії 1,83 333.33 

141 0.00 900.00 

ARCA sed) 510.42 

ARCI 0.00 600.00 

пі? -.33 791.67 
DOEL.F 68.95 18725.00 


SAL Dis 18725.000000 
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LINEAIR-FROGRAMMERINGSPROGRAM IMOFTI 
WILT U ПЕ GEGEVENS VIA ТАРЕ INLEZENPTYPE JA OF NEE 7 
CM LWAti = 25077R» LOALER USED 41400ВМЕЄЕ 


GEEF ПЕ VOLGENDE SPECIFICATIE VAN НЕТ PROBLEEM? 
AANTAL. BESLISSINGSVAKRIABELEN ? 9 


TOTAAL. AANTAL RESTRICTIES 213 
AANTAL. =< RESTRICTIES ? 11 
AANTAL = RESTRICTIES ? 2 
AANTAL. >= RESTRICTIES ? 0 


MOET DE DOELSTELLINGSFUNKTIE бЕМАХІМЕЕКІ OF GEMINIMEERD WORDEN 
ТҮРЕ МАХ OF MIN ? MAX 


GEEF ПЕ 9 | NAMEN VAN ШЕ BESLISSINGSVARIABELEN? 
Dii ,П12 »D21 »n22 ,ПЗі +032 »IDA41 ,»П42 iM 


GEEF ПЕ RIJREHORENIE COEFFICIENTEN VAN ПЕ DOELSTELLINGSFUNKTIE IN 
DE ZELEDE VOLGORDE? 

1,4.9,2.4,6.6,4.5»8.5,2,1,7»,-80 
GEEF ACHTEREENVOLGENS VOOR ELK VAN DE 13 RESTRICTIES EEN KEGEL MET DE 
NAAM VAN DE RESTRICTIE GEVOLGD DOOR DE TECHNISCHE COEFFICIENTEN EN HET 
RECHTERL ІП) 
ПЕ 11 «є RESTRICTIES 

1720М1/4»,8»,0»0»0»0»0»0»0,8000 

2770М2/0»,0,4,8»,0,0»0»0»0%11000 

3?2Z0N3/0»0»0»0»4»8»0»0»-10»1500 

4?20N4/0»0»0»0»0»0»4»8»0»6000 

О"ТМПЕМ/О,О0,0,0,0,0»0,0,1,400 

Ф?АКС1/-.25».25,0»0»0»0»0»0»0»0 

??АКС2/0»,0»-.2%8,.2%5,0,0»0»0»0У0 

gTARC3/0»0»0»0»0»0»-.753». 25»0»0 

9"AàÀRC4/.298»—-.75»0»0»0»0»0»0»05»0 

10?7АКС5/0»,0».25»-.275»,0»0»0»0»0»0 

117АКСь/0,0,0,0,0,0,,25,-,75,0,0 


DE 2 = RESTRICTIES 
17TIIE/1s1»1»1»1»1»1»1»—-1»5000 


2?DUTY/1»-2»1»-2»1»-2»1»-2»0»0 
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ECHTE VARIABELEN Dii Di2 021-022 .0581 09042 09204200 
SLACK VARIABELEN 70м1 ZÜN2 20М3 20М4 MDEM ARCI ARC2 ARCI ARCA 


ARCS АКСФ 
ARTIFICIAL VARIABELEN ТОМЕ ПМТҮ 


EINBOFLOSSING NA 13  ITERATIES: 
TABLEAU ZONDER DE KOLOMMEN ПЕК ARTIFICIAL VARTABELEN 


EASIS 
VARTA NIET-BASISVARIABELEN 
КЕШЕМ 
20М1 ARC ZÜN2 32 ZON3 ARCS 
П21 0.00 0.00 .15 0.00 0.00 1.60 
ARC2 0.00 0.00 .10 0.00 0.00 1.40 
151 -.34 0.00 -.34 2.00 ~, 29 0.00 
ARC A 0.00 1.40 0.00 0.00 0.00 0.00 
МПЕМ .21 0.00 .21 0.00 .21 0.00 
n22 0.00 0.00 .О5 0.00 0.00 "80 
42 0.00 -80 0.00 0.00 0.00 0.00 
п11 39 -1.60 .24 -2.00 .14 "1.60 
141 0.00 1.60 0.00 0.00 0.00 0.00 
ARCA "15 1.00 "449 1.25 ا‎ 1.00 
ПМ -.21 0.00 -.21 0.00 — 21 0.00 
п12 -.07 .80 wf 12 1.00 -.07 + 80 
DOELE 14.78 1.68 15.07 5.40 15.65 .88 


VERVOLG TABLEAU 


RASIS 
VARTA NIET-RASISVARIARELEN 
BELEN 
20м4 WAARDE BASIS 

П21 0.00 1650.00 

ARC2 0.00 1100.00 

n3i -.54 464.29 

. ARC3 .10 800.00 

МПЕМ .21 364.29 م‎ 
АКС1 .21 90.00 

Dn22 0.00 550.00 

п42 e 05 300.00 

111 .24 542.86 

ARCA “4195 535.71 

[M - 21 35,71 

012 "$2 828.57 
BOELF 15.20 16115.00 


SALDO= 16115.000000 
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LINEAIR-PROÜGRAMMERINGSPROGRAM IMOPTI 
WILT U ПЕ GEGEVENS VIA TAFE INLEZEN?TYPE JA OE NEE ^ 
CM LWAFL = 250776, LOANER USED 41400ВМЕЕ 


GEEF DE VOLGENDE SPECIFICATIE VAN HET PROBLEEM? 
AANTAL. BESLISSINGSVARIARELEN ? 9 


TOTAAL. AANTAL RESTRICTIES 713 
AANTAL. =< RESTRICTIES ? 11 
AANTAL = RESTRICTIES ? 2 
AANTAL >= RESTRICTIES ? 0 


MOET ПЕ DOELSTELLINGSFUNKTIE GEMAXIMEERD OF GEMINIMEERD WORDEN 
ТҮРЕ MAX OF MIN ? MAX 


GEEF ПЕ 9 NAMEN VAN ПЕ BESLISSINGSVARIABELEN? 
011 +D12 s021 »n22 ,ПІЗІ »D32 »I41 ,П42 e DN 


GEEF ПЕ BIJBEHORENDE COEFFICIENTEN VAN ПЕ ПОЕ STELL INGSFUNNTIE IN 
DEZELFDE VOLGORDE? 

1,4,9,2,4,6,6,4.5,8,5,3.1,7,-80 
GEEF ACHTEREENVOLGENS VOOR ELK VAN ПЕ 13 RESTRICTIES EEN REGEL МЕТ ПЕ 
NAAM VAN DE RESTRICTIE GEVOLGD LOOR DE TECHNISCHE COEFFICIENTEN EN HET 
RECHTERL II 

ПЕ 11 = RESTRICTIES 

1%20М1/4»,8»0У0»0»0»0»0»0,8000 

2?ZÜ0N2/0»0»4»8»0»0»0»0»0»10000 

3?20M3/0»0»0»0»4»8»0»0»-10»1500 

4?Z0N4A/0»0»0»0»0»0»45»8»0» 6000 

О?МПЕМ/О»0»0,0,0,0,0,0,1,400 

6Т”АКС1/-.25».25,0»0»,0»0,0»0»0»0 

72?АКС2/0,0»-.25».25»0,0,0»0»0»0 

втаксз/о0,0,0,0,0,0,-,75,,25,0,0 

9?АКС4/.25»-.25,0»0»0,0»0»0»0»0 

10?АКС5/0»,0».25,-.75,0,0»0»0»0»0 

117АКС6/0»0»0,0»0»90».25»-.75»0»0 


ПЕ 2 = RESTRICTIES 
17TTDUE/1»1»1»1»1»1»1»1»-1» 5000 


2?ПЫТҮ/1»-2,1»-2,1»-2,1»-2,0»0 
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ECHTE VARIABELEN D11 Di? D21 D22 D31 D32 DAL D42. DM 
SLACK VARIABELEN ZONI ZON ZOM3 ZONA MEM ARCI ARC2 ARC3 ARCA 
ARUS ARCÓ | 


ARTIFICIAL VARIABELEN ТІМЕ ПМТҮ 


EINDOFLOSSING NA 14 ITERATIES 
TABLEAU ZONDER ПЕ KOLOMMEN DER ARTIFICIAL VARIABELEN 


RASIS 

VARIA NIET-BASISVARIABELEN 

BELEN 

| 20М1 ARCO ZONA АКС1 ZOM3 32 
21 . 56 "1.60 .39 1.14 .14 -2.00 
ARC2 e 31 -1.40 e 51 1.00 .13 -1.79 
ARCS -.22 1.00 -.15 — Z éi -.09 1.25 
ARC3 0.00 | 1.40 0.00 0.00 0.00 0.00 
MODEM .21 0.00 .21 0.00 .21 0.00 
П31 -.54 0.00 -.94 0.00 "ому 2,00 
n22 -.18 .80 -.07 -.57 -.07 1.00 
DA? 0.00 - 80 0.00 0.00 0.00 0.00 
011 .04 0.00 . 0.00 -1.14 0,00 0.00 
n41 0.00 1.60 0.00 0.00 0.00 0.00 
ARCA .07 0.00 0.00 .71 0.00 0.00 
IM -.21 0.00 -.21 0.00 -.21 0.00 
D12 ell 0.00 0.00 «37 0.00 0.00 

DOELE 14.97 .80 15.20 .63 19.73 2.30 


VERVOLG TABLEAU 


RASIS 
VARIA NIET-RASISUARIABELEN 
REL.EN 
ZÜNA WAARDE BASIS 

n21 .24 1314.29 

ARC .21 837.50 

ARCS -.15 116.07 

ARCI .10 600,00 

МПЕМ ЖА 150.00 

D31 -.54 1000.00 

n22 -.12 592.86 

[11 0.00 285.71 

Dn41 .15 900.00 

ARCA 0.00 571.43 

ПМ „ei 250.00 

р12 0.00 857.14 
LOELF 15.33 942.86 


SAL кіш 942.857143 
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LINEAIRK-FROGRAMMERINGSFROGRAM ІМОРТІ 
WILT U ПЕ GEGEVENS VIA ТАРЫ INLEZEN?TYFE JA OF NEE °? 
CM LWUAti з 25077B» LOADER USED  41400BNEE 


GEEF ПЕ VOLGENDE SPECIFICATIE VAN HET FRORLEEM: 
AANTAL. BESLISSINGSVARIABELEN ? 9 


TOTAAL AANTAL RESTRICTIES 713 
AANTAL =< RESTRICTIES ? 11 
AANTAL. = RESTRICTIES ? 2 
AANTAL. += RESTRICTIES ? 0 


MOET ПЕ DOELSTELLINGSFUNKTIE БЕМАХІМЕЄКО OF GEMINIMEERD WORDEN 
ТҮРЕ MAX OF MIN ? MAX 


GEEF ШЕ 9 NAMEN VAN ПЕ BESLISSINGSVARIABELEN? 
D11 ; 


FOUTIEVE ІМУОЕК.ПЕ 9 NAMEN DIENEN TE BESTAAN UIT FRECIES VIER 
KARAKTERS(INCLUSIEF SFATIES) GESCHEIDEN DOOR KOMMA*S 
Dil ,012 »n21 +022 +031 ,П32 +041 »I42 v DM 
GEEF ПЕ BIJBEHORENDE COEFFICIENTEN VAN DE IIQELSTELLINGSFUNKTIE IN 
DEZELFDE VOLGORDE? 

1,4.9,2.4,6.6,4.5»8.5,2,1,7»,-80 
GEEF ACHTEREENVOLGENS VOOR ELK VAN ПЕ 13 RESTRICTIES EEN REGEL MET ПЕ 
NAAM VAN ПЕ RESTRICTIE GEVOLGD DOOR ПЕ TECHNISCHE COEFFICIENTEN ЕМ HET 
КЕСНТЕКІ II$ 
ПЕ 11 -z RESTRICTIES 

1?720N1/4»8»0»0»0»0»0»0»05»8000 

2?20N2/0»0»4»8»0»0»05»0»05»9350 

3?20N3/0»0»0»0»4»8»0»0»—10»1500 

4?20M4/0»0»0»0»0»0»4»8»0»6000 

S?MIIEM/0»0»0»0»0»0»0»0»15400 

Ф?”АКС1/-.75,.25,0»,0»,0»0,0»,0»0»0» 

7?TAÀRC2/0»0»—-.75»,259»0»0»0»0»0»0 

8?024ARC3/0»0»0»0»0»0»—.75».25»050 

9?АКС4/.25»-.75»,0,0»0»0»0»0»0»0 
10?АКС5/0»0».25»-.75,0,0»,0»,0»0,0 
117АКС6/0»0»0»0»,0»,0»,.25»-.275,0»0 


ПЕ. 2 = RESTRICTIES 
1?TDUE/1»1»1»1»1»1»1»1»-1» 5000 


2?ПМТҮ/1»-2,1»-2,1»-2,1»,-2»0,0 
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ECHTE VARTABELEN pii. 012 BOL .D22 "di (032:4:041 7 042 DM 


SLACK VARIABELEN 70мі ZON2 20М3 20М4 МПЕМ ARC1 ARC2 ARC3 ARCA 
ARCS ARCS 
ARTIFICIAL VARTABELEN TUWE ПМТҮ 


EINDOFLOSSING МА 14  ITERATIES 
TABLEAU ZONDER ПЕ KOLOMMEN ПЕК ARTIFICIAL VARIABELEN 


RASIS 
VARIA NIET-RASISUARIARBELEN 
BELEN 
20М1 АКС1 ZÜN2 20М3 32 20М4 
п21 . 56 1.14 227 .14 -2.00 .24 
ARC2 e Al 1.00 e 31 «13 “17% «2 
ARCS "22 „PE -.15 -.09 1.25 "Хх 
ARC A 0.00 0.00 0.00 0.00 0.00 .10 
МПЕМ 421 0.00 .21 .21 0.00 .21 
ARCA ۰07 .71 0.00 0.00 0.00 0.00 
п22 -.18 -,.97 -.07 - 07 1.00 -.12 
142 0.00 0.00 0.00 0.00 0.00 .05 
011 ۰04 “1.14 0.00 0.00 0.00 0.00 
П41 0.00 0.00 0.00 0.00 0.00 (15 
ПМ -.21 0.00 -.21 -.21 0.00 = 21 
D31 -.54 0.00 -.54 “429 2.00 -.54 
п12 e11 .97 0.00 0.00 0.00 0.00 
DOELE 14.97 ‚63 15.20 241 73 2,30 15,33 


VERVOLG TABLEAU 


BASIS | 
VARIA NIET-RASISUVARIARELEN 
BEL EN 
ARCO WAARDE BASIS 

21 -1.60 1058.93 

ARC2 -1.40 634,37 

ARCS 1.00 214.73 

ARCI 1.40 900.00 

MOEM 0.00 10.71 

ARCA 0.00 571.43 

n22 .80 939.29 

[142 ~. Во 500.00 

П41 1.60 900.00 

ПМ 0.00 389.29 

031 0.00 1348.21 

р12 0.00 857.14 


SAL 00: -8939.464286 
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Alternatief 
A B C D E 
ру] D11 ру] D11 ру] 
012 012 012 Ві; 012 
D21 021 021 021 021 
022 022 022 022 022 
D31 Рз] D31 D31 D31 
D41 Dun Dun D41 Dun 
Du? Du? 0.2 Б, В, 
ZON 1 ZON 2 DM DM DM 
MDEM MDEM MDEM MDEM MDEM 
ARC 1 ARC 1 ARC 1 ARC 2 ARC 2 
ARC 2 ARC 2 ARC 2 ARC 3 ARC 3 
ARC 3 ARC 3 ARC 3 ARC 4 ARC 4 
ARC 4 ARC 4 ARC 4 ARC 5 ARC 5 


- tabel 10.2. Basisvariabelen in de vijf alternatieve 
oplossingen - 


In tabel 10.2. zien wij, dat de verzameling basisvariabelen bij 
alternatief B dezelfde is als bij alternatief A. Alle schaduw- 
prijzen bij alternatief B zijn dus !") gelijk aan de schaduw- 
prijzen bij alternatief А. Een soortgelijke uitspraak kunnen 
мі) doen met betrekking tot de alternatieven р en E. 


Wij kunnen nu z als functie van y vrijwel geheel in een grafiek 
in beeld brengen (zie figuur 10.1). De alternatieven А, ві С; 

D еп Е уап tabel 10.1. komen overeen met de punten As В Є. р 

еп Е іп figuur 10.1. De hellingen van de door de punten A, B, 

С еп D getekende lijnstukjes zijn gelijk aan de met -1 vermenig- 
vuldigde waarden van u in tabel 10.1. 21). De punten С, D en 

E liggen nagenoeg ор één lijn (de schaduwprijzen van het natuur- 
gebied zijn in de alternatieven С, D en E nagenoeg gelijk). 
Teneinde de grafiek volledig te kunnen tekenen moeten wij nog 
het snijpunt van de lijn AB en de lijn DC berekenen. Wij kunnen 
dit onder andere doen door de benedengrens van de range of 
feasibility van het rechterlid van restrictie ZON 2 in alterna- 
tief B te berekenen, 


2 - welfare function 


fig. 10.1 


3000 y = aantal aren 
te sparen 
natuurgebied 


135 


Deze benedengrens is: 


mx CE, 21150, 2333,23, 205). 
ofwel 

-333,33 < Abu 
ofwel 

рю = 393,23. < Ma 
ofwel 

11166,67 « bj 


De coordinaten van punt S zijn dus 
y = 833,33, z - 18872 - 833,33 « 0,30 = 18622 


Met figuur 10.1. hebben wij de efficiency-curve van de beide 
doelstellingen gevonden. Alle punten op de lijn ASE zijn effi- 
ciente oplossingen. Tussen deze efficiënte oplossingen kan 
slechts op subjectieve gronden een keuze worden gemaakt. 


10.5. LINEAIRE PROGRAMMERING ALS HULPMIDDEL BĲ HET 
BESLUITVORMINGSPROCES TEN AANZIEN VAN DE 
RUIMTELIJKE ORDENING 


Het in paragraaf 10.2, gegeven voorbeeld is sterk gesimplifi- 
ceerd, De modellen, die in de praktijk bij het besluitvormings- 
proces ten aanzien van de ruimtelijke ordening worden gehanteerd 
zijn aanzienlijk groter. In de praktijk zullen in het algemeen 
naast het bouwen en slopen van woningen ook andere activiteiten 
in het model worden opgenomen. Daarbij moet dan worden gedacht 
aan het bouwen en/of slopen van voorzieningen (winkels, scholen 
en andere volgende dienstverlenende bedrijven en instellingen), 
aan het bouwen en/of slopen van stuwende dienstverlenende be- 
drijven en instellingen, aan het creëren van industrieterrein, 
openbaar groen op stadsniveau of stadsdeelniveau etc. Voorts 
zullen in het algemeen twee of meer perioden worden onderschei- 
den. Tenslotte zal er in het algemeen meer dan беп doelstelling 
zijn, die niet rechtstreeks in de doelstellingsfunctie tot uit- 
drukking kan worden gebracht (in het voorbeeld van paragraaf 
10.2. was er één doelstelling, die niet rechtstreeks іп de doel- 
stellingsfunctie tot uitdrukking kon worden gebracht namelijk 
het in standhouden van het natuurgebied in zone 2). Door deze 
toevoegingen wordt het model niet wezenlijk anders; wel wordt 
het model natuurlijk aanzienlijk groter. 
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Welke voor- en nadelen zijn verbonden aan het gebruik van line- 
aire programmeringsmodellen Бі) het besluitvormingsproces ten 
aanzien van de ruimtelijke ordening? 

Voordat wij deze vraag beantwoorden maken wij enkele algemene 

opmerkingen over het besluitvormingsproces ten aanzien van de 

ruimtelijke ordening. Enkele kenmerken daarvan zijn: 

- de besluitvorming vindt plaats op verschillende niveaus (Rijk, 
provincie, gemeente en soms nog een vierde niveau tussen de 
provincie en de gemeente, namelijk de agglomeratie of het 
gewest) 

- de besluitvorming op ieder niveau wordt gekenmerkt door het 
feit, dat er niet één "decision maker" is, die de doelstel- 
lingen vaststelt en het plan uiteindelijk vaststelt. Er zijn 
steeds meerdere groepen (bij een gemeente het college van 
burgemeester en wethouders, de gemeenteraad, enkele gemeente- 
lijke diensten zoals de afdelingen financien, openbare wer- 
ken en inspraakgroepen uit de bevolking), die aan de besluit- 
vorming meedoen | 

- de probleemstelling vraagt om een multidisciplinaire aanpak. 
Het door een gemeente geformeerd team voor het opstellen van 
een structuurplan kan bijvoorbeeld bestaan uit deskundigen op 
het gebied van architectuur en stedebouw, planologie, socio- 
logie, Ғіпапсіеп, verkeer en vervoer, natuur en landschap, 
milieu etc. 

- niet alle doelstellingen kunnen door middel van wegingsfacto- 
ren in dezelfde eenheid worden uitgedrukt (het instandhouden 
van natuurgebieden en het minimeren van reistijden zijn voor- 
beelden van twee doelstellingen, die niet in dezelfde eenheid 
kunnen worden uitgedrukt). 


Het gebruik van lineaire programmering als hulpmiddel bij het 
besluitvormingsproces ten aanzien van de ruimtelijke ordening 
heeft de volgende voor- en nadelen: 

- door gebruik te maken van een LP-model worden alle betrokkenen 
gedwongen zich te bezinnen op de doelstellingen. Indien geen 
gebruik wordt gemaakt van een mathematisch model, gaan de be- 
trokkenen al zeer snel denken aan bepaalde oplossingen in de 
vorm van kaartbeelden. Dit geeft dan aanleiding tot een onder- 
handelingssituatie, waarin sommige actoren de éne andere ac- 
toren een andere oplossing trachten te pousseren of te elimi- 
neren. Door gebruikmaking van een LP-model wordt de discussie 
allereerst gericht op de doelstellingen en pas daarna op de 
oplossingen. | | 

- het is mogelijk een groot aantal efficiente plannen op te 
stellen. Zo vermijdt men het anders bestaande gevaar, dat 
gediscussieerd wordt tussen plan A en plan B, terwijl plan C 
zowel plan A als plan B domineert. 
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- door middel van de gevoeligheidsanalyse van de coëfficiënten 
van de doelstellingsfunctie kan men nagaan welke graad van 
nauwkeurigheid voor bepaalde gegevens is vereist. Men kan bij 
het verzamelen van gegevens als volgt te werk gaan: inventa- 
riseer eerst welke gegevens direct beschikbaar zijn en maak 
ten aanzien van de ontbrekende gegevens schattingen. Voer deze 
gegevens inen pas de gevoeligheidsanalyse van de coëfficiënten 
van de doelstellingsfunctie toe. Besteed nu bij het onderzoek 
naar de nog ontbrekende gegevens de meeste aandacht aan die 
&egevens, waarvoor de oplossing het meest gevoelig is. 

Wij merken hierbij op dat de gevoeligheidsanalyse van een be- 
paalde coefficiënt van de doelstellingsfunctie geheel berust 
op de veronderstelling, dat alle overige gegevens ongewijzigd 
blijven (zie paragraaf 7.1.). Zijn van vele gegevens slechts 
grove schattingen bekend, dan hebben de conclusies die men 
uit de gevoeligheidsanalyse kan trekken dus slechts een zeer 
beperkte waarde. Daarbij komt nog, dat men bij de gevoelig- 
heidsanalyse van de coéfficiénten van de doelstellingsfunctie 
slechts één doelstelling in beschouwing neemt (de overige 
doelstellingen zijn dan in de vorm van restricties in het mo- 
del opgenomen). Ook hierdoor neemt de waarde van de conclusies, 
die men uit de gevoeligheidsanalyse zou kunnen trekken af. 

- men kan het LP-model gebruiken als eoördinatiemiddel in orga- 
nisatorische zin. Bij wijze van voorbeeld beschouwen wij het 
opstellen van een structuurplan voor een gemeente. Het opstel- 
len van het structuurplan kan men zien als een project, waar- 
aan verschillende personen, die dikwijls bij verschillende 
organisaties werkzaam zijn, deelnemen. Het opstellen van het 
LP-model kan in een dergelijke projectorganisatie de functie 
van coordinatiemiddel krijgen. Men verdeelt het project dan 
in fasen zoals vaststellen van doelstellingen, gegevens ver- 
zamelen eerste ronde, output eerste ronde bespreken, bijstel- 
len van doelstellingen, gegevens verzamelen tweede ronde, out- 
put tweede ronde bespreken etc. 

¬ men kan LP-modellen gebruiken als eoördinattemiddel in piano- 
logische zin. Om duidelijk te maken wat Wij daarmee bedoelen 
beschouwen wij het volgende voorbeeld. Stel men stelt een LP- 
model op voor een agglomeratie van 9 gemeenten. In het model 
worden 9 zones onderscheiden die overeenkomen met het grond- 
gebied van de 9 gemeenten. In het uiteindelijk gekozen plan 
vindt men dan onder andere de schaduwprijzen van de grond in 
ieder van de negen gemeenten. Deze schaduwprijzen kan men dan 
gebruiken als input-gegeven bij de LP-modellen, die men voor 
de afzonderlijke gemeenten kan opstellen. Op deze wijze komt 
een goede coordinatie tussen de betrokken gemeenten tot stand. 

- het gebruik van een LP-model kan heel goed worden gecombineerd 
met een goede inspraakprocedure door de bevolking. De inspraak 
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kan meer inhoud krijgen, doordat men met behulp van een LP- 
model gemakkelijk alternatieve plannen kan produceren. Deson- 
danks ervaren sommige bij de besluitvorming betrokken perso- 
nen het gebruik van een LP-model als "technocratisch". 


UITWERKING VAN ENKELE 
VRAAGSTUKKEN 


2-1. De eerste stap is het opschrijven van een tabel, waarin de 
voor dit probleem relevante data zijn vermeld: 


Capaciteitsbeslag Produkt Beschikbaar aantal 
in uren per stuk A B uren per week 
afdeling І 2 4 1800 
afdeling II 4 5 3000 


Van produkt А Кап men voor wat betreft afdeling I maximaal 
900 stuks, voor wat betreft afdeling II maximaal 750 stuks 
fabriceren. De bottleneck is dus afdeling II. Men kan maxi- 
maal 750 eenheden A fabriceren. 

Van afdeling I zijn dan nog beschikbaar 1800 - 750 ж 2 = 300 
uren. De behaalde dekking bedraagt 750 ж 200 = 7 150.000 ,-, 


Van produkt B kan men ten hoogste produceren 


min (SS , T = 450 eenheden. De bottleneck is dan afde- 


ling І. 
Van afdeling II zijn dan nog beschikbaar 3000 - 450 ж 5 = 750 
uren. De behaalde dekking bedraagt 450 ж 300 - 7 135.000,-. 


Het LP-model luidt: 
maximeer 


200x; + 300x, 


onder де voorwaarden 
2X| + 4X» « 1800 


4x] + 5x9 < 3000 


хі, X2 2 0 


X] 


II 
Het gebied van de feasible solutions is OABC. 
Wij tekenen in de grafiek lijn 01: 200х) + 300x, = 60.000. 
Door deze lijn zo ver mogelijk evenwijdig aan zichzelf van 
de oorsprong te verschuiven komen wij in punt B. Punt B komt 
dus overeen met de optimale oplossing. 

De coórdinaten van B vinden wij door oplossen van het stelsel 


I 


2xj + 4x2 = 1800 
аху + 5x2 = 3000 


; 7 500, x, = 200. De behaalde dekking is dan 


200 х 500 + 300 х 200 = 160.000 
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2-2. Пе gegevens, die voor dit probleem van belang zijn, staan 
vermeld in onderstaande tabel: 


 Capaciteitsbeslag staart- stoeltjes- Beschikbare capaciteit 
іп uren per stuk* klokken klokken in aantal uren per week* 
houtafdeling 3 1 160 
spuiterij 0,4 0,8 40 
montage 2 2 120 
bankwerkeri]j* 0 1 45 


* Voor de afdeling bankwerkerij is het capaciteitsbeslag en 
de beschikbare capaciteit niet uitgedrukt in uren maar in 
"aantal uurwerken". 


“2 


160 
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Het LP-model luidt: 
maximeer 


160х | + 200 х2 
onder de voorwaarden 


3x1 + X « 160 (1) 


0,4x, + 0,8x, = 40 (2) 
2x4 + 2х2 5 120 (3) 
X»? $ 55 (4) 

хі» X2 Š 0 


Het gebied уап de feasible solutions is OABCDE. De lijn 

160x, + 200x, = 8000 geeft de verzameling van oplossingen met 
een dekking van 8000. Deze lijn is in de figuur aangegeven als 
een onderbroken lijn. 

Het optimum is punt C. De cofodinaten van punt C vinden wij door 
de lijnen (2) en (3) met elkaar te snijden: 


0,4x4 + 0,8x9 = 40 
2x4 + 2x5 = 120 
Oplossen van dit stelsel levert x - 20, x - 40. De behaalde 
dekking is dan 160 * 20 + 200 * 40 - 11 200. 
2-3. Het gegeven LP-model is: 
maximeer 2x4, + 3x2 


onder de voorwaarden 


Х 1 + x) < 4 (1) 
-х + хә 2 1 (2) 
X15 х2 5-0 


Driehoek АВС is het gebied van 
de feasible solutions. Ver- 
schuift men de onderbroken 
lijn evenwijdig aan zichzelf 
zo ver mogelijk naar rechts- 
boven, dan ziet men dat punt 
C de optimale oplossing is. 
De coórdinaten van punt C 
zijn: x4 = 0, x2 = 4. De 
waarde van de doelstellings- 
functie in punt C is 12. 
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2-4. Het gegeven LP-model is: 
minimeer 2x4 + 3х2 


onder de voorwaarden 


ху t 2x5.» 8 (1) 
Х 1 + X9 > 6 (2) 
хі, хо 8 0 

х2 

6% 

4 


0 6 9 X1 


Het gebied van de feasible solutions is het gehele eerste kwa- 
drant met uitzondering van de vierhoek OABC. Dit gebied strekt 
zich dus uit tot het oneindige. Omdat we de doelstellingsfunctie 
moeten minimeren, moeten we lijn 2 zo ver mogelijk evenwijdig 
aan zichzelf naar de oorsprong toe verschuiven. Men ziet op die 
manier, dat punt B het optimum is. De coördinaten van B zijn: 

хр = 4, хә = 2. De waarde van de doelstellingsfunctie in punt B 
bedraagt 14. 


2-5. Het gegeven LP-model luidt: 
maximeer 10x, + 12 X» 


onder de voorwaarden 


Het gebied van de feasible solutions bestaat uit het lijnstuk AB. 
De optimale oplossing is punt В (x, = 20, хә = 5). De waarde van 
de doelstellingsfunctie in punt B bedraagt 260. | 


Opmerking: Restrictie (3) is niet van belang voor het bepalen van 
de grenzen van het gebied van de feasible solutions. Men noemt 
zo'n restrictie redundant (overbodig). 

2-6. Het gegeven LP-model is: 


maximeer 2х - X5 


onder de voorwaarden 


X1 wee x < 2 (1) 
X2 8-3 (2) 
X] + x 2 10 (33 


о 


X19 х2 ë 


(2) 


T 


4 10 


Het gebied уап de feasible solutions is het gearceerde gebied 
dat wordt begrensd door de х2-ав en de lijnen (3) en (1). Men 
moet de onderbroken lijn zover mogelijk naar rechts schuiven. 
Men kan in dit geval zo ver naar rechts schuiven als men maar 
wil, waardoor de waarde van de doelstellingsfunctie willekeurig 
groot kan worden. Het model heeft dus wel feasible solutions, 
maar niet een eindige optimale oplossing. Men noemt een dergelijk 
LP-probleem unbounded (onbegrensd). 

Als men in een praktijksituatie een dergelijk probleem ontmoet, 
dan betekent dat gewoonlijk dat er bij de modelformulering een 
fout is gemaakt. 
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Opmerking: Restrictie (2) is in dit probleem redundant. 


2-7. a. Het model luidt: 
maximeer 200 хі + 400 х2 


onder de voorwaarden 


EH + x, ك‎ 10 000 (1) 
3*1 + x9 5 7 000 (2) 
хі + 3x2 ك‎ 30 000 (3) 

кə 71 Ф Xx 2 1 000 (5) 


b. Vierhoek ABCD is het gebied van de feasible solutions. Punt C 
is het optimum. De coórdinaten van C vindt men door de lijnen 


(3) 


(5) 
7 000 (1) | 


(4) 


1 000 
X] 
м 13 233 
(2) еп (4) te snijden. Deze coördinaten zijn x4 = 4800, x2 = 


6400. De bruto-winst bedraagt in punt С 200 * 4800 + 400 ж 
6400 = 3.520.000. 


Opmerking: De restricties (1) en (3) zijn in dit probleem 
redundant. 


2-9, Hint: Noem х, de fractie van olie А en ха de fractie van 
olie B. Het heet bevat dan een fractie 1-хү-хә van өзге С. 


3-1. а. Na toevoeging van de verkooprestrictie luidt het LP-model: 
maximeer 200х) + 300 x2 


onder de voorwaarden: 


2x4 + 4хҙ 5 1800 (1) 
4x4 + 5x9 = 3000 (2) 
х | < 400 (3) 
Хі» X5 2 0 


Het gebied van de feasible solutions is het gearceerde 
gebied. Het optimum is punt D. De coórdinaten van punt D 
zijn хі = 400, za = 250. De waarde van de doelstellings- 
functie in punt D is 200 ж 400 + 300 ж 250 = 155 000. 


b. 


147 


450 


- 400 
< 
(2) (1) 


Om de schaduwprijs van de derde restrictie te berekenen ver- 
vangen wij deze door хү <S 401. Het optimum is dan het punt p 
dat ligt op de lijn 2X4 + 4x2 = 1800 èn op de lijn хі = 401. 
De coördinaten van D' zijn хү = 401, хә = 249,5. De waarde 
van de doelstellingsfunctie in punt D' bedraagt 200 х 401 + 
300 * 249,5 = 155 050. Doordat de verkooprestricties met één 
eenheid is verruimd is de waarde van de doelstellingsfunctie 
toegenomen met 155 050 - 150 000 = 50. De schaduwprijs van de 
verkooprestrictie is dus 7 50,-. 


Opmerking: Als de verkooprestrictie met één eenheid wordt ver- 
ruimd kan men één eenheid van produkt A extra verkopen. Men 
moet dan ook één eenheid van produkt A extra produceren; omdat 
afdeling I reeds vol bezet is, moet men van produkt B een hal- 
ve eenheid minder produceren. Per saldo neemt de behaalde dek- 
king toe met 200 - 1 ж 300 = 50 en dit is de schaduwprijs van 
de verkooprestrictie. 


De waarden, die de doelstellingsfunctie z aanneemt als men het 
optimum steeds bepaalt voor verschillende waarden van b3 zijn 
af te lezen uit onderstaande tabel. 


160 000 


я 


450 
449,5 
449 


200,5 
200 
200 


200 


2 


135 000 
135 050 
135 100 


159 950 
160 000 
160 000 


160 000 


helling 0 
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3*2. а, Het optimum is punt А 
š (7,0). De waarde van de 
2 doelstellingsfunctie in 
punt А is 14. 


b. Vervang restrictie II door x + 2x < 8. Het optimum wordt dan 
punt А" (8,0). De waarde van de doelstellingsfunctie in punt 
А" bedraagt 16. De schaduwprijs van restrictie II is dus f 2,-. 


helling 0 
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3-3. Het LP-model is: 
maximeer 200х1 + 300 x» 


onder de voorwaarden 


ECH ^ 4x2 « 1800 І 
4x4 + 5х2 2 3000 ІІ 
X15» x,20 
X? Voor c, = 200 wordt de richting 
van de lijnen met constante waar- 
den van de doelstellingsfunctie 
gegeven door (1). Het optimum is 
punt B (500,200). De waarde van 
(1) de doelstellingsfunctie z in punt 
B is 
ps z = 200*500 + 300*200 = 160 000. 


Als c4 daalt van 200 naar 160, draait de onderbroken lijn een 
eindje tegen de wijzers van de klok in om punt B en gaat over in 
lijn (2). Het optimum blijft dan punt B. Wel verandert de waarde 
van de doelstellingsfunctie in punt B; deze bedraagt nu 

z = 160 * 500 + 300 ж 200 = 140 000. 


Als сі verder daalt tot 120, draait de onderbroken lijn verder 
tegen de wijzers van de klok in tot lijn (3). Punt C is nu het 
optimum еп z = 200 х 0 + 300 х 450 = 135 000. 


Wij bepalen пи de grenzen, waarbinnen c, kan variëren zonder dat 
het optimum verschuift van punt B naar een ander punt. De lijn 
суху + 300x = constant is evenwijdig met lijn I als 

c4 : 2 = 300 : 4, dus als сі = 150. 
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De lijn суху + 300х) = constant is evenwijdig met lijn II als 
ср á = 300 : 5, dus als су = 240. 
De gezochte grenzen zijn dus 150 < c4 < 240. 


4-1. a. De restricties luiden als volgt: 


3xj + хә + 2x34 < 2000 manuren 
X] + 2X2 + 2x3 < 1400 machine-uren I 
2х] + 2X2 + хз < 1500 machine-uren ІІ 


Теп aanzien van de eerste restrictie merken wij nog het уо1- 
gende op. Het feit, dat uitbreiding of inkrimping van het 
aantal personeelsleden niet mogelijk is, betekent niet dat 
men de eerstgenoemde restrictie zou moeten schrijven als een 
gelijkheidsrestrictie. Men kan immers besluiten een bepaald 
aantal manuren onbenut te laten. Wel moet men de directe 
loonkosten als vaste kosten beschouwen. Dit houdt in, dat 
men de behaalde dekking I zal willen maximeren. De doelstel- 
lingsfunctie is dus 


50x; + 30х2 + 60x3 


en men wenst deze functie te maximeren. 


De oplossing is gegeven in onderstaande tabel: 
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b. In de optimale oplossing bedraagt de dekking voor vaste kos- 
ten en winst f 48.000,- per dag. De directe loonkosten Бе- 
dragen 250 ж 160 = f 40.000,- per dag, de kosten voor sala- 
rissen van niet direct produktief personeel, de afschrijvin- 
gen, de rente en de overige vaste kosten bedragen tesamen 
f 10.000,-. Het verlies van de onderneming bedraagt dus 
f 2000,- per dag. 

De schaduwprijs van arbeid vinden wij in de z-rij in de kolom 
van хц. Deze bedraagt dus f 10,- per direct manuur. Neemt 

men 1 man extra in dienst dan neemt dus de dekking voor di- 
recte loonkosten en overhead toe met f 80,- per dag. De di- 
recte loonkosten nemen toe met f 160,-, zodat het verlies 
toeneemt met f 80,- per dag. 

Zou men | man kunnen ontslaan, dan neemt het verlies af met 

f 80,- per dag. i 


4-72. 

a. 3; 7-400х1. * 300. = 0 
2x1 + Su k NA = 1800 
аха 7 OND + xy = 3000 


0 


1800 
3000 


135000 
450 
750 

160000 
200 
500 


b. Zij хз = het aantal eenheden, dat men van het produkt C wil 
gaan produceren. De spelingsvariabelen noemen wij nu x, еп 
X5. р 
Het LP-model is dan te schrijven als: 


„200: 9200 0 
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2 - 200x; = 300х; - 250x3 = 0 
2X] + 4х0 + 4х3 + x, - 1800 
Wy €; ка” Ок * Xs = 3000 


"ZUG NK 430 


і 
2 
n 
1 
3 


оДоодм 


OPMERKING 

Bekijken wij bovenstaande tabel, dan zien wij dat de eerste 
drie tableaus gelijk zijn aan de drie tableaus van het pro- 
bleem bij vraag a) met dien verstande dat er nu беп kolom 
bijgekomen is. 

Bekijken wij nu speciaal het derde tableau, dan zien wij dat 
de reduced revenue van produkt C in dit derde tableau gelijk 
is aan 50. Er geldt: 


50 = 250 = 4 ж 335 - 2 ж 331 


ofwel: de reduced revenue van produkt C in tableau 3 is 
gelijk aan de dekkingsbijdrage per eenheid van produkt C ver- 
minderd met de opportunity costs, die zouden ontstaan door 
Eén eenheid van dat produkt in de oplossing te introduceren. 
Deze opportunity costs zijn gelijk aan het aantal uren van 
afdeling I, die men zou moeten prijsgeven om één eenheid 

van produkt C te kunnen maken vermenigvuldigd met de scha- 


duwprijs per uur van afdeling І іп tableau 3 plus het aantal 
uren van afdeling II, die men zou moeten prijsgeven om één 
eenheid van produkt C te kunnen maken vermenigvuldigd met 

de schaduwprijs per uur van afdeling ІІ in tableau 3. 
Bekijken wij nu nog eens het laatste tableau van het pro- 
bleem van vraag a), dan zien wij dat dit tableau ook reeds 
de gegevens bevat, die wij nodig hebben om de reduced revenue 
van het nieuwe produkt C te kunnen berekenen. De enige gege- 
vens, die wij hiervoor nodig hebben zijn de schaduwprijzen 
per uur van de afdelingen I en II en de gegevens die over 
produkt C bij de formulering van vraag b) werden gegeven. 
Dàt produkt C in de optimale oplossing zou worden opgenomen, 
hadden wij dus ook kunnen afleiden zonder het probleem voor 
vraag b) opnieuw met de simplex-methode op te lossen. 

Wel moeten wij het probleem opnieuw nu voor 3 produkten met 
behulp van de simplex-methode oplossen om te bepalen hoevéél 
stuks van produkt C wij zullen moeten gaan maken en tevens 
om te bepalen welk produkt (A of B) nu de optimale oplossing 
verlaat. 


c. De reduced revenue van produkt B in het laatste tableau is 
- 7 25%. 
Door een eenheid уап produkt В іп de oplossing te introduce- 
ren neemt dus de behaalde dekking af met f 25,-. 


4-3. a. Het LP-model is: 


maximeer 160х4 + 200х2 


onder de voorwaarden 


3x, + хә 5 160 (1) ` houtafdeling 
0,4x4, + 0,8x, = 40 (2)  spuiterij 
2x, + ر×2‎ S 120 (3) montage 
хо 5 45 (4)  bankwerkerij 
хү, хо 2 Ü ` 


De oplossing is gegeven іп de tabel ор р. 153. 


b. Bij de houtafdeling hoort spelingsvariabele x3. In de optimale 
oplossing geldt x3 = 60. In de houtafdeling zijn dan 60 uren 
niet benut. | 
Bij de spuiterij hoort spelingsvariabele хл. De schaduwprijs 
van één arbeidsuur in de spuiterij bedraagt 7 100,-. 


Basis 

2 
*3 
XL 
^5 
Zë 

2 
-3 

2 X4 
55 

x2 

2 

"3 

3 X4 
43 

жә 

2 
3 

4 X1 
X6 
“2 


c. De nieuwe stoeltjesklok vraagt 2 uur arbeid van de houtafdeling, 
0,8 uur arbeid van de spuiterij en 2 uur van de montage. De 
reduced revenue van de nieuwe stoeltjesklok wordt dus 
240 - 2x0 - 0,84100 - 2460 = 40. 

Deze reduced revenue is positief. Het is dus aantrekkelijk om 


het nieuwe type stoeltjesklok in het produktieprogramma op te 
nemen. 


d. De oplossing is gegeven in de tabel op p.156. 


1 -160 -200 -240 Don MN es 
3 1 2 1 

1 
E 
PE 


x4 


4-4. Het LP-model is: 


maximeer 160х | + 200x» + 240x3 + 60х, 


onder de voorwaarden 


3x4 + х2 + 2X4 « 160 
0,4x4 + 0,4х2 + 0,8х: > 40 
2x4 + 2x, + 2X4 < 120 
хә t X3 + X, 2.439 


Хү, Хә, Хз» ХА e H 


De oplossing is gegeven іп de tabel op р.157. 


% 
FS 
"6 
24 
х8 
2 
^2 
2 X6 
mr 
ж 
2 
e 
3 


1 -160 -200 -240 -60 0 
DX 1: 2 а. 
0 3104-08: 0:8. 50 
0 0 
0 


oo о 


Opmerking. Deze oplossing is nagenoeg gelijk aan de oplossing van 
vraagstuk 4-3. In het optimum worden 20 staartklokken en 40 stoel- 
tjesklokken nieuwe stijl geproduceerd, net als bij de oplossing 
van vraagstuk 4-3. Voorts worden nog 5 uurwerken per week ver- 
kocht; in de oplossing van vraagstuk 4-3 werden alleen uurwerken 
voor eigen gebruik geproduceerd en was er in de bankwerkerij een 
onbenutte capaciteit van vijf uurwerken per week. 


b. De schaduwprijs van één arbeidsuur in de spuiterij bedraagt 
BeF PS, 
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4-5. Noem x het aantal hectaren, die ingezaaid worden met gewas 
А. Eén hectare ingezaaid met gewas A levert 2 ton verkoopbaar 
produkt А van gewas А. Eén hectare ingezaaid met gewas A geeft 

dus een opbrengst уап 2%800 = f 1600,-. De variabele kosten bedra- 
gen voor één hectare ingezaaid met gewas A f 100,-. De contribu- 
tiemarge bedraagt dus voor één hectare ingezaaid met gewas A 

f 1500,-. De contributiemarge voor één hectare ingezaaid met gewas 
B bedraagt 1,8 « 475 - 55 = 800. Een hectare ingezaaid met gewas C 
geeft een contributiemarge van 126 *.250 - 35 1165. 

De doelstellingsfunctie is dus 1500xj + 800х2 + 1165x353. 

Er zijn twee randvoorwaarden: 


хі, X) У 4 9 60  (grond, in hectare) 
30x. + 10x, + 15x3 5 1000 (arbeid, in uren) 
terwijl voorts moet gelden x4, x2, x3 2 0. 


De oplossing is gegeven in onderstaande tabel. 


Tableau 


-1500 -800 -1165 


1 
0 
0 


Opmerking. Een andere modelformulering verkrijgt men door xq te 
definiëren als het aantal ton van produkt A dat de boer wenst te 
verbouwen (en evenzo x4 en хә als het aantal tonnen van В respec- 
tievelijk C). Het model luidt dan: 


maximeer 750х, + 444,44х2 + 728, 12x3 
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onder de voorwaarden 


X x x 

1 2 3 

GU TIENES RE 
Fri + کل‎ хә + کا‎ x3 S 1000 


5-1. a. Wij beschouwen een hoeveelheid van 100 ml van het nieu- 
we haarwater. Stel deze 100 ml bestaat uit 
xjml A + xml B + хзпі C. Er geldt dus: 
| їх» ? 4X3 100 


Verder moet gelden 


X^ * X59 + 2х3 > 40 


X] * GA 60 


Het model luidt dus: 
minimeer 


0, Зх) + 0,2x» Ф 0,5хз 


onder de voorwaarden 


тес хе + хз = 100 І 
X1;*X» WAX; > An ІІ 
X] ? ҡас "60 III 


Х1» X29 X3 z 0 


Aan de tweede en derde restrictie voegen wij de surplusvari- 
abelen xi en xs toe. Aan alle drie de restricties voegen wij 
een kunstmatige variabele toe (xg, x; en xg). Het model luidt 


пи: 
minimeer 
0,3xj + 0,2х2 + 0,5х3 + Mxg + Mx, + Mxg 
onder de voorwaarden 
xi * X> + X3 * T - 100 
хі + КЖ” MES 7 Xš - жо = 40 


X] * X3 "e ae + хӯ = 60 


160 А 


De oplossing is gegeven іп onderstaande tabel: 


b. De schaduwprijs van restrictie III vinden wij in de z-rij in 
de kolom van xs. Deze bedraagt 0,1 cent per mg. Door te eisen 
dat 100 ml van het nieuwe boorwater 61 mg anti-haaruitvalmid- 
del moet bevatten, stijgen de grondstofkosten van 26 cent per 
100 ml tot 26,1 cent per 100 ml. 


Opmerking. Men kan het onder a. geformuleerde model vereenvoudi- 
gen door хз = 100 - x, - x, te substitueren in de restricties II 


еп ІІІ en іп de doelstellingsfunctie. De doelstellingsfunctie 
wordt dan: 
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minimeer 0,3x, + 0,2х2 + 0,5(100-x,-x5) 
en dit kan men ook schrijven als 
maximeer 0,2x + 0,3х2 


De restricties luiden dan 


хі + X4 5 160 ІЗ 
X9 5 40 ІІІ 
Х4» хә 2 0 


terwijl de restrictie хз 2 0 overgaat іп 100=x;=xə š 0, 
dus x1 + x2 S 100. Men ziet nu direct dat restrictie II overbodig 
15. Het vereenvoudigde model luidt dus: 


maximeer 0,2x, + 0,3x> 
onder de voorwaarden 
хү + хо 5 100 
ко 5. 40 


5-2. Het model luidt: 
minimeer 4x] + x2 + хз + 4X4 
onder de voorwaarden 
4x1 + 2X4 + x, 2 50 
2X4 *: Xo! $. X3 > 40 
2%; + 2x; š 10 
X], Хә, X4, X, 2 0 


Na toevoeging van surplusvariabelen en kunstmatige variabelen 
luidt het probleem: 


minimeer 4x4 + хә + X3 * Mxg * Mxg * Мх10 
onder de voorwaarden 
4x1 T UA Ki = Xs + ха = 0 
24 Kg + хз 7X6 % хо = 
2X4 + 2х, “ka "хау 


De oplossing is gegeven in de tabel ор р.162. 


5-3. Laat men de restricties (1) еп (3) weg, dan luidt het model: 
maximeer 200x, + 400x, 


onder de voorwaarden: 


Ze + хә < 7 000 
3x4 + 4x2 S 40 000 
-іхі + хә 2 100 


X45 xj 20 
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Na toevoeging уап spelingsvariabelen en een surplusvariabele en 
een kunstmatige variabele luidt het model: 


maximeer 200х | + 400x» -Mxg 


onder de voorwaarden 


| + хә + x3 - 7 000 
Зх | ^ 4х, + X4 - 40 000 
-іхіз ха > Ka е kaja: 1:000 


De oplossing is gegeven in onderstaande tabel: 


-200 -400 
0 


6/5 -1/20 


164 


6-1. Het duale probleem is: 
minimeer 18001 + 3002 


onder de voorwaarden 


2u4 + hus 2 200 (A) 
hug + 502 > 300 (В) 
04» 02 2-0 


De duale variabelen uy respectievelijk u2 Кап men opvatten als de 
prijzen, die een buitenstaander zou willen betalen voor het huren 
van één arbeidsuur van afdeling I respectievelijk afdeling II. 


Restrictie (A) kan men als volgt interpreteren: als de prijzen 

u, еп ua zodanig zijn, dat voldaan is aan restrictie (A), dan is 
het voor de oorspronkelijke ondernemer niet meer aantrekkelijk om 
produkt A te fabriceren. Een soortgelijke interpretatie geldt 
voor restrictie (B). 

De doelstellingsfunctie geeft het totale bedrag, dat de buiten- 
staander moet betalen voor het huren van alle produktiecapaciteit. 
De buitenstaander wenst dit bedrag te minimeren. 


Voegt met aan bovenstaand model surplusvariabelen en kunstmatige 
variabelen toe, dan luidt het: 


minimeer 1800u, + 300 u> + Mus + Mug 
onder de voorwaarden 

241 + 609 " u$ + 05 = 200 

дар + 5u2 ў + ug 72300 


De oplossing is gegeven іп de tabel op р. 165. 


6-2. Het oorspronkelijke probleem is: 
minimeer 3x1 + 5х2 

onder de voorwaarden 

6 

8 


ІМ 


хі 7..х2 


ІМ 


X1 + 2х2 


X15» x) е 


о 


In primale vorm geschreven is dit: 


max imeer "E Т 5X2 


-1800 -3000 
0 0 
2 


-500 


«2/3 ЫЗ 
0 1 0 5/6 .—243 4 


-200 


"x; з X5 Š -6 

"Zi » 2х, 2 -8 

X1, X2 2 0 
Het duale probleem luidt: 
minimeer "60; ~ 8и, 


onder de voorwaarden 


c 

c 
N 

І 

о 


Dit kan ook worden geschreven als 


maximeer би! + 8u, 
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onder de voorwaarden 
di * 42 5 3 
3 


IA 


Ui + 209 


о 


ЧІ» u2 É 


De oplossing is gegeven in onderstaande tabel: 


2 
3 
WA 


6-3. Het duale probleem is: 


maximeer 50u, + 40u, + 10u3 


onder de voorwaarden 


hu, + 20, < A 
u2 $1 

20} + ua + Zu: 5 1 
u1 + Zu: < 4 


ці» Ugs UZ Z 0 


De oplossing is gegeven іп de tabel ор р.167. 
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2 

u4 
u5 
u6 
u7 


6-4. a. De primale vorm is: 
maximeer -10х; - láx, - хз 
onder de voorwaarden 
TX] + х2- 2x4 5-16 
-5x4 + 43.7 3X4 S 142 


хі + X5 + 4X4 S.. 240 
OE MR 4X4 5 vn 
X1» X9, X420 
b. Het duale probleem is: 
minimeer -16u, - 12u, + 20u3 = 20u, 


onder de voorwaarden 


ІМ 


-10 
-14 
w D 


-u4 ki 502 + ug us HA 


ІМ 


ug * чо T ug Ron Uu, 


ІМ 


-2u, Те 3и, d hug be 4u, 


| 
о 


ЧІ» u»; Ugs u, ë 
Men kan dit probleem ook schrijven als: 


maximeer 16u, + 120, - 20u4 * 20u, 


onder de voorwaarden 


чі + 5ч2 е9 ko us 510 
MEA > NN + "uq S Fé 
244 Ф 3u5 T" hug + бл, = 5 


Чі» Uy, وولا‎ HA zd 


7-1. а. De schaduwprijs уап één manuur bedraagt f 10,-. Door één 
man te ontslaan zou men het verlies kunnen reduceren tot 7:80)" 
per dag. Het totale verlies bedraagt f 2000,- per dag. Door 
25 man te ontslaan kan men dus inderdaad het verlies tot nul 
reduceren, mits de range waarbinnen de schaduwprijs van arbeid 
gelijk blijft aan f 10,- per uur voldoende groot is. 

Deze range wordt bepaald door: 


max 7) < Ab, < nin» = 


ofwel 
-600 < Ab; < 4665 
ofwel 
1400 < bj < 24665 
Door het ontslag уап 25 man zou het aantal directe manuren 
dalen tot 1800. Dit ligt binnen de boven berekende range, 


zodat men mag concluderen, dat de schaduwprijs van arbeid 
f 10,- blijft als men 25 man zou ontslaan. 


b. De range van cj wordt bepaald door 
nax (j) « Aci < nin, = 
ofwel | 
20 * Ae < 40 
ofwel 


30 « сүт 90 
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Aangezien с, - 40 binnen het interval 30 < су < 90 ligt, ver- 
andert de optimale oplossing niet. Wel neemt de behaalde дек- 
king af met 300 х 10 = f 3000,-. Het verlies neemt hierdoor 
toe van f 2000,- per dag tot f 5000,- per dag. 


c. De reduced revenue van produkt B in het laatste simplex- 
tableau is -f 20,-. Pas als de dekkingsbijdrage van produkt B 
groter is dan f 50,- per eenheid, wordt produkt B in het opti- 
male produktiepakket opgenomen. De verkoopprijs van produkt B 
zou dus tenminste f 120,- moeten bedragen. 


7-2. а. De reduced revenue van produkt B is -f 25,-. Als de dek- 
kingsbijdrage van produkt B toeneemt met f 30,-, wordt produkt 
B dus in het optimale produktiepakket opgenomen. 


b. De range van C3 wordt bepaald door 


max, 229) < Ac3 < min (>) 
ofwel 

«30 < Ас; < 190 
ofwel 


200 « c4 < 400 


De nieuwe waarde с3 = 225 ligt binnen bovenstaand interval. 

De samenstelling van het optimale produktieprogramma verandert 
dus niet als de dekkingsbijdrage van produkt C daalt van 

LSU" et f 225,-. 

De in totaal behaalde dekking daalt daardoor wel, en wel met 

100 х 25 = f 2500,-. Als de dekkingsbijdrage van produkt C 

f 225,- is, bedraagt de in totaal behaalde dekking dus 

7 14.000,-. ' 


c. De schaduwprijs van een arbeidsuur in afdeling II bedraagt 
f 25,-. Om te bepalen of de schaduwprijs constant blijft als 
de capaciteit van afdeling II toeneemt van 3000 tot 3500 uur, 
berekenen wij de range of feasibility: 


-700 7100 
4732. 7925 2176 


“2500 < Ab» « 600 
ofwel 
900 < b» « 3600 


Als b toeneemt van 3000 tot 3500 blijft de schaduwprijs dus 
constant; de in totaal behaalde dekkingsbijdrage neemt door de 


vergroting van de capaciteit van afdeling II toe met 500 ж 25 
u$ 12:5090,-. 

Als bj toeneemt van 3000 tot 4000 kan men met de gegevens uit 
het in de opgave gegeven laatste simplex-tableau niet berekenen 
met welk bedrag de in totaal behaalde dekkingsbijdrage toeneemt. 


8-1. a. Het LP-model is: 
maximeer  1500x, + 800x» + 1165x3 


onder de voorwaarden 


IA 


60 
Ы; 


Zi + x % Жа 


20x, + 10x» + 15X4 


IA 


Dit probleem moeten wij oplossen voor verschillende waarden 
van bj, die voldoen аап 0 < b, < 1500. 

Voor bas 1000 is dit probleem al opgelost in vraagstuk 4-5 
Er geldt dan z = 76 600, terwijl de schaduwprijs van één 
arbeidsuur f 67,- bedraagt. De range of feasibility wordt 
bepaald door 


-20 -40 
175 < 452 < “175 
ofwel -100 < Ab; < 200 ofwel 900 < b2 < 1200. 


Wij kiezen nu een waarde van b2 groter dan 1200, maar kleiner 
dan 1500, bijvoorbeeld bə = 1300. Door het probleem opnieuw op 
te lossen vinden wij z = 90 000, de schaduwprijs van een uur 
arbeid is nul en de range of feasibility is b2 > 1200. 

Wij kiezen nu een waarde voor b2, die voldoet aan 0 < b2 < 900, 
bijvoorbeeld b2 = 800. Bij b2 = 800 vinden wij z = 62 600, de 
schaduwprijs voor één uur arbeid is f 73,- еп de range of 
feasibility is 600 < b2 < 900. 

Wij kiezen nu een waarde voor b2, die voldoet aan 0 < b» « 600, 
bijvoorbeeld b2 = 500. Bij b2 = 500 vinden wij z = 40 000, de 
schaduwprijs van één uur arbeid is f 80,- en de range of 
feasibility is 0 < b» < 600. 


Samenvattend hebben wij nu gevonden: 


БАТЫ schaduwprijs | range of feasibility 


0 < bz < 600 
600 < ba < 900 
900 < bz < 1200 
b» » 1200 
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Om de grafiek te kunnen tekenen bepalen мі) de waarde van 2 
voor die waarden van by, die overeenkomen met de grenzen van 

de ranges of feasibility, dus voor b2 = 0, 600, 900 en 1200. 
Voor b» = 0 geldt z = 0. 

Voor b2 = 600 kunnen wij de waarde van z op twee manieren bere- 
kenen, namelijk ten eerste door uit te gaan van de waarde van 

2 voor bas 500 en ten tweede door uit te gaan van de waarde 
van z voor bj = 800 (dit geeft dus een controle op mogelijke 
rekenfouten). Wij vinden dan 


z = 40 000 + 100 х 80 48 000 
еп 2 - 62 600 - 200 х 73 48 000. 


Na soortgelijke berekeningen voor b2 = 900 en bə = 1200, hebben 
wij gevonden: 


Aangezien geldt bj = 1500 - y hebben wij nu ook de functie 
2(у) gevonden. Deze is geschetst in onderstaande figuur. 


2 


helling 
is 0 


90 000 


69 900 


helling is -73 


48 000 


helling is -80 


900 1500 aantal uren 
vrije tijd 


300 600 


b. De verzameling van efficiënte oplossingen is het gebroken lijn- 
stuk BCDE. 
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. a. Bij het vertalen van een praktijkprobleem in een LP-model 
is het meestal verstandig te beginnen met het aangeven van de 
beslissingsvariabelen, die men in het model wenst op te nemen. 
Daarbij zijn dikwijls verschillende mogelijkheden. 
In dit geval kan men als beslissingsvariabelen kiezen: 
C; = het aantal in periodëé i te produceren caravans 
I; = de toeneming van het aantal werknemers aan het begin 
van periode i 
K; = de afneming van het aantal werknemers aan het begin 
van periode i. 
Heeft men eenmaal vastgesteld hoeveel caravans men in periode 
i wil produceren, dan ligt de voorraad aan het einde van perio- 
de i vast. De voorraad aan het einde van periode i behoeft men 
dus niet op te nemen als beslissingsvariabele. Hetzelfde geldt 
voor "het aantal werknemers gedurende periode 1". 


Wij formuleren nu eerst de doelstellingsfunctie. 


Loonkosten: Het aantal werknemers gedurende de drie perioden 
bedraagt: 


perioge 1 : 200 + kr оку 
“periode 2 : 200 Я L, ЖЖ, $ Ig," Ко 
periode 3 : 200 + Яр « КУ WIET Бо * I4 — K4 
De loonkosten bedragen dus in totaal: 
10 000(200%1,-К|) + 10 000(200-1,-К|%12-К2) + 
+ 10 000(200%1,-К|%12-К2%1;-Қ4) 
Dit kan men herschrijven als 


10 000(314 + 21) + Iş = 3K, - 2К; - Кз) + 6 000 000 


Voorraadkosten: De voorraad caravans bedraagt aan het einde van 
periode 1 : 200 + C, - 600 
periode 2 : 200 $ Ç I = 600 + C5 ^ 800 
periode 3 : 200 + C, - 600 + С» - 800 + C3 - 2000 
De voorraadkosten bedragen: 
180 (200+C4-600) + 180 (200+C4-600+C5-800) + 
+ 180 (200+C4-600+C2-800+C3-2000) 
Dit kan men herschrijven als 


180(3C, + 2C, + C4) - 864 000 
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Afvloeiingskosten еп wervingskosten. Deze bedragen 
нон ни Wervingskosten 


15 000(K, + K> Ф Кз) ` ea 000 (I, F 12 ғ 13) 


Пе doelstellingsfunctie luidt dus (de constante term laten wij 
weg): 


minimeer 540C4 + 360C, + 180C3 + 31 000 I, + 21 000 1» + 
4 11 000 Із 15 000 КІ - 5 000 K, + 5 000 K3 
Wij formuleren nu de restricties. Er zijn twee typen restric- 


ties, namelijk capaciteitsrestricties en vraagrestricties. 
De capaciteitsrestricties zijn: 


САР! : 80C, < (200-1,-К|).580 
CAP3 : 80C3 < (200*1,-K,*I5-K,*I4-K4).500 


De vraagrestricties luiden: 
РЕСІ : 200+C4-600 > 0 
DEG Ze 200*C4-6004C5-800 > 0 
DEC3 : 200+C4-600+C-800+C3-2000 > 200 


Na herschrijven luiden de restricties: 
САРІ : 80C, - 5801, + 580K, = 116 000 
CAP2 : 80C, - 5801, - 5801, + 580K, + 580K, < 116 000 
CAP3 : 80С3 - 5801) - 5801, - 5001; + 580K, + 580K, + 
+ 500K4 < 10 000 
РЕСІ : Су 2 400 
DEC2 : C4 + C 2 1200 
DEC3 : C4 + C + C4 2 3400 


IV 


. а. Wij kiezen als beslissingsvariabelen: 

Сі = het aantal in periode і te produceren caravans 

Dj = het aantal in periode i te produceren direktieketen 

Aj = het aantal werknemers in éénploegendienst in periode i 
Bi = het aantal werknemers in tweeploegendienst in periode i 
Ili = de toeneming van het aantal werknemers aan het begin van 


periode i 
Кі = de afneming van het aantal werknemers aan het begin van 
periode і 


Oi = het aantal gewerkte overuren in periode i 


. De loonkosten bedragen inclusief ploegentoeslag: 


De 


De 
De 


De 


mi 


Soo — = 


. De 


10 000(А|%А2%А4) + 12 500(B,*B5*B3) 

kosten уап ontslag en van werving уап personeel ton: 
15 000(K4*K2*K3) + 1000(I4*I5*13) 

overwerkkosten bedragen: 30(01%02%03) 

voorraadkosten bedragen: 


540С + 360С2 + 180C3 + 48004 + 32002 + 160D3 — een constante 


doelstellingsfunctie is: 

nimeer 
10 000(А|%А2%Аҙ) + 12 500(B,*B5*B4) + 15 O00(K,*K5*K4) + 
+ 1000(I4+19+13) + 30(0,*05*04) + 540C, + 360С; + 180C3 + 
+ 4800, + 320D, + 1600: 


bain: 
capaciteitsrestricties: de produktiecapaciteit wordt begrensd 
door het aantal beschikbare arbeidsuren. Noem deze capaci- 
teitsrestricties САРІ, САР2 еп CAP3. 
vraagrestricties: men moet tenminste zoveel caravans en 
direktieketen produceren dat aan de vraag kan worden voldaan. 
Noem deze restricties DEC1, DECH, DECH, DED1, DED2 en DED3. 
overwerkrestricties: het aantal overgewerkte uren per регіо- 
de mag niet groter zijn dan 20% van de in die periode in de 
normale werktijd beschikbare uren. Noem deze restricties 
ОУЕ1, OVE2 еп OVE3. 
ploegendienstrestricties: ten hoogste 1/3 deel van het aantal 
werknemers is bereid in tweeploegendienst te werken. Noem 
deze restricties РІ,01, PLO2 еп PLO3. 
ruimterestricties: in de fabriek kunnen ten hoogste 300 per- 
soneelsleden tegelijk werken. Noem deze restricties RUI!, 
RUI2 en RUI3. 
definitierestricties, die A;, Bi, І; еп K; aan elkaar rela- 
teren. Noem deze restricties УАМ1, VAW2 en VAW3. 


restricties luiden als volgt: 


САР! : 80C, + 60D, = 580A, + 580B, + 0; 
CAP2 : 80С» + 6002 < 580Ag + 5808, + 0» 
CAP3 : 80C3 Lë 60D4 5 500A4 + 50083 * 03 


DECI : 200 + C, = 600 > 0 
DEC2 : 200 + C, - 600 + С; - 800 > 0 
DEC3 : 200 + C, - 600 + Cj - 800 + C} - 2000 > 200 


IV 
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V 


DED1 : 100 + Dr - 725 2 0 
DED2 : 100 + D - 725 + р; - 750 2 0 
DED3 : 100 + D4 —25 4 D» - 250 A D3 - 800 > 100 


ІМ 


OVE! 3 04 Š 116A, 
OVE2 : 02 116A» 
OVE3 : 04 S 10043 


IA 


PLOI : Ву = 1/3(A4*B4) 
PLO2 : Вә < 1/3(A5*B3) 
PLO3 : Вз 5 1/3 (АзтВз) 


RUI| : Ар + ÀB, = 300 
RUIZ : A, + 1B, = 300 
RUI3 : Аҙ + ]B4 = 300 
VAI 1 EBS 267 T1: Қ. 
VAW2 : A9 + В» = 
VAW3 : A3 + Вз = Aj + B) + I} - Кз 
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6) 


7) 


VOETNOTEN 


Voor een goed overzicht leze men Ackoff en Sasieni, hoofd- 
stuk | 


Zie bijvoorbeeld Kramer en de Smit, paragraaf 6.4. en de 
Leeuw paragraaf 5.4. 


De іп paragraaf 1.4. beschreven modelcyclus vertoont ook ge- 
lijkenis met de in de literatuur van het organisatie-onder- 
zoek voorkomende besluitvormingscyclus. Zie bijvoorbeeld 
Botter, paragraaf 3.3.2. 


Voor een goed overzicht van verschillende soorten modellen 
zie men bijvoorbeeld Bertels en Nauta. 


De zin en de gebruiksmogelijkheden van dergelijke modellen 
zijn beschreven in Douma, 19/3, Een soortgelijk model voor 
de middellange termijnplanning bij een gemeente is beschre- 
ven in Douma, S.W. еп F.J. van de Linde. 


Voor de lezer die het begrip opportunity cost niet kent ge- 
ven wij de volgende toelichting. De opportunity cost van het 
aanwenden van een bepaald produktiemiddel voor een bepaalde 
activiteit wordt gevormd door het verloren gaan van opbreng- 
sten, die verkregen hadden kunnen worden door het aanwenden 
van dat produktiemiddel voor andere activiteiten. Indien een 
tandarts op een werkdag gaat vissen, bestaat de opportunity 
cost van het aanwenden van zijn arbeid voor de activiteit 
vissen uit de inkomsten, die hij had kunnen verkrijgen door 
zijn arbeid voor een andere activiteit (het behandelen van 
patienten) aan te wenden. 


Wij bewijzen dit als volgt: Stel er is een oplossing anders 
dan xj = 40, х) = 40, хз = 0, х, = 20, х5 = 0 die een hogere 
waarde van z heeft. Iedere feasible solution moet voldoen 
aan het stelsel vergelijkingen in tabel 4.5. (immers dit 
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stelsel is equivalent met het oorspronkelijke stelsel). Stel 
deze andere oplossing is (x, = хі, Kä xl, de хі, ж" 
y. 7 xi). Nu moet tenminste ofwel x ^ 0 ofwel Ж. 70 zijn 
(immers als xl - 0 en xt - 0, dan is хі - 40 еп e - 40 еп 
*: = 20 zoals dan uit de tweede, derde en vierde rij van het 
tableau van tabel 4.5. gemakkelijk is af te leiden, еп dan 
hebben wij dus dezelfde oplossing). Maar als ofwel x} > 0 
ofwel хі > 0 (of zelfs beiden 0) dan Кап 

= 3000 - 15х1 e EE nooit groter dan 3000 zijn. 


8) Deze keuze is niet erg fundamenteel. Wij kunnen ook een ап- 
dere variabele kiezen, mits het element van de gekozen vari- 
abele in de z-rij negatief is. De ervaring leert echter, dat 
het aantal iteraties dat nodig is om de optimale oplossing 
te vinden in het algemeen kleiner is wanneer men de door ons 
geformuleerde regel hanteert dan wanneer men willekeurig 
kiest. 


9) Wij kunnen ook schrijven: 
maximeer 
"IR 905 
onder de voorwaarden 
AT Фо хо 6 
ki; * 2849 2 3 
X], X2. 2 0 
De vorm waarin wij de doelstellingsfunctie schrijven (als 
maximalisatie of minimalisatie-probleem) doet niet ter zake. 


Wel zullen wij de restricties steeds zo schrijven, dat de 
constanten in het rechterlid alle > 0 zijn. 


10) Een andere formulering van het model verkrijgt men door de 
restricties LAB І en OVE І te vervangen door: 
0,05А) + 0,08В) - O, < 350 
en 
01 < 175 


en de restricties LAB 2, LAB 3, OVE 2 en OVE 3 door soort- 
gelijke restricties te vervangen. De beslissingsvariabelen 


11) 


12) 


13) 


U}, U, en Оз zijn іп dat geval overbodig. Het model wordt 
daardoor iets kleiner. 

De in paragraaf 9.2, gegeven modelformulering kan ook worden 
gebruikt, indien ook aan het niet-gebruiken van beschikbare di” 
rect produktieve manuren kosten verbonden zijn. In dat geval 
krijgen de variabelen Оу, 0» en Оз een positieve coefficient 
in de doelstellingsfunctie. Bij de kosten van het niet-ge- 
bruiken van beschikbare direct produktieve manuren moet men 
niet denken aan de directe loonkosten. Het aantal werknemers 
is in deze probleemstelling immers een gegeven. De directe 
loonkosten worden dus beschouwd als vaste kosten en komen 
derhalve in het model niet voor. Bij de kosten van het niet- 
gebruiken van beschikbare direct productieve manuren moet 
veeleer worden gedacht aan het feit, dat onderbezetting de- 
motiverend kan werken en als zodanig "psychische kosten" kan 
veroorzaken. 


Aan de breuken met in de teller de coëfficiënten van de niet- 
basisvariabelen in de rij van de doelstellingsfunctie en in 
de noemer de coëfficiënten van de niet-basisvariabelen in de 
03-гіі voegen wij nu in afwijking van de in paragraaf 7.1. 
beschreven procedure geen minteken toe. De reden is, dat ons 
probleem een minimeringsprobleem is, terwijl in paragraaf 7.1. 
een maximeringsprobleem werd besproken. 


De methode die wij in hoofdstuk 10 beschrijven is ontwikkeld 
door het Israel Institute of Urban Studies te Tel-Aviv. In 
Nederland is deze methode geïntroduceerd door Raadgevend 
Bureau Berenschot. De methode staat in Nederland bekend on- 
der de naam Urbanics. 


In de "welfare-coefficients" van een bepaalde activiteit komt 
de bijdrage van die activiteit tot de welfare-function tot 
uitdrukking. Bij het definieren van een dergelijke collec- 
tieve welfare-function doen zich talrijke problemen voor. 

Een van deze problemen wordt gevormd door het feit, dat men 
het nut van een bepaalde oplossing voor een bepaald individu 
niet kan vergelijken met het nut van diezelfde oplossing 

voor een ander individu. In de welvaartseconomie gebruikt 
men daarom het begrip Pareto-optimum. Een oplossing is Pare- 
to- optimaal, indien er geen andere oplossing bestaat die voor 
tenminste één individu beter en voor alle andere individuen 
niet slechter is. 

In de in dit voorbeeld gehanteerde welfare-coefficients komt 
tot uitdrukking de marktwaarde van de te bouwen woningen ver- 
minderd met de bouwkosten en de grondkosten en eveneens ver- 
minderd met de gekapitaliseerde transportkosten. Bij de 
grondkosten moet worden gedacht aan de waarde van de aanwen- 


14) 


15) 
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dingsmogelijkheden van de grond, die door het bouwen уап 
woningen verloren gaan, niet aan de marktwaarde als bouwgrond. 
In de transportkosten is onder andere de waarde van de reis- 
tijden begrepen. Men zie voor de specificatie van de doel- 
stellingsfunctie ook Ben-Shahar. 


Voor de motivering van het woordje "dus" verwijzen wij naar 
de eerste drie alinea's van paragraaf 7.2. 


In paragraaf 3.2. hebben wij bij de behandeling van het pa- 
rametrisch programmeren z geschetst als functie van een der 
parameters in het rechterlid. Het ging in paragraaf 3.2. net 
als hier om een maximeringsprobleem. Beschouwen wij figuur 
3.3. dan constateren wij dat 2 = f(bj) monotoon stijgend en 
convex (bol) is. Men kan bewijzen, dat bij een maximerings- 
probleem z als functie van één der parameters in het rechter- 
lid altijd monotoon stijgend en convex is. De eigenschap dat 
2 monotoon stijgend is betekent, dat bij een maximeringspro- 
bleem een toeneming van de capaciteit van één der restricties 
nooit kan leiden tot een afneming van z. De eigenschap dat 

Z convex is, betekent dat de marginale opbrengst van een 
vergroting van de capaciteit nooit kan toenemen bij vergro- 
ting van de capaciteit. 

In figuur 10.1. schetsen wij z niet als functie van één der 
parameters in het rechterlid, maar als functie van y. Omdat 
in het rechterlid van restrictie ZON 2 voorkomt het getal 
12000-y, is z als functie van y monotoon afnemend en convex. 
Aan de waarde van u2 in tabel 10.1. moeten wij om dezelfde 
reden een -teken toevoegen om de helling van de lijnstukjes 
in figuur 10.1. te vinden. 
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Over het boek 


Lineaire programmering is een bekende optimaliseringstechniek. 
Voor het gebruik van lineaire programmering in de praktijk 
waren tot voor kort grote  computerconfiguraties nodig. 
Tegenwoordig kan men voor niet te grote LP-problemen ook 
gebruik maken van een mierocomputer. Daarmee is het gebruik 
van LP binnen het bereik van zeer velen gekomen. 

In dit boek ligt het accent op een intuitieve, niet-wiskundige 
benadering. Veel aandacht wordt gegeven aan de economische 
betekenis van de bij lineaire programmering gebruikte 
begrippen en aan de betekenis van deze begrippen voor de 
besluitvorming. 

De belangrijkste leerdoelen van dit boek zijn de student te 
leren een probleem te formuleren als een lineair 
programmeringsprobleem, dat probleem vervolgens op te lossen 
met behulp van de (тісго)сотриїег еп de verkregen resultaten 
te interpreteren. 


